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I . INTRODUCTION 


Although  nonparametr ic  procedures  are  powerful 
tools  to  the  analyst,  they  are  currently  underused  and 
often  avoided  by  potential  users-  Perhaps  one  reason 
for  this  is  the  difficulty  in  generating  the  exact 
distributions  of  the  test  statistics,  even  for  moderate 
sample  sizes.  Consequent ly ,  tables  of  these 
distributions  are  only  available  for  very  small  sample 
sizes  and  normal  theory  based  approximations  must  then 
be  used. 

The  purpose  of  this  thesis  is  to  make  a  variety  of 
nonparametr ic  procedures  quick,  easy  and  accurate  to 
apply  using  menu  driven  computer  programs  in  APL.1 
These  programs  use  enumeration,  recursion,  or 
combinatorial  formulas  to  generate  the  exact  null 
distribution  of  the  various  nonparametr ic  test 
statistics.  This  allows  hypothesis  testing  and 
confidence  interval  estimation  to  be  based  on  exact 
distributions  without  the  use  of  tables.  For  larger 
sample  sizes,  the  normal.  F,  and  T  distributions  are 

1APL  was  chosen  because  it  is  an  interactive 
language  that  is  especially  powerful  at  performing 
calculations  dealing  with  rank  order  statistics  and 
vector  arithmetic.  Menus  are  not  included  in  the 
workspace  designed  for  the  mainframe. 


used  to  approximate  the  distributions  of  the  test 
statistics  with  three  decimal  place  accuracy. 

Section  II  addresses  workspace  design  issues,  to 
include,  workspace  requirements  and  assumptions 
regarding  its  use.  Section  III  discusses  the  methods 
used  to  assess  the  accuracy  of  different  asymptotic 
approximations,  and  the  sample  sizes  required  for  an 
approximation  to  yield  three  decimal  place  accuracy. 
Section  IV  gives  background  information  and  discusses 
programming  methodology  for  nonpar ametr ic  tests  based 
on  single  and  paired  sample  data.  In  Section  V, 
nonparametric  tests  for  two  or  more  independent  samples 
are  considered.  Section  VI  discusses  nonparametric 
tests  for  association;  and.  Section  VII  deals  with 
nonparametric  simple  linear  regression.  Section  VIII 
recommends  other  nonparametric  tests  that  may  be  added 
to  the  workspace  and  areas  for  further  work. 

To  show  application  of  nonparametric  statistical 
methods  to  Department  of  Defense  problems,  several 
military  examples  are  worked  in  Appendix  C. 


II. WORKSPACE  DESIGN  ISSUES 

This  section  presents  a  brief  overview  of  the 
design  considerations  used  in  developing  the  APL 
workspace  for  both  the  mainframe  and  microcomputer. 

A.  EQUIPMENT  AND  SOFTWARE  REQUIREMENTS 

The  microcomputer  must  be  an  IBM  PC  or  AT 
compatible,  equipped  with  512  kilobytes  of  RAM  and  the 
APL*PLUS/PC  system  software,  release  3.25  or  later,  and 
IBM’s  DOS,  version  2.00  or  later.1  The  8087  math 
coprocessor  chip  is  not  required  to  run  this  software, 
but  will  increase  the  computational  speed. 

B.  KNOWLEDGE  LEVEL  OF  THE  USER 

The  user  is  expected  to  have  had  some  exposure  to 
APL  and  a  working  knowledge  of  nonpar ame tr ic 
statistics.  Familiarity  with  microcomputers  or  the 
Naval  Postgraduate  School  mainframe  computer  is 
assumed . 


1  The  APL  system  software  requires  144  kilobytes  of 
RAM  while  the  NONPAR  workspace  requires  an  additional 
190  kilobytes. 


11 


c. 


SELECTION  OF  TESTS 


The  nonparametric  tests  chosen  for  this  workspace 
are  some  of  the  more  widely  known,  and  are  considered 
basic  material  for  any  nonparametric  statistics  course. 
More  information  about  the  tests  can  be  found  in  any  of 
the  textbooks  that  are  referred  to  in  this  document. 

D.  MENU  DISPLAYS 

The  microcomputer’s  workspace  is  designed  around 
the  use  of  menus.  This  was  accomplished  using  the 
software  package  PC  TOOLS  from  STSC .  These  menus  are 
designed  to  guide  the  user  through  the  selection  of  the 
tests  without  an  excessive  amount  of  prompting.  The 
main  menu  displays  the  choices  available  in  the 
workspace,  while  the  test  menus  give  the  background 
information  and  options  available  for  each  test.  Help 
menus  to  provide  additional  information  about  the 
tests  are  also  available. 

E.  ORGANIZATION  OF  WORKSPACE  DOCUMENTATION 

Separate  documentation  is  included  for  the 

microcomputer’s  and  mainframe  computer’s  workspaces 
(see  Appendices  A  and  B.  respectively).  These 
appendices  explain  the  organicat con  and  operation  of 
the  workspaces.  Appendix  C,  which  provides  example 
problems  for  each  nonparametric  test,  is  applicable  to 
both  workspaces. 


r 


III. GENERAL  SAMPLE  SIZE  CONSIDERATIONS  AND  ASYMPTOTIC 

APPROXIMATIONS 

In  this  thesis,  the  term  alpha  value  is  used  in 
a  general  sense,  and  refers  to  the  probability  of 
rejecting  a  true  null  hypothesis.  The  term  P-vaiue 
refers  to  the  probability  that  a  test  statistic  will 
exceed  (or  not  exceed  in  the  lower-tailed  test)  the 
computed  value,  when  the  hypothesis  being  tested  is 
true  . 

For  selected  values,  the  exact  cumulative 
distribution  functions  (C.D.F.)  of  the  test  statistics 
are  compared  with  those  obtained  from  normal  based 
asymptotic  approximations.  The  results  of  the  compari¬ 
sons  are  used  as  a  basis  for  assessing  the  accuracy  of 
the  approximations.  In  those  cases  where  more  than  one 
asymptotic  approximation  has  been  suggested  in  the 
literature,  the  accuracy  of  each  approximation  is 
compared  over  a  range  of  desired  C.D.F.  values  and 
sample  sizes.  From  the  results,  the  most  consistently 
accurate  approximation,  and  the  sample  size  for  which 
~hac  approximation  provides  an  least  three  lecimal 
place  accuracy  is  determined. 

Once  the  accuracy  comparisons  were  completed  for  a 
specific  nonparametr ic  test,  microcomputer  capabilities 

13 


-  •  %"  »/  V  •/ 


were  considered.  In  some  cases,  generation  of  the 
exact  distribution  up  to  the  desired  sample  size  took 
too  long  or  was  not  possible  on  the  PC.  When  this 
occurred,  the  mainframe  computer  was  used  to  generate 
the  required  distributions  with  the  results  stored  in 
numerical  matrices  for  quick  recall  by  the 
nonpar amet r ic  test  programs. 


IV. TESTS  FOR  LOCATION  BASED  ON  SINGLE  AND  PAIRED- 
SAMPLE  DATA 

The  tests  assume  that  the  data  consists  of  a 
single  set  of  independent  observations  X^  or  paired 
observations  (Xi.Yj_),  i  =  l,2,...,N,  from  a  continuous 
distribution.  For  the  single  and  paired-sample  cases, 
the  null  hypotheses  are  concerned  with  the  median  of 

the  X j_  and  the  median  of  the  differences  Xj_-  , 

respectively.  The  tests  considered  are  the  Ordinary 
Sign  Test  and  the  Wilcoxon  Signed-Rank  Test. 

A.  ORDINARY  SIGN  TEST 

The  Sign  test  can  be  used  to  test  various 
hypothesis  about  the  population  median  (or  the  median 
of  the  population  of  differences).  Confidence 

intervals  for  these  parameters  can  also  be  constructed. 
As  a  final  option,  nonparametric  confidence  intervals 
for  the  quantiles  of  a  continuous  distribution  are 

offered . 

1 .  Computation  of  the  Test  Statistic 

For  single-sample  data,  the  test  statistic  X 
is  computed  as  the  number  of  observations  Xi  greater 


than  the  hypothesized  median  Mg.  For  the  paired-sample 
case,  K  is  the  number  of  differences  X^-  Y*  that  exceed 
Mg .  All  observations  Xi  (or  Xi~  Yi)  that  are  equal  to 


size 


decreased 


Mg  are  ignored  and  the  sample 
accordingly.  As  long  as  the  number  of  such  ties  is 
small  relative  to  the  size  of  the  sample,  the  test 
results  are  not  greatly  affected.  Gibbons  [Ref.  2:pp. 


1  08]  . 


2.  The  Null  and  Asymptotic  Distribution  of  K 


The  null  distribution  of  K 


binomial 


with  p  =  .5.  In  Table  1,  the  exact  values  of  the 
C.D.F  are  compared  with  the  corresponding  approximate 
values  using  a  normal  approximation  with  and  without 
continuity  correction. 


TABLE  1.  C.D.F.  COMPARISONS  FOR  THE  SIGN  TEST 

PSOBCX  i  k];  702  SAMPLE  SIZE  EQUAL  TO  24. 

TEST  STAT.  VALUE  I  51  61  71  81  9  I  10  l  11 


EXACT  C.D.F. 


I  .00331  I  .01133  I  .03196  I  .07579  I  .15373  I  .27063  I  .41941 


E2202;  N0BMAL  I  .00117  I  .00417  I  .01134  I  .02456  I  .04339  I  .06352  I  .07736 

- 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

EBBQB;  N0BM.  Id/CC  I  *.00063  1  *.00104  I  *.00114  I  *.00073  \  .00001  I  .00043  I  .00023 


PBOBCK  i  k] ;  FOB  SAMPLE  SIZE  EQUAL  TO  25. 

.  - r  ■  —  i - 1 - 1 - 1 - 1 - 1 - 

TEST  STAT.  VALUE  I  61  71  81  9  I  10  I  11  I  12 


EXACT  C.D.F. 


I  .00732  I  .02164  I  .05333  I  .11476  I  .21213  I  .34502  I  .50000 


22P.02:  MOBMAL 


. O026o  .  30774  '  .  01795  . 0o400  . .07O. 


£3202;  >10 BM.  jL  CC  i  *.00033  >*.00111  I  *.00032  l  *.00031  I  .00032  I  .00044  I  .00000 


As  can  be  seen,  for  sample  sizes  greater  than 
25,  a  normal  approximation  with  continuity  correction 


is  accurate  to  at  least  three  decimal  places. 

3 •  Hypothesis  Testing 

P-values  are  computed  for  three  basic 
hypotheses  comparing  the  median  of  the  population  or 
the  median  of  the  population  of  differences  M  with  a 
hypothesized  median  Mg.  P-values  are  taken  from  the 
cumulative  distribution  of  the  binomial  for  the 
following  tests  of  hypothesis. 

a.  One-sided  Tests 

(1)  H0 :  M  =  Mg  Versus  HI:  M  <  Mg.  The 

P-value  equals  Pr[K  ±  k] ,  where  k  is  the  computed 
value  of  the  test  statistic. 

(2)  H0:  M  ■  Mg  Versus  HI:  M  >  Mg.  The 

P-value  equals  Pr[K  j>_  k]  . 

b.  Two-sided  Test. 

(1)  H0:  M  =  Mg  Versus  HI :  M  *  Mg.  The 

P-value  equals  twice  the  smaller  value  of  a( 1 )  or 

a(2),  but  does  not  exceed  the  value  one. 

For  sample  sizes  greater  than  25,  a  normal 
approximation  with  continuity  correction  is  used. 

1 .  Confidence  Interval  Estimation 

Confidence  intervals  for  the  population 
median  are  based  on  the  ordered  observations  in  the 
sample.  For  paired-sample  data,  confidence  bounds  are 


V*. 


£■: 

>  -ji 


obtained  from  the  ordered  differences  of  the  pairs  of 
data.  A  100(1-  a  )%  confidence  interval  is  determined 
in  the  following  manner.  Let  k  be  the  number  such  that 
Pr[K  ±  k]  _<_  (  a /2) .  Then,  the  (k+1)th  and  (N-k)th 
order  statistics  constitute  the  end  points  of  the 
confidence  interval.  Gibbons  [Ref.  1:  pp.104]. 

For  computing  confidence  intervals  when 
sample  size  N  is  greater  than  25,  a  normal 
approximation  with  continuity  correction  is  used. 

Also  included  under  this  test  is  an  option  to 
generate  nonpar ametr ic  confidence  intervals  for  any 
specified  quantile  given  a  sample  size  N  from  a 
continuous  distribution.  The  end  points  of  the 
intervals  are  sample  order  statistics. 

B.  WILCOXON  SIGNED-RANK  TEST 

The  signed-rank  test  requires  the  added  assumption 
that  the  underlying  distribution  is  symmetric.  This 
test  uses  the  ranks  of  the  differences  X*-  M0  (or  Xj_- 
Yi~  Mg)  together  with  the  signs  of  these  differences  to 
determine  the  test  statistic.  Confidence  intervals  for 
the  median  can  also  be  constructed. 

1 .  Computation  of  the  Test  -Statistic 

For  single-sample  data,  the  test  statistic  W 
is  computed  as  follows. 


Let  Zi  * 


N*l 


j  1  if  Xf  M0  >  0 
|  0  if  Xj.-  M0  <_  0 


and  let  ri  ■  rank(|Xi-  M0  |  )  .  Then,  W  =  £]Zj.ri. 

1*1 

For  paired  sample  data,  W  is  calculated  in 
the  same  manner,  except  the  differences  to  be  ranked 
are  the  paired-differences  minus  the  hypothesized 
median.  Zero  differences  are  ignored  and  the  sample 
size  is  decreased  accordingly.  When  ties  occur  between 
ranks,  the  average  value  of  the  ranks  involved  are 
assigned  to  the  tied  positions.  It  has  been  shown  that 
a  moderate  number  of  ties  and  zero  differences  has 
little  effect  on  the  test  results.1 

2 .  The  Null  and  Asymptotic  Distribution  of  W 

The  exact  null  distribution  of  W  is  given  by: 

Pr [W  *  w]  -  uN(w)/2N,  w  -  0,1 ,2 . N(N+1 )/2,  where 

un(w)  is  the  number  of  ways  to  assign  plus  and  minus 
signs  to  the  first  N  integers  such  that  the  sur  e 

positive  integers  equals  w.  It  can  be  shown  'ee 
Gibbons  [Ref.  1:pp.  112])  that  uji(w),  for  successive 
values  of  N,  can  be  computed  using  the  recursive 
relationship : 

un(w)  =  ujj_  i  (  w-N  )  f  un_i(w) 


^or  more  information  on  the  effects  that  zeros  and 
tied  ranks  have  on  the  Wilcoxon  Signed-Rank  Test,  see 
Pratt  and  Gibbons  [Ref.  3]. 


- 


Exact  C.D.F.  values  were  compared  with  those 
obtained  using  the  following  asymtotic  approximations: 
student’s  T  with  (N-1 )  degrees  of  freedom  (T), 
student's  T  with  continuity  correction  (TC),  normal 
(Z),  normal  with  a  contin  lity  correction  (ZC),  the 
average  of  T  and  Z  as  suggested  by  Iman  [Ref .  4-]  ,  and 
the  average  of  TC  and  ZC . 

As  can  be  seen  in  Table  2  below,  the  average 
of  TC  and  ZC  gives  the  most  consistently  accurate 
results  with  three  decimal  place  accuracy  when  the 
sample  size  exceeds  9. 


TABLE  2.  C.D.F.  COMPARISONS  FOR  THE  WILCOXON 
SIGNED-RANK  TEST 


PROBE y  i 

»];  FOR  SAMPLE  SIZE 

EQUAL  TO 

9. 

TEST  STftT.  VALUE 

1  3 

- 

1  5 

~ 1 - 

1  6 

1 

_L- 

8  1 

9 

12 

1  14 

i. 

EXACT  C.D.F. 

1  .00977 

1  ,01953 

i 

l  .02734 

i 

1 

.04383  1 

. .  - . .  . .  -  - 1- 

.06445 

.12500 

1  .17969 

L  ... 

ESSOE; 

NORMAL 

1  ‘.00067 

1  .00046 

L 

1  .00204 

i 

1 

.00591  1 

i 

.00958 

.01324 

1  .02272 

L 

EEROS; 

NORM.  y/CC 

1  ‘.00243 

1  ".0024? 

_  L 

1  ".00167 

i  - 

"T" 

1 

.00023  1 

i 

.00269 

.00693 

1  .00806 
l 

ESSOE; 

T  DIST 

1  .00518 

1  .00608 

1  00673 

i 

1 

4- 

.00701  1 

L 

.00817 

.00326 

1  .00813 

i 

ESSOE; 

T  U/CC 

I  .00355 

1  .00276 

i 

1  .00233 

"T* 

1 

-U 

.00012  1 

‘.00010 

*.00437 

1  ".00725 

i_  .. 

ESEOS; 

ALE  T/Z 

1  .00225 

1  .00327 

1  .00433 

.00646  1 

.00883 

.01325 

1  .01545 

i  i  i 


TABLE  2.  (Continued) 
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PRCBCU  i  ud;  FOR  SAMPLE  SIZE  EQUAL  TO  10. 


TEST  STAT.  VALUE 

-r 

I 

5 

7 

1  8 

i 

10 

12 

15 

1 

17 

j-  -  -  - 

EXACT 

C.3.F. 

i 

.00977 

.01355 

1  .02441 

.04199 

.06543 

.11621 

.16113 

ERROR; 

NORMAL 

I 

1 

*.00115 

.00023 

.0009? 

i  .  .  _ 

.00476 

.00837 

.01490 

i  _ 

.01333 

ERROR; 

NORM.  U/CC 

i 

*.00270 

*.00219- 

I  *.00198 

.00043 

.00229 

.00558 

.00710 

ERROR: 

T  0I3T 

i 

.00399 

.00512 

.00525 

. 00665 

. 00661 

.00661 

.00631 

ERROR; 

T  U/CC 

1 

1 

.00249 

.00243 

.00132 

.00151 

'.00052 

“.00376 

“.00571 

L 

ERROR: 

AVE  T/2 

1 

.00142 

.00267 

.00312  1 

.00571 

.00749 

. 01075 

.01234 

ERROR; 

AVI  TC/ZC 

i 

JL 

*.00010 

_ _ _ 

.00012 

*.ooooa 

‘ 

.00097 

.00039 

. 00091 

.00070 

3 •  Hypothesis  Testing 

P-values  are  computed  for  three  basic 
hypotheses  comparing  the  median  of  the  population  or 
the  median  of  the  population  of  differences  M  with  a 
hypothesized  median  Mg  as  shown  below, 
a.  One-sided  Tests 

(1)  H0 :  M  *  Mg  Versus  HI:  M  <  Mg.  The 
P-value  equals  Pr  [W  <_  w]  ,  where  w  is  the  computed 
value  of  the  test  statistic  W. 

(2)  H0 :  M  =  Mg  Versus  Hi :  M  >  Mg.  The 


?-value  equals  PrfW  >_  w]  • 


b.  Two-sided  Test 

( 1 )  H0 :  M  =  Mg  Versus  HI :  M  t  Mg.  The 
P-value  equals  twice  the  smaller  value  of  a( 1 )  or 
a(2),  but  not  exceeding  the  value  one. 


For  sample  sizes  greater  than  9,  an  average 
of  the  normal  and  student’s  T  approximations,  each  with 
continuity  correction,  is  used.  Computations  of  the  P- 
value  for  each  of  the  alternative  hypotheses  are: 

a.  HI  :  '  <  M0 

Let  Pzc  =  Pr  C  Z  <_  (w  +.5  -  Uw)  /  and 

|w  -  uw  |  -  .5 

N<7w2  [  | w  -  fiw\  -  .5]2  *5 

N  -  1 

where  Z  is  standard  normal,  T(*j-i)  has  a  student’s  T 
distribution  with  (N-i)  degrees  of  freedom,  uw  = 
M(  N+ 1  )  /  4  and  aw2  =  (N(N+1  )( 2N+1  )/24  )  )  .  Then,  the 

P-value  for  the  test  is  (Pzc  +  (1  -  Ptc))/2  w  is 
less  than  vw  and  ( Pzc  +  pTC)/2*  otherwise.  The  above 
formulas  are  obtained  from  those  given  by  Iman  [Ref.  4] 
after  inclusion  of  a  continuity  correction. 

b.  HI :  M  >  M0 

The  P-value  equals  ((1  -  Pzc)  +  pTC)/2  if 
w  is  less  than  Vw  and  ((1  -  Pzc)  +  (1  pTc))/2» 

otherwise.  The  computation  of  Pzc  and  PTC  is  similar 

to  the  above  except  the  sign  of  the  ccntinuit'/ 
correction  is  changed. 

c .  HI  :  M  #  M0 

The  P-value  equals  twice  the  smaller 
value  of  a  or  b  above,  but  not  exceeding  the  value  one. 
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4 .  Confidence  Interval  Estimation 

For  single-sample  data,  the  confidence 
interval  for  the  population  median  is  based  on  the 
ordered  averages  of  all  pairs  of  observations  (Xs  + 
Xj  )/2  such  that  i  <_  3  .  A  1 00  (  1  -  a  )at>  confidence 
interval  is  determined  in  the  following  manner.  Let  w 
be  the  number  such  that  Pr  [w  <_  w]  <_  (  a  12).  Then,  the 
(w+1 )th  and  (m-w)th  order  statistics,  where  m  = 
N(N+1 )/2  or  the  total  number  of  paired-averages, 
constitute  the  end  points  of  the  confidence  interval. 
A  confidence  interval  for  paired-sample  data  is 
computed  in  the  same  manner,  except  the  end  points  are 
taken  from  the  paired-averages  of  the  differences  Xj_- 
Yi.  Gibbons  [Ref.  1:pp.  114-118]. 

For  computing  confidence  intervals  when 
sample  sizes  are  greater  than  9,  a  normal  approximation 
with  continuity  correction  is  used. 


V. TESTS  BASED  ON  TWO  OR  MORE  SAMPLES 


The  tests  assume  that  the  data  consists  of 
independent  random  samples  from  two  or  more  continuous 
distributions.  The  general  null  hypothesis  is  that  the 
samples  are  drawn  from  identical  populations.  The 
Mann-Whitney  and  Kruskal-Wal 1 is  tests  are  considered. 


A.  MANN-WHITNEY  TEST 

The  Mann-Whitney  test  is  based  on  the  distribution 
of  the  test  statistic  U,  which  can  be  used  to  compare 
the  equality  of  the  population  medians  or  variances 
for  two  samples.1  The  Mann-whitney  test  with  a 
modified  ranking  scheme  can  be  used  to  test  for 
equality  of  variances  if  the  population  means  or 
medians  are  assumed  to  be  equal  (Conover  [Ref.  5:pp. 
229-230]).  If  the  medians  differ  by  a  known  amount,  the 
data  can  be  adjusted  before  applying  the  test.  A 
confidence  interval  for  the  difference  in  the  medians 
of  the  two  populations  can  also  be  estimated. 

1 .  Computation  of  the  Test  Statistic 

For  “he  comparison  if  population  medians,  -.he 


test  statistic  U  is  computed  from  the  combined  ordered 


1  The  test  statistic  U  and  the  method  used  to 
compute  it  are  taken  from  Gibbons  [Ref.  1:pp.  140-141]. 


arrangement  of  observations  X^  and  Y-j  ,  i  =  1,2 . N; 

j  *  1  ,2,  .  .  .  ,M.  Let  rj_  =  rank(X^)  in  the  combined 

N 

ordered  sample  and  Rx  =  ri*  Then , 

i=  1 

U  =  Rx  -  M(M+1 ) / 2 . 

For  testing  the  equality  of  variances,  the 
computation  of  U  is  similar  except  for  the  method  of 
assigning  ranks  to  the  ordered  sample.  This  method 
ranks  the  smallest  value  1 ,  largest  value  2,  second 
largest  value  3,  second  smallest  value  4,  and  so  on  , 
by  two’s,  until  the  middle  of  the  combined  ordered 
sample  is  reached.  For  either  test,  tied  ranks  for  the 
combined  sample  are  assigned  the  average  value  of  the 
ranks  involved.  A  moderate  number  of  ties  has  little 
effect  on  the  test  results  . 

2 .  The  Null  and  Asymptotic  Distribution  of  U 

The  exact  null  distribution  of  U  is 
determined  using  a  recursion  algorithm  due  to  Harding 
[Ref.  6], 

Exact  C.D.F.  values  were  compared  with 
approximate  values  obtained  from  the  following 

asymtotic  distributions:  student’s  T  with  (n-2)  degrees 
o  f  f r eedom  wne re  n  =  N  +  M ,  the  total  number  of 


observations  in  both  samples  (T),  student’s  T  with 
continuity  correction  (TC),  normal  (Z),  normal  with 
continuity  correction  (ZC),  the  average  of  T  and  Z 
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( Iman  [Ref.  7]),  and  the  average  of  TC  and  ZC.  The 
results  for  various  sample  sizes  are  given  in  Table  3. 


TABLE  3.  C.D.F.  COMPARISONS  FOR  THE  MANN-WHITNEY 
TEST 


PROSC LI  i  u3;  FOR  SAMPLE  31 ZZ3  N  EQUAL  TO  3  AMD  rt  EQUAL  TO  3. 

i  '  ;  .  ?  ;  r  i  .  '  “I  11  —  1 — 

TEST  3TAT.  VALUE  I  14  I  17  I  13  I  21  I  23  I  27  1  23 


EXACT  C.D.f.  I  .00333  I  .01333  I  .02515  I  .04336  !  .06736-  I  .12304  I  .17005 


ER20R;  M02MAL  l '.00026  I  .00100  I  .00163  I  .00441  i  .00632  I  .01243  I  .01511 

- 1 - 1 - i - 1 - 1 - i - i - 

ER20P.;  NORM.  U/CC  I  *.00146  I  ”.00114  I '.00033  I  .00026  I  .00130  t  .00354  I  .00436 


ERROR;  T  3137  I  .00237  I  .00337  i  .00372  I  .00464  :  .00525  !  .00675  I  .00773 


ERROR;  T  U/CC  I  .00113  I  .00103  I  .00035  I  .00007  I *.00073  I *.00253  I  *.00323 

— - - : - 1 - ; _ ; _ i _ ; _ 1  ■ 

£2202;  ACE  T/Z  I  .00105  I  .00213  I  .00270  i  .00453  I  .00603  I  .00353  I  .01145 


£2202;  AVE  TC/ZC  I  *.00014  I  *.00003  I  .00004  I  .00017  !  .00026  I  .00043  I  .00054 


PR0BCU  i  u] ;  E02  SAMPLE  SIZES  N  EQUAL  TO  7  AND  M  EQUAL  TO  12. 


TEST  STAT.  VALUE 

1  14 

{- 

1  1? 

\ 

1  19 

1 

1  21 
i 

1  24 

L 

27 

t 

30 

EXACT  C.D.f. 

1  .00353 

1 

1  .01792 

i 

1  .02732 

i 

i  .04156 

1 

1  .07111 

L 

.11342 

.17012 

ERROR; 

NORMAL 

1  ".00045 

L- 

1  .00062 
f 

1  .00137 
{ 

1  .00359 

1 

1  .00701 

i 

.01097 

.01437 

ERROR; 

NORM.  U/CC 

1  ".00152 

L 

1  ".00123 

1  ".00079 

t 

i  *.00003 

1 

l  .00154 

l— 

.00322 

.00453 

ERROR; 

T  DIST 

.00199 

t. 

.00292 

. 

1  .00356 

r 

1  .00422 

i 

1  .00527 

-L  - 

.00643 

.00796 

ERROR; 

T  U/CC 

.00094 

f 

.00092 

- 

1  .00063 

1 

I  .00026 

! 

1  ".00066 

i — 

*.00173 

*.00262 

ERROR; 

A  VS  T/Z 

.00077 

.00177 

1  .00272 

1  .00331 

1  .00614 

.00370 

.01142 

ERROR: 

4vz  ::/z: 

*.  00023 

*.00013 

i  '.  00005 

'  .0001: 

1 

. 00072 

■  .. 

. JO 096 

a 
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TABLE  3.  (Continued) 


PSQ3CU  1  uj;  702  3AMPLZ  3 1 213  M  IfiUAL  TO  5  AMD  .1  2QUAL  TO  17. 


■a 

mm 

1  13 

20 

i 

i  23 

T  “  1  r 

1  27  1 

30 

Bam 

1  .00963 

1 

1  .01333 

i  .02316 

.04245 

i  .07013 

l  .  1244*0  l 

.17335 

22202;  NORMAL 

i  \  00077 

1  .00033 

L 

l  .00170 

l 

.00349 

1  .00633 

1  .01210  1 

!  1 

.01564 

22202;  NORM.  WOO 

1 *.00130 

!  -.00130 

i  ".  00037 

.00007 

:  .00133 

:  . 00444  i 

.00573 

22202;  T  OIST 

1  . 00120 

-1- 

!  .00222 

!  .00239 

.00330 

!  .00544 

:  .00774  ! 

.00355 

£2202;  T  WOO 

1  ,0001? 

1 

1  .00023 

1  .00023  1 

i  t 

.00021 

1  .00007 

! 

r.  00025  1 

‘.00051 

22202;  AVI  T /I 

I  .0002? 

!  .00120 

:  .00235  ; 

.00363 

i  .00616 

.00932  i 

.01250 

22202;  AVI  TC/IC 

i  *.  00036 

l 

i  -.00063 

;  ".  00022  . 

’ _ L 

.00014 

*  .00037 
_! _ 

.00209  i 

_ i _ 

•  00264 

P20BCLI  < 

u2 ;  FOR 

aAt-PH  31 II 3  M  iOUAL  TO  3 

AND  M  ZQ.UAL  TO  27. 

n 

1  10 

MM 

15 

1  19 

. 

23  1 

26 

IXACT  C.D.r. 

1  .00764 

1  .01630 

x 

1  .02512  1 

4  i 

.04236 

1  .07734 

•12660  1 

.17433 

£2202;  N02MAL 

1  ".00265 

1  '.00099 

L 

1  .00072  1 

i  I 

.00339 

1  .00374 

.01342  1 

» 

.01630 

£2202;  NORM.  WOO 

*.00363 

\ 

i  *.00254 

! 

C.00133  ( 

4  i 

.00089 

1  .00405 

.00665  1 
[ 

.00331 

£2202;  T  DI3T 

*.00121 

L 

I  .00031 

V 

1  .00173  1 

.00414 

1  .00741 

X  l 

.01026  1 

i 

.01253 

£2202;  T  WOO 

*.00213 

L- 

1  -.00123 
4- 

(-.00042  1 
4-  \ 

.00097 

1  .00250  1 

i  » 

.00323  1 

i 

.00390 

£2202;  AV£  T/Z 

*.00133 

1  -.00034 

t 

1  .00122  1 

X  L 

.00402 

1  .00303  1 

i  .  t 

.01134  1 

X 

.01466 

ERROR;  Ml  TC/ZC 

. 

*.00291 

i  -.00132 

l 

l ".00037  i 

X  !_ 

.00093 

i  .00327  l 

1  _  1 

.00496  l 

i 

.00610 

As 

can  be 

seen 

from  th 

e  tab 

1  e  s  .  T 

he  ave 

r  age 

of  ZC  ana  TC  gives  'he  most  consistently  accurate 


results.  For  sample  sizes  NxM  >  80,  nearly  three 


decimal  place  accuracy  is  obtained  in  all  cases. 


P-values  are  computed  for  three  basic 
hypotheses  comparing  the  medians  or  variances  of  the 
two  populations  as  shown  below. 

a.  One-sided  Tests 

( 1 )  H0 :  Mx  =  My  Versus  HI :  Mx  <  My  or 
H0  :  Vx  =  Vy  Versus  HI:  Vx  >  Vy .  The  P-value  equals 
Pr  [U  <_  u]  ,  where  u  is  the  observed  value  of  the  test 
statistic . 

(2)  H0 :  Mx  -  My  Versus  HI :  Mx  >  My  or 
H0 :  Vx  =  Vy  Versus  HI:  Vx  <  Vy.  The  P-value  equals 
Pr  [U  _>  u]  • 

b.  Two-sided  Test 

( 1 )  H0 :  Mx  =  My  Versus  HI :  Mx  #  My  or 
H0:  Vx  *  Vy  Versus  HI:  Vx  #  Vy.  The  P-value  equals 
twice  the  smaller  value  of  a( 1 )  or  a(2),  but  not 
exceeding  the  value  one. 

For  sample  sizes  NxM  greater  than  80,  the 
average  of  the  normal  and  student’s  T  approximations, 
each  with  continuity  correction,  is  used.  Computations 
of  the  P-value  for  each  alternative  hypothesis  are: 


«« 

? 

« 

> 


a.  HI  :  Mx  <  My  or  Vx  >  VY 

Let  Pzc  *  Pr[  Z  £  (u  +.5  -  ^u)/  *u3  and 


let  Ptc  *  Pr  T(n-2)  1 


|u  -  /iul  -  .5 _ 

(N+M-I)<yu2  C  |u-/iu|-.5]2 


N+M-2 


N+M-2 


where  Z  is  standard  normal,  T(n_2)  has  a  student’s  T 
distribution  with  (n-2)  degrees  of  freedom,  =  NxM/2 
and  <ru2  *  (N(M)(N+M+1  )  )/1 2.  Then  the  P-value  for  the 
test  is  (Pzc  +  O  -  PTC ) ) /2  for  u  less  than  Uu  and 
(PZC  +  pTC ) / 2 ,  otherwise.  The  above  formulas  are 
obtained  from  those  given  by  Iman  [Ref.  7]  after 
inclusion  of  the  continuity  correction. 

b.  HI :  MX  >  My  or  VX  <  Vy 

The  P-value  equals  ((1  -  Pzc)  +  ?tc)/2  if 
u  is  less  than  A*u  and  (O  “  pZC)  +  (1  -  Ptc))/2, 
otherwise.  The  computation  of  Pzc  and  P^c  is  similar 
to  the  above  except  the  sign  of  the  continuity 
correction  is  changed. 

c.  HI :  Mx  t  My  or  Vx  *  Vy 

The  P-value  equals  twice  the  smaller 
value  of  a  or  b  above ,  but  not  exceeding  the  value  one . 

4 .  Confidence  Interval  Estimation 

Confidence  intervals  for  the  difference  in 
medians,  (My  -  Mx ) ,  are  based  on  the  ordered 
arrangement  of  the  differences  (Yj-  X^),  j  ■  1,2,...,M; 


i  *  1  , 2  ,  .  .  .  ,  N  for  all  1  and  j  .  A  100(1  -  a  )% 
confidence  interval  is  determined  in  the  following 
manner.  Let  u  be  the  number  such  that  Pr  [U  <_  u]  £ 
(  a  /2).  Then,  the  (u+1 )th  and  (m-u)th  order 
statistics,  where  m  =  NxM  or  the  total  number  of 
possible  differences,  constitute  the  end  points  of  the 
confidence  interval. 

For  computing  confidence  intervals  when 


sample  sizes  NxM  are  greater  than  80,  a  normal 


approximation  with  continuity  correction  is  used. 


B.  KRUSKAL-WALLIS  TEST 


The  Kruskal -Wal 1 is  test  is  a  nonpar ametr ic  analog 


of  the  one-way  classification  analysis  of  variance  test 


for  equality  of  several  population  medians.  Gibbons 
[Ref .  1  :pp .  99]  . 


1 .  Computation  of  the  Test  Statistic 


Calculations  of  the  test  statistic  H  center 


around  the  ordered  arrangement  of  the  combined  samples 


from  which  the  sum  of  ranks  for  each  sample  is  derived. 


Let  X^-j  ,  j  =  1  , 2  ,  .  .  .  ,  ni  and  i=1,2,...,k,  be  independent 
random  samples  from  k  populations.  Let  r^-j  =  rank(Xij), 


ni  * 

Ri  =  2J  r  i  j  »  and  N  ni  •  Then , 

j  =  1  i  =  1 


H  =  (  12/(N(N+1  )  )  £  (RiS/m)  -  3(  N  +  1  ) 

(i  =  1  ) 


! 


If  ties  occur  in  the  combined  sample,  they 
are  resolved  by  assigning  the  average  value  of  the 
ranks  involved.  A  correction  based  on  the  number  of 
observations  tied  at  a  given  rank  and  the  number  of 
ranks  involved,  is  included  in  the  calculations.  A 
complete  description  of  the  correction  factor  is  given 
in  Gibbons  [Ref.  2:pp.  178-179]. 

2 .  The  Null  and  Asymptotic  distribution  of  H 

The  null  distribution  of  H  is  generated  by 
enumeration.  Each  possible  permutation  of  ranks  is 
listed  for  the  combined  sample,  and  the  corresponding  H 
value  computed.  The  frequency  distribution  of  H  is  the 
total  number  of  occurrences  of  each  distinct  H  value. 
The  H  values  are  arranged  in  increasing  order  while 
maintaining  the  frequency  pairings.  The  null 
distribution  is  obtained  by  dividing  the  cumulative 
frequencies  by  ni ! n2 ! . . • n^ ! /N ! . 

Due  to  computer  limitations,  generation  of 
the  exact  distribution  of  H  was  only  possible  for  k  *  3 
populations  with  n  *  4  observations  in  each,  and  4 
populations  with  3  observations  in  each.  Most  of  the 
distributions  were  generated  on  the  mainframe  computer 
and  saved  in  matrices  for  quick  recall  by  the  Kruskal- 
Wallis  test  program. 


Exact  C.D.F.  values  were  compared  with  the 
corresponding  approximate  values  using  the  following 


distributions:  chi-square  with  (k-1 )  degrees  of  freedom 
(C),  F  distribution  with  (k-1 )  and  (N-k)  degrees  of 
freedom  (F),  and  F  with  (k-1)  and  (N-k-1 )  degrees  of 
freedom  (FI).  The  chi-square  distribution  uses  the 
Kruskal-Wallis  H  statistic,  while  the  F  and  FI 
distributions  use  a  modified  H  statistic,  HI  = 
( (N-k)H)/(k-1  )( ( N— 1  )-H) ;  see  Iman  and  Davenport  [Ref. 
8].  As  caui  be  seen  in  Table  4,  FI  gives  the  most 
consistently  accurate  estimates. 


TABLE  4.  C.D.F.  COMPARISONS  FOR  THE  XRUSXAL- 
W ALL IS  TEST 


probzr  >  h]:  j  croup  OP  3  samples  consisting  op  *,  AND  3  DBS. 


TEST  5TAT.  VALUE 

! 

7.1439 

1 

3.7121  !  5.1318 

!  5.5985 

1 

5.0530 

I  4.2121 

j 

3.5985 

EXACT  C.2.F. 

i 

.  C0970 

1 

.01905  l  .02961 

1  .04866 

i 

.07310 

1  .12918 

i 

. 17784 

ERROR ;  CB I SQUARE 

r 

.01840 

r 

.01382  j". 01585 

i“. 01220 

1  “ 

.00184 

[  .00746 

| 

.01241 

ERROR;  P  BIST 

i 

.00304 

i 

.00736  I  .00836 

j  .01113 

i 

.01820 

1  .01696 

I 

.00990 

ERROR;  P  N/~l  DP 

i 

.00034 

i 

.00426  I  .00403 

j  .00544 

i 

.01138 

1  .00895 

i 

.00172 

PROBZB  i  hi 

FOR  A  CROUP  OP  3  SAMPLES  CONSISTING  OP  4,4 

,  AND  4  OBS. 

TEST  ST AT.  VALUE 

1 

7.6538 

i 

6.9615  |  6.5000 

1  5.6923 

i 

4.9615 

1  4.2692 

i 

3.5769 

EXACT  C.D.P. 

1 

.00762 

i 

.01939  I  .02996 

1  .04866 

i 

.08000 

j  .12190 

i 

.17299 

ERROR;  CBISQOARE 

r 

.01416 

r 

.01139  (“. 00882 

i“. 00941 

j- 

.00368 

j  .00361 

i 

.00577 

ERROR;  P  DIST 

T 

.00290 

i 

.00839  I  .01204 

(  .01100 

i 

.01272 

1  .01225 

i 

.00263 

ERROR;  P  »/~l  DP 

i 

.00149 

i 

.00600  I  .00390 

1  .00647 

i 

.00708 

1  .00592 

I- 

.00384 

PROBZB  2h]| 

FOR  A  GROUP 

OP  4  SAMPLES  CONSISTING 

OP 

3,  3, 

2,  AND  2 

OBS. 

TEST  STAT.  VALUE 

1 

7.6364 

i 

7.1818  I  7.0000 

1  6.5273 

1 

6.0182 

1  5.3818 

1 

4.8727 

EXACT  C.D.P. 

1 

.01000 

! 

.01921  1  .02921 

i  .04921 

j 

.07984 

(  .12984 

1 

.17952 

ERROR;  CRT SQUARE 

r 

.04416 

r 

.04712  (“.04269 

(“.03939 

j“ 

.03089 

[“.01604 

j" 

.00183 

ERROR :  P  DIST 

i 

. 0023U 

j 

.00250  1  .00730 

!  .00880 

I 

.01094 

1  .31123 

1 

.00902 

ERROR:  T  V/~:  OF 

— 

. 001“2 

.  ** 

. ~ ; Q i u  7 

" .  00007 

- 

J  C  2  7  5 

■  “.00570 

— 

. 00843 

PROBZB  i  hi: 

FOR  A  GROUP 

OP  4  SAMPLES  CONSISTING 

OP 

3,  3, 

3,  AND  2 

OBS. 

TEST  STAT.  VALUE 

j 

8.0152 

i 

7.6364  |  7.1315 

1  6.7273 

j 

6 . 1970 

I  5.4697 

1 

4.9697 

EXACT  C.D.P. 

j 

.00961 

i 

.01831  |  .02974 

j  .04948 

j 

.07805 

1  .12740 

1 

.17571 

ERROR;  CBI SQUARE 

j' 

.03809 

r 

.03584  (“.03748 

(“.03164 

i“ 

.02436 

(“.01306 

j 

.00168 

ERROR ;  P  DIST 

1 

.00213 

j 

.00481  (  .00441 

1  .00920 

i 

.01185 

1  .01036 

j 

.01214 

ERROR ;  P  N/~l  DP 

_T 

.00133 

r 

.00019  |“. 00269 

t  .00030 

_j_ 

.00098 

l“.  00230 

i~ 

.00091 
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A  final  accuracy  comparison  between  the  C  and  Fl 


approximations  was  conducted  by  computer  simulation  for 
5  populations  with  8  observations  each.  Initially, 
30,000  permutations  of  the  4-0  ranks  were  randomly 
generated  (no  tie  ranks  allowed),  and  the  H  statistic 


calculated  for  each  permutation.  Then  the  empirically 
determined  percentiles  Hp  for  selected  values  of  p 
between  .01  and  .18  were  compared  with  the 


approximations  given  by  the  C  and  FI  distributions. 
The  results  are  shown  in  Table  5.  It  can  be  seen  that 


the  FI  approximation  compares  well  with  the  simulated 


results,  giving  three  decimal  place  accuracy,  while 


the  C  approximation  is  less  accurate. 


TABLE  5.  C.D.F.  COMPARISONS  FOR  THE  KRUSKAL- 
W ALL IS  TEST  USING  COMPUTER  SIMULATION 


PROBZB  ah];  BASED  OB  10000  GENERATED  B'S  FOR  5  SAMPLES  OF  8  OBS.  EACB. 


TEST  ST AT.  VALUE 

1 

12.229 

1 

11.065 

1 

10.248 

1 

9.212 

j 

8.129 

1 

7.030 

i  6.232 

C.D.F.  VALUE 

1 

.01000 

1 

.02000 

1 

.03000 

1 

.05000 

i 

.08000 

1 

.13000 

j  .18000 

ERROR;  CBISQUARE 

r 

.00573 

j". 00584 

j". 00646 

r 

.00601 

r 

.00696 

i". 00432 

I". 00246 

ERROR ;  F  H/~l  DF 

i 

.00081 

u 

.00231 

1 

.00259 

t* 

.00350 

1 

.00161 

1 

.00109 

I". 00103 

PROBZB  ih];  BASED  OB 

20000  GENERATED  B'S 

FOR  5  SAMPLES  OF 

8 

OBS.  EACB. 

TEST  STAT.  VALUE 

i 

12.315 

1 

11.054 

1 

10.184 

j 

9.163 

i 

8.129 

1 

7.034 

|  6.220 

C.D.F.  VALUE 

i 

.01000 

1 

.02000 

1 

.03000 

j 

.05000 

j 

.08000 

j 

.13000 

|  .18000 

ERROR;  CBISQUARE 

r 

.00516 

r 

.00596 

j “ . 00  74  5 

r 

.00716 

r 

.00696 

r 

.00413 

(".00334 

ERROR;  E  V/"i  DF 

i 

. 30126 

1 

.  00221 

1 

.00156 

I 

. 30235 

! 

.00161 

i 

.30130 

1  ".00199 

PR03Z3  2  1 : ;  3 AS 

ED  1 N 

30000  GENERATED  B’S 

FOR  5  SAMP 

LES  IF 

i 

03S.  SACS. 

TEST  STAT.  VALUE 

i 

12.305 

1 

10.976 

1 

10.147 

1 

9.179 

] 

8.168 

i 

7.072 

I  6.265 

C.D.F.  VALUE 

1 

.01000 

1 

-4*- 

.02000 

i 

.03000 

1 

.05000 

i 

.08000 

i 

.13000 

1  .18000 
-+- - - - 

[  ERROR •,  CBISQUARE  I". 00522  I". 00684  f. 00802  00677  I". 00563  I". 00213  I". 00020 

'  ERROR t  F  u'/~\  DF  I  .00121  I*. 00143*1  . 00111*  I  *. 00273  |  .00301  |  .00344  |  .00142 
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M- 


P-values  for  the  test  H0 :  the  population 


medians  are  all  equal  versus  Hi :  at  least 

population  medians  are  not  equal,  are  computed 
Pr[H  2.  h]  *  where  h  is  the  value  of  the  observed 


two 
as : 
test 


statistic . 


The  test  statistic  is  computed  by  comparing 
each  observation  (  X j_ .  Y- )  with  all  other  observations 
( X j , Y j )  in  the  sample.  If  the  changes  in  X  and  Y  are 
of  the  same  sign,  sgn(Xj  -  Xj.)  »  sgn(Yj  ~  Y±),  the  pair 
(xi»Yi)  and  (Xj,Yj)  is  "concordant"  and  a  +1  is  scored. 
If  the  signs  are  different,  the  pair  is  "discordant" 
and  a  -1  is  scored.  Any  ties  between  either  the  X’s  or 
the  Y’s  scores  a  zero  for  that  pair.  The  sum  of  all 
scores  divided  by  the  total  number  of  distinguishable 
pairs,  (N(N-1))/2,  gives  B.  If  zeros  are  scored,  the 
denominator  is  reduced  by  a  correction  factor  which  is 
based  on  the  number  of  observations  tied  at  a  given 
rank  and  the  number  of  ranks  involved  in  each  of  the  X 
and  Y  samples.  A  complete  description  of  the 
correction  for  ties  is  given  in  Gibbons  [Ref.  2:pp. 


°  w  »-■  »-•  '' 


289] .  The  value  of  B  ranges  between  1 ,  indicating 
perfect  concordance,  and  -1,  for  perfect  discordance. 
Gibbons  [Ref.  1:pp.  209-225]. 

2 .  The  Null  and  Asymptotic  Distribution  of  B 

The  null  distribution  of  B  is  derived  from 
the  following  recursive  formula  given  in  Gibbons  [Ref. 
1  :  pp .216] . 

u( N+1 , P )  *  u(N,P)  +  u(N,P-1)  +  u( N , P-2 )  +  ...+  u(N,P-N) 


where  u(N,P)  denotes  the  number  of  P  concordant 
pairings  of  N  ranks.  This  formula  is  used  to  generate 
the  frequency  with  which  ohe  possible  values  of  P 
occur.  Division  by  N!  results  in  the  probability 
distribution  of  P.  Since,  B  =  (4P/(N(N-1 ) ) )-1 ,  the 
null  distribution  of  B  is  easily  determined. 

Exact  C.D.F.  values  were  compared  with  those 
obtained  using  a  normal  approximation,  with  and  without 
a  continuity  correction  factor  ( CC  =  6/N(N2-1),  pro¬ 
posed  by  Pittman  [Ref.  11]  for  the  Spearman’s  R  test). 
The  results  for  various  sample  sizes  are  provided  in 
Table  6.  As  can  be  seen,  for  sample  sizes  greater  than 
13,  a  normal  approximation  with  continuity  correction 
provides  three  decimal  place  accuracy. 


K5 


V.vV 


I 


TABLE  6.  C.D.F.  COMPARISONS  FOR  KENDALL'S  B 


PR03C3  >  W;  ?02  SAMPLE  3122  EQUAL  TO  13. 


7IST  3TAT.  VALUE  I  0.5123  !  0.4615  I  0.4103  i  0.3590  I  0.3333  !  0.2564  I  0.2303 


XACT  C.D.F. 


22208;  NORMAL 


.00743  I  .01524 


.02363  i  .04999  I  .06443  I  . 12593  -i  .15309 


.00014  I  .00121  I  .00313  I  .00620  I  .00303  I  .01473  I  .01703 


■.00013  I  .00073  I  .00223  !  .00497  I  .00653  !  .01223  1  .014’ 


PR0BC8  2  b] ;  FOR  SAMPLE  3IZ2  EQUAL  TO  14. 


ACT  C.9.F. 


TEST  3TAT.  VALUE  I  0.4725  !  0.4236  !  0.4066  I  0.2626  I  0.2367  !  0.2527  I  0.2033 


.00964  I  .01773  I  .02353  !  .03373  I  .07353  i  .11656  I  .16541 


E220E;  NORMAL 


.00035  I  .00140  I  .00213  I  .00431  I  .00339 


ERROR;  NORM.  W/CC  I  .00003  I  .0003 


.00742  !  .01057  !  .01371 


othesis  Testin 


P-values  for  tests  of  no  association  between 
X  and  Y  are  computed  for  three  types  of  alternative 
hypotheses.  Because  the  distribution  of  B  is  symmetric, 
all  probabilities  can  be  taken  from  the  upper  tail 
using  the  absolute  value  of  b,  the  observed  value  of 
the  test  statistic.  Linear  interpolation  is  used  when 
b  lies  between  tabulated  values.  The  P-values  are 
ocmcuted  is  follows. 


a.  One-Sided  Alternatives 


The  one-sided  alternative  tested  depends 
on  the  sign  of  b.  A  positive  b  will  automatically  test 


for  direct  association  or  concordance,  while  a  negative 
b  will  test  for  indirect  association  or  discordance. 
The  P-value  equals  Pr[B  _>  lb|]. 

b.  Two-Sided  Alternative 

The  P-value  equals  twice  the  probability 
computed  for  the  one-sided  hypothesis. 

For  sample  sizes  greater  than  12,  a  normal 
approximation  with  continuity  correction  is  used.  The 
approximate  P-values  are  then: 

1  -  ?r[  Z  <_  (  (  |  b  j  -  CC)  -  ^b)/  <rb]  ,  where  Z 
is  standard  normal,  CC  is  the  continuity  correction, 
=  0,  and  =  ( 4N  +  10)/9N(N-l),  for  the  one¬ 

sided  test  and  twice  this  P-value  for  the  two-sided 
test . 

B.  SPEARMAN’S  R 

The  Spearman’s  R  Test  requires  the  added 
assumption  that  the  underlying  bivariate  distribution 
is  continuous.  The  test  measures  the  degree  of 
correspondence  between  rankings,  instead  of  the  actual 
variate  values,  and  can  be  used  as  a  measure  of 
association  between  X  and  Y.  Gibbons  [Ref.  1:pp.  226]. 

'  .  Computation  of  the  Test  Statistic 

The  test  statistic  R  is  computed  in  the 
following  manner.  Let  r^  =  rank(Xi)  and  si  =  rank(Yi) 
and  Di  =  r^  -si-  Then, 
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N 

6  Y  Di2 

R  =  1  -  _ i  =J _ 

"n(n2  -  1 ) 

where  M  is  the  size  of  the  sample.  If  ties  occur  in  X 
or  Y,  they  are  resolved  by  assigning  the  average  value 
of  the  ranks  involved.  A  correction  factor,  based  on 
the  number  of  observations  tied  at  a  given  rank  and  the 
number  of  ranks  involved,  is  included  in  the 
calculations.  A  complete  description  of  the 

correction  factor  is  given  in  Gibbons  [Ref.  2:pp.  279]  . 
The  value  of  R  ranges  between  1 ,  indicating  perfect 
direct  association,  and  ,  for  perfect  indirect 
association.  Gibbons  [Ref.  l:pp.  226-235]. 

2 .  The  Null  and  Asymptotic  Distribution  of  R 

The  null  distribution  of  R  for  a  given  sample 
size  N  is  generated  by  enumeration.  The  method,  as 
presented  in  Kendall  [Ref.  9],  involves  generation  of 
an  N  by  N  array  of  all  possible  squared  differences 
between  any  two  paired  ranks  of  X  and  Y.  All  N! 
permutations  of  N  ranks  are  used  to  index  values  from 
the  array.  The  sum  of  these  indexed  values  for  each 
permutation  gives  rise  to  M'  sum  di*  squared  iifferencss 
which  are  then  converted  to  the  R  statistic.  The 


frequency  distribution  of  R  is  the  total  number  of 
occurrences  of  each  distinct  value  of  R  divided  by  N! . 


Due  to  mainframe  computer  memory  limitations 
in  the  APL  environment,  generation  of  the  distribution 
of  R  was  limited  to  sample  sizes  of  7  or  less.1  Using 
tables,  provided  by  Gibbons  [Ref.  2:pp.  417-418]  to 
supplement  computer  computations,  a  numerical  matrix, 
called  PMATSP,  was  created  to  store  the  cumulative 
distributions  of  R  for  sample  sizes  less  than  1 1 .  This 
matrix  allows  for  quick  recall  of  cumulative 
probabilities  by  the  Spearman’s  R  Test  program. 

Exact  C.D.F.  values  were  compared  with  those 
obtained  using  a  student’s  T  approximation  with  (N-2) 
degrees  of  freedom  (see  Glasser  and  Winter  [Ref.  13]), 
and  a  normal  approximation.  Both  normal  and  T 
approximations  were  computed  with  and  without  a 
continutity  correction  factor,  CC  =  6/N(N2-1 )  (Pittman 
[Ref.  11]).  From  the  results  presented  in  Table  7,  the 
most  consistently  accurate  approximation  is  given  by 
the  T  distribution  with  a  correction. 

3 .  Hypothesis  Testing 

P-values  for  tests  of  no  association  between 
X  and  Y  can  be  computed  for  three  types  of  alternative 
hypotheses.  Because  the  distribution  of  R  is  symmetric, 
ail  probabilities  are  taken  from  the  upper  tail  using 

1  The  memory  capacity  of  the  mainframe  computer  in 
the  APL  environment  is  limited  to  2.5  megabytes. 
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TABLE  7 


C.D.F.  COMPARISONS  FOR  SPEARMAN’S  R 


PS0BC2  2 

r 3;  £02 

SAMPLE  SITE 

EQUAL  TO 

3. 

TEST  3TAT.  VALUE 

1  0.7333 

1  0.7167 

1  0.6667  1 

0.6000  1 

0.5333 

1  0.4333  1 

1  » 

0.3500 

EXACT  C.O.E. 

!  .  00361 

1  .01343 

1  .02344  1 

.04340  ! 

.  0737  6 

1  .12436  ! 

.17323 

£2202;  NORMAL 

1  *.00475 

l 

!  '.00230 

1 

1  *.  00022  i 

.00356  ! 

.00305 

l  .01479  1 

;  1 

.01313 

£2202:  N02M.  U/CC 

1  *.00553 

1  *.00413 

T 

1  *.00135  1 

.00123  1 

.00433 

1  .01023  1 

.01236 

1  .00235 

1  .00352 

1  .00451  1 

.00453  I 

.00415 

1  .00233  1 

.00138 

£2202;  T  W CC 

1  .00153 

1  .00203 

1 

1  .00251  1 

1  » 

.00175  1 

1 

.00042 

1  *.00212  1 

L 

*.00473 

PP.03C2  i 

p2:  j 02 

:AMPL£  SIZE 

EQUAL  70 

10. 

TEST  3TAT.  VALUE 

1  0.7455 

1  0.6364  1 

0. 5636  1 

0.4303 

0.4061  ! 

0.3333 

aXAc i  c • 0 • i ■ 

!  .00370 

1  .01343 

1  • 

I  .02722  ! 

.04314  ! 

.07741 

■Raw 

.17437 

£220?.;  NORMAL 

1  *. 00336 

i  *.  00220 

! '.  00031  ! 

mm 

.00700 

.01216  ! 

.01572 

12202;  NO EM.  U/CC 

1  *.00457 

1 

1  ".00327 
| 

i  ".00210  1 

l  1 

.00095  l 

1 

.00451 

.00367  i 
|- 

.01123 

12202;  T  DIST 

1  .00204 

1 

1  .00236 

I 

1  .00326  1 

f  1 — 

.00323  1 

L 

.00253 

.00153  1 

L 

.00107 

12202;  T  U/CC 

1  .00144 

1 

1  .001*35 

1 

1  .00134  1 

j _ [ _ 

.00114  r 

_ 1 _ 

'.00035 

! 

*.00230  1 

*.00361 

the  absolute 

value 

of  r  , 

the  observed 

value 

of  the 

test 

statistic.  The  P-values  are  computed  as  follows. 


a.  One-Sided  Alternatives 

The  one-sided  alternative  tested  depends 

on  the  si^n  of  r.  A  positive  r  will  test  for  direct 
3.S  c>oc  ..on,  .1.  ^  ne  ^  el  t  i  vs  .p  o  0  s  u  5  ^  ^  ^  i  r  0  c  * 


association . 


The  P-value  equals  Pr  [R  _>  |r|]. 


b.  Two-Sided  Alternative 


The  P-value  equals  twice  the  probability 
computed  for  the  one-sided  hypothesis. 

For  sample  sizes  greater  than  13,  an 
approximation  based  on  the  student’s  T  distribution 
with  (N-2)  degrees  of  freedom  and  continuity 
correction,  is  used.  The  P-values  are: 

1  -  Pr  [  T(N_2)  1  (  (  |  r  |  -  CC)  -  ur )/  <7r], 

where  T(n-2)  denotes  the  T  distribution  with  (N-2) 
degrees  of  freedom,  CC  is  the  continuity  correction. 

=  0.  and  <rr2  =  (1  -(  |  r  | -CC  )2  )  /  (  n-2  )  ,  for  the  one¬ 
sided  test,,  and  twice  this  P-value  for  the  two-sided 
test.  Gibbons  [Ref.  l:pp.  218]. 


VII .NONPARAMETR I C  SIMPLE  LINEAR  REGRESSION1 
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Nonparametr ic  Linear  Regression  assumes  that  the 
data  consists  of  independent  pairs  of  observations  from 
a  bivariate  distribution  and  that  the  regression  of  Y 
on  X  is  linear.  The  program  estimates  linear 
regression  parameters  based  on  the  data  samples.  It 
then  allows  the  user  to  input  X  values  to  predict  the  Y 
values.  Hypothesis  testing  and  confidence  interval 
estimation  for  the  slope  of  the  regression  equation  is 
offered.  If  the  estimated  slope  lies  outside  the 
confidence  interval,  an  alternate  regression  equation 
is  offered  with  an  opportunity  to  input  X  values  to 
predict  the  corresponding  Y  values. 


A.  COMPUTATION  OF  THE  ESTIMATED  REGRESSION  EQUATION 

The  least  squares  method  is  used  to  estimate  A  and 
B  in  the  regression  equation  Y*  =  A  +  BX*  +  ei 
( i= 1 , 2 , . . .N ) ,  where  ei  (unobservable  errors)  are 
assumed  to  be  independent  and  identically  distributed. 
A  and  B  are  computed  from  the  following  equations: 


Except  for  program  design  considerations,  the 
information  and  concepts  provided  in  the  section  are 
paraphrased  from  Conover  [Ref.  12:pp.  263-271]. 


N  N  N 

N  £  xiYi  -  Exi£Yi 

i  =  1 _ i  =  1  i  =  1 

N  /  N  \2 

»S"i2  -  ( E  xi) 

i  =  1  'i  =  1  f 

N  N 

£  Yi  -  B  £  Xi 

i=J _ i=J _ 

N 

B.  HYPOTHESIS  TESTING 

P-values  for  testing  hypotheses  about  the  slope  of 
the  regression  equation  are  based  on  the  Spearman’s 
rank  correlation  coefficient  R  between  the  Xi  and  U-j_  = 
Yi  -  3gXj_,  where  Bg  is  the  hypothesized  slope.  The 

appropriate  one-sided  test  of  hypothesis,  H0 :  B  =  Bg 

versus  HI:  B  <  B@  or  HI:  B  >  Bg ,  is  automatically 
chosen  based  on  the  sign  of  the  computed  test  statistic 
r  (positive  r  tests,  HI :  B  >  Bg;  negative  r  tests,  HI: 

B  <  Bg ) .  The  P-value  is  computed  as:  Pr  [R  j>  |r|]. 

P-values  for  two-sided  tests,  H0:  B  =  Bg  versus  HI :  B  # 
Bg ,  are  also  presented. 

For  sample  sizes  N  greater  than  10,  P-values  are 
approximated  using  a  T  distribution  with  (N-2)  degrees 
of  freedom  and  continuity  correction. 

C.  CONFIDENCE  INTERVAL  ESTIMATION 

100(1-  q  )%  confidence  bounds  for  the  slope 
parameter  B  are  determined  as  follows.  The  n  possible 


slopes,  Sj_j  =  (Yi-Yj  )/( Xj_-Xj  ) ,  are  computed  for  all 
pairs  of  data  (Xi,Yj_)  and  (X-j,Yj)  such  that  i  <  J  and 
Xi  #  Xj  and  rearranged  in  increasing  order  to  give  s( 1 ) 
<_  s(2)  <_  ...  <_  s(n).  Let  w  be  the  (1-  a/2)  percentile 
of  the  distribution  of  Kendall’s  statistic  with  sample 
size  n.1  Let  d  be  the  largest  integer  less  than  or 
equal  to  (n-w)/2  and  u  the  smallest  integer  greater 
than  or  equal  to  (n+w)/2  +  1.  Then  and  S^u)  are 
the  desired  lower  and  upper  confidence  bounds, 
respectively . 

For  sample  sizes  larger  than  13,  a  normal 
approximation  with  continuity  correction  is  used  to 
estimate  the  confidence  intervals. 

If  the  slope  of  the  estimated  regression  equation 
does  not  lie  within  the  computed  confidence  interval, 
the  program  automatically  calculates  a  new  regression 
equation  where  the  slope  is  the  median  of  the  two-point 
slopes  Sij  and  the  intercept  is  the  difference  of  the 
medians  of  the  X  and  Y  samples.  My-  Mx«2 


"'Kendall’s  statistic  is  defined  here  as  Nc  -  N^, 
where  Nc  is  the  number  of  concordant  pairs  of 
observations  and  is  the  number  of  discordant  pairs. 
Conover  [REF.  1 2  :  pp  .  256]. 

^This  procedure  is  recommended  by  Conover  [REF. 
1 2 :pp .  256] . 
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VIII .AREAS  FOR  FURTHER  WORK 


To  create  a  more  versatile  and  powerful  software 
package,  the  NONPAR  workspace  could  be  expanded  to 
include  some  or  all  of  the  following  nonparametr ic 
tests:  tests  for  randomness  based  on  runs,  Chisquare 
and  Kolmogor ov-Smi rnov( K-S )  Goodness-of-f i t  tests, 
Chisquare  and  K-S  general  two  sample  distribution 
tests,  Chisquare  test  for  independence.  and  the 
Friedman  test  for  association. 
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APPENDIX  A 


DOCUMENTATION  FOR  THE  MICROCOMPUTER  WORKSPACE 

1  .  *  General  Information 

This  appendix  describes  the  organization  and 
operation  of  the  IBM-PC  (or  compatible)  version  of  the 
workspace.  Appendix  C  continues  from  where  this 
appendix  leaves  off,  to  walk  the  user  through  each  test 
by  working  practical  examples. 

Before  proceeding  any  further,  the  user  should 
refer  to  section  II  (Workspace  Design  Issues)  for 
general  information  about  workspace  requirements  and 
assumptions  regarding  its  use. 

To  get  started,  enter  the  APL  environment  in  the 
usual  manner  and  load  the  NONPAR  workspace. 

2 .  Workspace  Menus 

This  workspace  is  designed  around  the  use  of 
menus.  They  guide  the  user  through  the  selection 
process  of  choosing  a  nonparametric  test  and  a  test 
option.  Three  types  of  menus  are  used;  the  main  menu, 
test  menus,  and  help  menus. 

a.  The  Main  Menu 

Within  moments  of  loading  the  NONPAR 
workspace,  the  main  menu  will  appear.  It  is  titled 
Nonparametric  Statistical  Tests.  This  menu  presents 
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general  information  about  the  workspace.  Its  primary 
purpose  is  to  list  the  choices  of  nonpar ametric  tests 
available  and  provide  an  option  which  allows  the  user 
to  exit  the  main  menu  into  APL  to  copy  data  into  the 
workspace  or  return  to  DOS.  Each  test  choice  is  listed 


with  some  information  about  the  test's  area  of 
application.  To  make  a  selection  from  the  menu,  move 
the  cursor  (using  the  cursor  keys)  to  highlight  the 
desired  choice,  and  press  enter.  As  a  reminder  to  the 
user.  a  footnote  at  the  bottom  of  the  screen  describes 
the  procedure  for  entering  a  choice.  Once  a  test  has 
been  selected  from  the  main  menu,  a  sub-menu 
appropriate  to  the  test  appears.  To  exit  from  any  menu 
back  to  the  main  menu,  press  the  Escape  key. 
b.  Test  Menus 

The  title  of  the  test  menu  is  the  name  of  the 
nonparametric  test  chosen.  The  text  portion  of  the 
menu  gives  a  general  overview  of  the  test,  to  include, 
the  method  used  to  compute  the  test  statistic,  and  a 
description  of  the  various  options  that  may  be 

exercised.  The  third  section  consists  of  the  list  of 

test  options  available.  These  options  include 
returning  to  the  main  menu  or  choosing  the  help  menu. 
Test  menus  may  have  options  listed  in  single  or 

multiple-paged  formats.  The  comment  in  the  final  block 
of  the  menu  lets  the  user  know  if  a  certain  menu  is 
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multiple-paged  or  not.  To  make  a  selection  from  a 
multiple-paged  menu,  use  the  page-up  or  page-down  key 
to  locate  the  desired  option.  Proceed  with  the  scroll 
keys  to  highlight  the  choice,  and  press  enter.  Once  a 
test  option  is  entered,  the  user  is  prompted  to  input 
the  data  required  to  run  the  test.  When  the  option  for 
more  information  is  selected,  the  help  menu  is 
displayed. 

c.  Help  Menus 

The  title  of  the  help  menu  usually  begins 
with  the  words  "More  Information  About..."  followed  by 
the  title  of  the  nonpar ametr ic  test.  The  text  portion 
of  the  menu  explains  the  test  and  its  options  in 
greater  detail.  No  choices  are  offered  in  the  menu. 
To  return  to  the  test  menu,  press  any  key. 
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APPENDIX  B 


DOCUMENTATION  FOR  THE  MAINFRAME  COMPUTER  WORKSPACE 

1 .  General  Information 

This  appendix  describes  the  organization  and 
operation  of  the  mainframe  computer  workspace.  To  load 
a  copy  of  the  NONPAR  workspace  from  the  APL  library, 
enter  the  APL  environment  and  type:  )LOAD  9  NONPAR. 
Within  a  few  moments  the  variables  LIST  and  DESCRIBE 
are  displayed  on  the  screen.  These  variables  provide  a 
description  of  the  workspace. 

2 .  The  NONPAR  Workspace 

The  NONPAR  workspace  consists  of  seven 
programs  which  call  several  subprograms  during  their 
execution.  The  exact  syntax  for  each  test  and  its 
corresponding  nonparametric  test  name  is  given  in  the 
following  format: 

SYNTAX:  Nonparametric  Test  and  Application. 

a.  SIGN:  Ordinary  Sign  Test  for  Location  in  Single 
and  Paired-sample  Data. 

b.  WILCOX:  Wilcoxon  Signed-rank  Test  for  Location 
in  Single  and  Pair ed-samo le  Data. 

c.  MANNWHIT:  Mann-Whitney  Test  for  Equal  Medians  or 
Variances  in  Two  Independent  Samples. 

d.  KRUSKAL:  Kruskal -Wal 1 i s  Test  for  Equal  Medians  in 
K  Independent  Samples. 
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e.  KENDALL:  Kendall’s  B;  Measure  of  Association  for 

Paired-sample  Data. 

f.  SPEARMAN:  Spearman’s  R;  Measure  of  Association 

Between  Rankings  of  Paired  Data. 

g.  NPSLR:  Nonparametr ic  Simple  Linear  Regression; 

Least  Squares. 

The  list  presented  above  can  be  displayed  at  any 
time  by  typing:  LIST. 

For  each  test  program,  there  exists  a  ROW  variable 
that  gives  a  full  description  of  the  test  and  the 
various  options  that  may  be  exercised.  To  display  any 
of  the  HOW  variables,  just  enter  the  test  program’s 
name  with  the  suffix  HOW  appended  (i.e.  SIGNHOW). 

A  test  is  run  by  entering  the  program's  name.  The 
user  is  immediately  prompted  to  input  data.  Enter 
numerical  data  separated  by  spaces  or  as  a  variable  to 
which  the  numbers  have  been  previously  assigned. 
Several  of  the  tests  require  a  considerable  amount  of 
prompting  before  all  the  necessary  data  has  been 


APPENDIX  C 


WORKSPACE  FAMILIARIZATION  THROUGH  PRACTICAL  EXAMPLES 

1 .  General  Information 

This  appendix  applies  to  both  the  mainframe  and 
microcomputer  workspaces.  Its  purpose  is  to  acquaint 
the  user  with  the  organization  of  the  programs  and  the 
type  of  prompts  to  be  expected. 

Extensive  error  checking  has  been  included  in  the 
programs  to  ensure  that  the  data  is  of  the  proper  form. 
Should  a  program  become  suspended,  clear  the  state 
indicator  by  entering:  )RESET,  check  over  the  data  for 
errors,  and  restart  the  program.  330  kilobytes  of 
computer  memory  are  needed  to  load  APL  and  the  NONPAR 
workspace;  to  avoid  filling  up  the  remaining  workspace 
area,  the  user  should  minimize  data  storage  in  the 
NONPAR  workspace.  To  exit  a  program  at  any  time,  press 
the  Control  and  Escape  keys,  simul taneously . 

2 .  Practical  Examples 
a.  Sign  Test 

( 1 )  Description  of  Problem  1 .  A  Sinclair  mine 
is  manufactured  to  have  a  median  explosive  weight  of 
not  less  than  16  ounces.  The  explosive  weights  of  15 
mines,  randomly  selected  from  the  production  line,  were 


recorded  as  follows: 


16.2  15.7  15.9  15.8  15.9  16  16.1 


15.8  15.9  16  16.1  15.7  15.8  15.9  15.8. 

(a)  Is  the  manufacturing  process  packing 
enough  explosives  in  the  mines? 

(b)  What  range  of  values  can  be  expected 
for  the  median  of  the  explosive  weights  9Q%  of  the 
time . 

(2)  Solution .  To  see  if  the  manufacturing 
process  is  meeting  the  specifications,  we  test  the 
hypothesis  Ho:  M  =  16  versus  HI :  M  <  16. 

( 3 )  Workspace  Decision  Process. 

(a)  Microcomputer:  Choose  the  Sign  Test 
from  the  main  menu,  and  the  option.  Single  Sample; 
Test  H0 :  M  =  Mo  versus  HI:  M  <  Mo,  from  the  test  menu. 
Skip  to  the  Program  Interaction  section  below. 

(b)  Mainframe:  Enter  SIGN  at  the  keyboard 
and  receive  the  prompt: 

DID  YOU  ENTER  THIS  PROGRAM  FOR  THE  SOLE 
PURPOSE  OF  GENERATING  CONFIDENCE  INTERVALS  FOR  A 
SPECIFIED  SAMPLE  SIZE  AND  QUANTILE?  (Y/N). 

Enter  N  ( If  Y  is  entered,  the  user  will 
go  directly  to  this  last  option  of  the  test).  The  next 
prompt  is: 

THE  NULL  HYPOTHESIS  STATES  -  THE 
POPULATION  MEDIAN  (M)  IS  EQUAL  TO  THE  HYPOTHESIZED 
MEDIAN  (Mo);  H0 :  M  =  Mo.  WHICH  ALTERNATIVE  DO  YOU  WISH 


TO  TEST?  ENTER:  1  FOR  HI:  M  <  Mo:  2  FOR  HI:  M  >  Mo: 
3  FOR  HI :  M  *  Mo. 

Enter  1.  The  next  prompt  is: 

ENTER:  1  FOR  SINGLE-SAMPLE  PROBLEM:  2 

FOR  PAIRED-SAMPLE  PROBLEM. 

Enter  1 . 

(4)  Program  Interaction.  The  prompt  is: 

ENTER  THE  DATA  (MORE  THAN  TWO  OBSERVATIONS 
ARE  REQUIRED). 

Enter  the  data  separated  by  spaces  or  as  a 
variable  to  which  the  data  has  been  previously 
assigned.  The  next  prompt  is: 

ENTER  THE  HYPOTHESIZED  MEDIAN. 

Enter  16.  The  following  is  dislayed. 
COMPUTATIONS  ARE  BASED  ON  A  SAMPLE  SIZE  OF: 

13. 

THE  TOTAL  NUMBER  OF  POSTIVE  SIGNS  IS:  3. 

THE  P-VALUE  FOR  H0 :  M  ■  16  Versus  HI: 
M  <  1 6  IS:  .0461 . 

Consider  a  significance  level  of  .05. 
Since  the  P-value  of  .0461  is  less  than  .05,  we  reject 
H0:  M  =  16  in  favor  of  Hi :  M  <  16  and  conclude  that  the 

manufacturing  process  is  not  packing  enough  explosives 
in  the  Sinclair  mine.  The  next  prompt  is: 

WOULD  YOU  LIKE  A  CONFIDENCE  INTERVAL  FOR 


THE  MEDIAN?  (Y/N). 


Enter  Y  ( If  N  is  entered,  the  progam  asks 
if  confidence  intervals  for  a  quantile  are  desired). 
The  next  prompt  is: 

ENTER  THE  DESIRED  CONFIDENCE  COEFFICIENT; 
FOR  EXAMPLE:  ENTER  95,  FOR  A  95%  CONFIDENCE  INTERVAL. 

Enter  90.  The  following  is  displayed. 

A  90%  CONFIDENCE  INTERVAL  FOR  THE  MEDIAN  OF 
THE  POPULATION  IS:  (  1  5 . 8  MEDIAN  1  16  ). 

The  next  prompt  is: 

WOULD  YOU  LIKE  CONFIDENCE  INTERVALS  FOR  A 
SPECIFIED  QUANTILE?  (Y/N). 

To  see  the  form  of  the  results,  we  generate 

confidence  intervals  for  the  30th  quantile.  Sample 

size  is  automaticly  set  at  the  number  of  data  points 
entered  earlier.  Enter  Y  (If  N  is  entered,  the 
mainframe  program  ends;  or,  the  Sign  test  menu 
reappears).  The  next  prompt  is: 

ENTER  DESIRED  QUANTILE;  FOR  EXAMPLE:  ENTER 
20,  FOR  THE  20TH  QUANTILE. 

Enter  30.  The  following  is  displayed. 

ORDER  STATISTICS  j_  COEFFICIENTS 
3  8  !  .823160 

2  9  !  .949190 

1  0  .991600 

*****  THIS  TABLE  GIVES  CONFIDENCE  COEF¬ 
FICIENTS  FOR  VARIOUS  INTERVALS  WITH  ORDER  STATISTICS  AS 
END  POINTS  FOR  THE  30TH  QUANTILE. 


The  mainframe  program  ends.  The  menu- 
driven  microcomputer  program  pauses  for  input  from  the 
keyboard  by  prompting: 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  the  Sign  test  menu 

reappears . 

b.  Wilcoxon  Signed-rank  Test 

( 1  )  Description  of  Problem  2.  A  special  train¬ 
ing  program  is  being  considered  to  replace  the  regular 
training  that  Radio  Telephone  Operators  receive.  In 
order  to  evaluate  the  effectiveness  of  the  new  training 
program,  proficiency  tests  were  given  during  the  third 
week  of  regular  training.  Twenty-four  trainees  were 
chosen  at  random  and  grouped  into  twelve  pairs  based  on 
proficiency  test  scores.  One  member  of  each  pair 
received  specialized  training  while  the  other  member 
received  regular  training.  Upon  graduation,  the 
proficiency  tests  were  given  again  with  the  following 
results . 

Specially  Trained  Group  (X):  60  50  55  71 

45  59  64  49  61  54  47  70 

Regularly  Trained  Group  ( Y ) :  40  46  60  55 

49  57  5!  55  45  59  40  35 

(a)  Does  the  special  training  program 
ensure  higher  scores? 


(b)  By  what  range  of  values  can  the 
scores  of  the  two  groups  be  expected  to  differ  95^  of 
the  time? 

(2)  Solution .  To  test  the  hypothesis  that  the 
special  training  program  raises  profficiency  scores,  we 
test  H0:  M(X-Y)  »  0  versus  His  M(X-Y)  >  0. 

( 3 )  Workspace  Decision  Process 

(a)  Microcomputer:  Choose  the  Wilcoxon 

Signed-rank  Test  from  the  main  menu,  and  the  option. 
Paired-sample;  Test  H0 :  M  =  Mo  versus  HI :  M  >  Mo,  from 
the  test  menu.  Skip  to  the  Program  Interaction  section 
below. 

(b)  Mainframe:  Enter  WILCOX  at  the  key¬ 
board  and  receive  the  prompts: 

THE  NULL  HYPOTHESIS  STATES  -  THE 
POPULATION  MEDIAN  (M)  IS  EQUAL  TO  THE  HYPOTHESIZED 
MEDIAN  (Mo);  H0 :  M  *  Mo.  WHICH  ALTERNATIVE  DO  YOU  WISH 
TO  TEST?  ENTER:  1  FOR  HI :  M  <  Mo;  2  FOR  HI :  M  >  Mo; 
3  FOR  HI :  M  ^  Mo. 

Enter  2.  The  next  prompt  is: 

ENTER:  1  FOR  SINGLE-SAMPLE  PROBLEM;  2 
FOR  PAIRED-SAMPLE  PROBLEM. 

Enter  2. 

(4)  Program  Interaction.  The  prompt  is: 

ENTER  X  DATA  (MORE  THAN  TWO  OBSERVATIONS 


ARE  REQUIRED). 


Enter  the  X  data  separated  by  spaces.  The 


next  prompt  is: 


ENTER  Y  DATA  (NUMBER  OF  Y  ENTRIES  MUST 


EQUAL  NUMBER  OF  X  ENTRIES). 


Enter  the  Y  data.  The  next  prompt  is: 


ENTER  THE  HYPOTHESIZED  MEDIAN  FOR  THE 


DIFFERENCES  OF  THE  PAIRED  DATA. 


Enter  0.  The  following  is  dislayed. 


COMPUTATIONS  ARE  BASED  ON  A  SAMPLE  SIZE  OF: 


THE  TOTAL  SUM  OF  POSITIVE  RANKS  IS:  60.5. 


THE  P-VALUE  FOR  COMPARING  THE  MEDIAN  OF  THE 


POPULATION  OF  DIFFERENCES  TO  THE  HYPOTHESIZED  MEDIAN, 


H0 :  M(X-Y)  =  0  Versus  HI:  M(X-Y)  >  0,  IS:  .0505. 


Consider  a  significance  level  of  .05 


Since  the  P-value  of  .0505  is  greater  than  .05,  we  do 


not  reject  the  null  hypothesis  that  the  two  training 


cources  are  equally  effective.  However,  due  to  the 


closeness  in  values,  the  choice  of  rejecting  or  not 


rejecting  the  null  hypothesis  is  strictly  a  Judgement 


call.  The  next  prompt  is: 


WOULD  YOU  LIKE  A  CONFIDENCE  INTERVAL  FOR 


THE  MEDIAN?  (Y/N). 


Enter  Y  (If  N  is  entered,  the  mainfr; 


progam  ends;  or,  the  Wilcoxon  test  menu  reappears). 


The  next  prompt  is 


J  Uii uv wm  A Va.  A 


yyiL 


ENTER  THE  DESIRED  CONFIDENCE  COEFFICIENT; 
FOR  EXAMPLE;  ENTER  95,  FOR  A  9556  CONFIDENCE  INTERVAL. 

Enter  95.  The  following  is  displayed. 

A  9556  CONFIDENCE  INTERVAL  FOR  THE  MEDIAN  OF 
THE  POPULATION  OF  DIFFERENCES  IS; 

(  -1  <_  MED  I  AN  ( X-Y )  <_  16.5  ). 

The  mainframe  program  ends.  The  menu- 
driven  microcomputer  program  pauses  for  input  from  the 
keyboard  by  prompting; 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  the  Wilcoxon  test  menu 

reappears . 

c.  Mann-Whitney  Test  for  Equality  of  Medians 

(1 )  Description  of  Problem  3.  A  group  of  Army 
and  Navy  officers  were  given  the  Defense  Language 
Aptitude  test.  From  the  results,  14  Army  and  17  Navy 
officers’  scores  were  randomly  selected.  These  scores 
are  listed  below. 

Army  (X);  35  30  55  51  28  25  16  63  60  44  20 

42  47  38. 

Navy  (Y):  54  26  41  43  37  34  39  50  46  49  45 

33  29  36  38  42  34. 

(a)  Is  there  sufficient  evidence  to  claim 
that  Navy  officers  score  higher  on  this  test  than  Army 
officers? 


(b) 

By  what 

range 

of  values 

can 

the 

scores 

between 

the  two 

groups  be 

expected  to 

differ 

90£ 

of  the 

time 

• 

(2) 

Solution . 

To 

see 

if  Navy  officers  score 

higher 

on 

the 

exam , 

we 

test 

H0 :  Mx  = 

My  versus 

HI :  Mx  <  My. 

( 3 )  Workspace  Decision  Process. 

(a)  Microcomputer:  Choose  the  Mann- 

Whitney  Test  from  the  main  menu,  and  the  option.  Test 
H0:  Mx  =  My  versus  HI :  Mx  <  My,  from  the  test  menu. 

Skip  to  the  Program  Interaction  section  below. 

(b)  Mainframe:  Enter  MANNWHIT  at  the 

keyboard  and  receive  the  prompts: 

DO  YOU  WISH  TO  COMPARE  THE  MEDIANS  OR 
VARIANCES  OF  THE  POPULATIONS?  ENTER:  1  TO  COMPARE 
MEDIANS;  2  TO  COMPARE  VARIANCES. 

Enter  1.  The  next  prompt  is: 

THE  NULL  HYPOTHESIS  STATES  -  THE 
MEDIANS  OF  X  AND  Y  ARE  EQUAL;  Mx  =  My.  WHICH 
ALTERNATIVE  DO  YOU  WISH  TO  TEST?  ENTER: 

1  FOR  HI :  Mx  <  My;  2  FOR  Mx  >  My;  3  FOR  Mx  #  My . 

Enter  1 . 

(4)  Program  Interaction.  The  prompt  is: 

ENTER  X  DATA  (MORE  THAN  ONE  OBSERVATION  IS 


REQUIRED) . 


Enter  the  X  data  separated  by  spaces.  The 
next  prompt  is: 

ENTER  Y  DATA. 

Enter  the  Y  data.  The  following  is 

displayed. 

THE  SUM  OF  THE  X  RANKS  IS:  224.  THE  U 
STATISTIC  EQUALS:  119. 

THE  P-VALUE  FOR  H0 :  Mx  =My  versus  HI : 

Mx  <  My  IS:  .5078. 

We  do  not  reject  the  hypothesis  of  equal 
population  medians  and  conclude  the  median  of  all  Army 
scores  is  equal  to  the  Navy’s.  The  next  prompt  is: 

WOULD  YOU  LIKE  A  CONFIDENCE  INTERVAL  FOR 
THE  SHIFT  IN  LOCATION  (My  -  Mx)?  (Y/N). 

Enter  Y  ( If  N  is  entered,  the  mainframe 
program  ends;  or,  the  Mann-Whitney  test  menu 
reappears).  The  next  prompt  is: 

ENTER  THE  DESIRED  CONFIDENCE  COEFFICIENT; 
FOR  EXAMPLE:  ENTER  95,  FOR  A  95*  CONFIDENCE  INTERVAL. 

Enter  95.  The  following  is  displayed. 

A  95*  CONFIDENCE  INTERVAL  FOR  THE  SHIFT  IN 
LOCATION  BETWEEN  POPULATIONS  X  AND  Y  IS: 

(  -10  <_  My-Mx  <_  1 0  ). 

The  mainframe  program  ends.  The  microcom¬ 
puter  program  pauses  for  input  from  the  keyboard  by 
prompting: 


PRESS  ENTER  WHEN  READY. 


Press  Enter  and  the  Mann-Whitney  test  menu 

reappears . 

d.  Mann-Whitney  Test  for  Equality  of  Variances 

(1)  Description  of  Problem.  Referring  to 
problem  3  in  section  c(1).  Is  there  sufficient 
evidence  to  claim  that  Army  scores  vary  more  than  Navy 
scores? 

(2)  Solution.  To  see  if  Army  scores  vary 
more,  we  test  H0 :  Vx  =  Vy  versus  HI :  Vx  >  Vy. 

( 3 )  Workspace  Decision  Process. 

(a)  Microcomputer:  Choose  the  Mann- 

Whitney  Test  from  the  main  menu,  and  the  option,  Test 
H0:  Vx  *  Vy  versus  HI:  Vx  >  Vy,  from  the  test  menu  and 
receive  the  prompt: 

ENTER  THE  DIFFERENCE  OF  THE  MEANS  OR 
MEDIANS  (Mx  -  My). 

Because  we  believe  the  population 
medians  to  be  approximately  equal.  We  enter  0.  Skip  to 
the  Program  Interaction  section  below. 

(b)  Mainframe :  Enter  MANNWHIT  at  the 

keyboard  and  receive  the  prompts: 

DO  YOU  WISH  TO  COMPARE  THE  MEDIANS  OR 
VARIANCES  OF  THE  POPULATIONS?  ENTER:  1  TO  COMPARE 
MEDIANS;  2  TO  COMPARE  VARIANCES. 

Enter  2.  The  next  prompt  is: 


•fe  4. a ''t 


THE  TEST  TO  COMPARE  VARIANCES,  REQUIRES 
THE  TWO  POPULATION  MEANS  OR  MEDIANS  TO  BE  EQUAL.  IF 
THEY  DIFFER  BY  A  KNOWN  AMOUNT,  THE  DATA  CAN  BE  ADJUSTED 
BEFORE  APPLYING  THE  TEST.  ENTER  THE  DIFFERENCE  OF 
MEDIANS  (Mx  -  My)  OR  900  TO  QUIT. 

We  enter  0.  The  next  prompt  is: 

THE  NULL  HYPOTHESIS  STATES  -  THE 
VARIANCES  OF  X  AND  Y  ARE  EQUAL;  Vx  =  Vy.  WHICH 
ALTERNATIVE  DO  YOU  WISH  TO  TEST?  ENTER: 

1  FOR  HI :  Vx  <  Vy;  2  FOR  Vx  >  Vy ;  3  FOR  Vx  *  Vy . 

Enter  2. 

(4-)  Program  Interaction.  The  prompt  is: 


REQUIRED)  . 


ENTER  X  DATA  (MORE  THAN  ONE  OBSERVATION  IS 


Enter  the  X  data  separated  by  spaces.  The 


next  prompt  is: 


displayed. 


ENTER  Y  DATA. 

Enter  the  Y  data.  The  following  is 


THE  SUM  OF  THE  X  RANKS  IS:  166.  THE  U 


STATISTIC  EQUALS:  61  . 


THE  P-VALUE  FOR  H0 :  Vx  =  Vy  versus  HI: 
Vx  >  Vy  IS:  .0112. 

Consider  a  significance  level  of  .05. 
Since  a  P-value  of  .0112  is  less  than  .05,  we  reject 
the  null  hypothesis  of  equal  variances  in  favor  of 
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Vx  >  Vy  and  conclude  that  Army  scores  do  vary  more  than 
Navy  scores. 

The  mainframe  program  ends.  The 

microcomputer  program  pauses  for  input  from  the  key¬ 
board  by  prompting: 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  the  Mann-Whitney  test  menu 

reappears . 

e.  Kruskal -Wallis  Test 

(1  )  Description  of  Problem  4.  During  a  recent 
Monster  Mash  involving  four  Navy  SEAL  Teams,  one  of  the 
events  consisted  of  the  number  of  pushups  a  man  could 
do  in  2  minutes.  Eight  men  were  chosen  randomly  from 
each  Team.  The  following  scores  were  recorded. 

SEAL  1:  90  96  102  85  65  77  88  70. 

SEAL  2:  64  79  99  95  87  74  69  97. 

SEAL  3:  101  66  93  89  71  60  76  98. 

SEAL  4:  72  78  73  81  83  92  94  86. 

Are  the  different  Seal  Teams  considered  to 
be  equally  fit? 

(2)  Solution .  To  see  if  the  Seal  Teams  are 
equally  fit.  we  test  the  hypothesis  that  all  the 
population  medians  are  equal . 
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(a)  Microcomputer:  Choose  the  Kruskal- 

Wallis  Test  from  the  main  menu;  and,  once  the  test  menu 
is  displayed,  press  Enter. 

(b)  Mainframe:  Enter  KRUSKAL  at  the 

keyboard. 

(4)  Program  Interaction.  The  prompt  is: 

ENTER  THE  NUMBER  OF  POPULATIONS  TO  BE 
COMPARED  (MUST  BE  GREATER  THAN  TWO). 

Enter  4.  The  next  prompt  is: 

ENTER  YOUR  FIRST  SAMPLE. 

Enter  the  SEAL  1  data  separated  by  spaces. 
The  next  prompt  is: 

ENTER  YOUR  NEXT  SAMPLE. 

Enter  the  SEAL  2  data.  The  next  prompt  is: 
ENTER  YOUR  NEXT  SAMPLE. 

Enter  the  SEAL  3  data.  The  next  prompt  is: 
ENTER  YOUR  LAST  SAMPLE. 

Enter  the  SEAL  4  data.  The  following  is 

displayed. 

THE  H  STATISTIC  EQUALS:  .1335. 

THE  P-VALUE  FOR  H0 :  THE  POPULATION  MEDIANS 
ARE  EQUAL  versus  HI :  AT  LEAST  TWO  POPULATION  MEDIANS 


ARE  NOT  EQUAL  IS: 


.98893. 


We  do  not  reject  the  null  hypothesis  that 
the  population  medians  are  equal  and  conclude  that  the 
SEAL  Teams  are  equally  fit. 

The  mainframe  program  ends.  The  microcom¬ 
puter  program  pauses  for  input  from  the  keyboard  by 
prompting : 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  the  Kruskal-Wal 1 is  test 
menu  reappears. 

f .  Kendall ’ s  B 

(1)  Description  of  Problem  5.  In  order  to 
determine  if  cold  weather  affects  target  marksmanship. 
Naval  Special  Warfare  recorded  small  arms  marksmanship 
scores  and  corresponding  air  temperatures  for  a  period 
of  one  year.  20  men  were  chosen  at  random,  and  their 
scores  averaged  for  different  air  temperatures.  The 
average  score  for  each  air  temperature  is  shown  below. 

Air  temperature  (X):  50  55  20  50  65  55  30 

52  40  60. 

Average  scores  (Y):  210  200  165  165  260 

215  175  191  180  235. 

Can  it  be  said  that  colder  temperatures 
have  an  ei'fect  on  marksmanship  scores?  Is  that  effect 
positive  or  negative? 


(2)  Solution .  We  test  the  null  hypothesis 
that  no  association  exists  between  cold  temperatures 
and  marksmanship. 

( 3 )  Workspace  Decision  Process. 

(a)  Microcomputer:  Choose  Kendall’s  B 

Test  from  the  main  menu;  and,  once  the  test  menu  is 
displayed,  press  Enter. 

(b)  Mainframe:  Enter  KENDALL  at  the 

keyboard. 

(4)  Program  Interaction.  The  prompt  is: 

ENTER  X  DATA  (MORE  THAN  TWO  OBSERVATIONS 

ARE  REQUIRED). 

Enter  the  X  data  separated  by  spaces.  The 
next  prompt  is: 

ENTER  Y  DATA  (NUMBER  OF  Y  ENTRIES  MUST 
EQUAL  NUMBER  OF  X  ENTRIES). 

Enter  the  Y  data.  The  following  is 

displayed . 

KENDALL’S  B  EQUALS:  .7817. 

THE  P-VALUE  FOR  H0 :  NO  ASSOCIATION  EXISTS 
versus:  HI:  DIRECT  ASSOCIATION  EXISTS  IS:  .00045. 

THE  P-VALUE  FOR  THE  TWO-SIDED  TEST  OF 
HYPOTHESIS  IS:  .0009. 

Since  the  P-value  for  the  one-sided  test 
equals  .00045,  we  reject  the  null  hypothesis  that  no 


association 


exists 


between 


temperatures 


and 


marksmanship  in  favor  of  direct  association.  We 
conclude  that  colder  temperatures  tend  to  cause  lower 
marksmanship  scores. 


The  mainframe  program  ends.  The  microcom¬ 
puter  program  pauses  for  input  from  the  keyboard  by 
prompting: 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  Kendall’s  B  test  menu 

reappears . 

g.  Spearman’s  R 

( 1 )  Description  of  Problem  6.  When  fitness 
reports  are  written,  officers  of  the  same  grade  are 
ranked  against  each  other  based  upon  their  demonstrated 
level  of  performance.  Last  marking  period,  the 
Commanding  and  Executive  Officers  separately  ranked  9 
Ensigns  as  shown  below. 

Ensigns 

ABCDEFGHI 

CO  ( X ) :  641  528379 

XO  (Y):  56341  9728 

Does  any  association  exist  between  the  two 
sets  of  rankings? 

(2)  Solution .  We  test  the  null  hypothesis 
that  no  association  exists. 
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(a)  Microcomputer:  Choose  Spearman’s  R 

Test  from  the  main  menu;  and,  once  the  test  menu  is 


displayed,  press  Enter. 

(b)  Mainframe:  Enter  SPEARMAN  at  the 

keyboard. 

(4)  Program  Interaction.  The  prompt  is: 

ENTER  X  DATA  (MORE  THAN  TWO  OBSERVATIONS 
ARE  REQUIRED). 

Enter  the  X  data  separated  by  spaces.  The 
next  prompt  is: 

ENTER  Y  DATA  (NUMBER  OF  Y  ENTRIES  MUST 
EQUAL  NUMBER  OF  X  ENTRIES). 

Enter  the  Y  data.  The  following  is 

displayed. 

SPEARMAN’S  R  EQUALS:  .5500. 

THE  P-VALUE  FOR  H0 :  NO  ASSOCIATION  EXISTS 
versus:  HI:  DIRECT  ASSOCIATION  EXISTS  IS:  .0664. 

THE  P-VALUE  FOR  THE  TWO-SIDED  TEST  OF 
HYPOTHESIS  IS:  .1528. 

Consider  a  significance  level  of  .05. 
Since  a  P-value  of  .0664  exceeds  .05,  we  do  not  reject 
the  null  hypothesis  that  no  correspondence  exists 
between  the  two  sets  of  rankings. 
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The  mainframe  program  ends.  The  microcom¬ 
puter  program  pauses  for  input  from  the  keyboard  by 
prompting: 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  the  Spearman’s  R  test  menu 


reappears . 

h.  Nonparametric  Simple  Linear  Regression;  Least 
Squares 

(1)  Description  of  Problem  7.  Battery-powered 
Swimmer  Proplusion  Units  are  sometimes  used  to  aide 
swimmers  during  long  underwater  swims.  Recent  tests 
have  shown  that  a  nearly  linear  relationship  exists 
between  water  temperature  and  battery  life  for  these 
units.  The  following  17  data  points  were  randomly 
selected  from  the  test  results. 


Water  temperature 

(X) 

Battery  life  (Y) 

70 

3 

65 

2.75 

50 

1  .8 

40 

1  .2 

60 

2.4 

55 

1  .9 

52 

1  .75 

50 

1  .7 

43 

1  .6 

40 

1  .  1 

72 

2.75 

55 

2 

48 

1  .5 

35 

.9 

70 

3.3 

68 

3 

57 

2.3 

(a)  Find  the  fitted  regression  equation. 

i 


(b)  For  the  following  water  temperatures, 
predict  the  battery  life  of  the  units:  61  52  46  36. 

(c)  Can  we  determine  with  any  certainty 
if  the  slope  of  the  regression  line  equals  .05. 

(d)  What  range  of  values  could  be  used 
as  the  slope  of  the  estimated  equation  line  9 of  the 
time? 

(2)  Solution.  To  determine  the  estimated 
regression  equation,  we  use  nonparametr ic  lirear 
regression . 

( 3 )  Workspace  decision  process. 

(a)  Microcomputer:  Choose  Nonparametr ic 
Simple  Linear  Regression  from  the  main  menu;  and,  once 
the  test  menu  is  displayed,  press  Enter. 

(b)  Mainframe:  Enter  NPSLR  at  the 

keyboard. 

(4)  Program  Interaction.  The  prompt  is: 

ENTER  X  DATA  (MORE  THAN  TWO  OBSERVATIONS 

ARE  REQUIRED). 

Enter  the  X  data  separated  by  spaces.  The 
next  prompt  is: 

ENTER  Y  DATA  (NUMBER  OF  Y  ENTRIES  MUST 
EQUAL  NUMBER  OF  X  ENTRIES). 

Enter  the  Y  data.  The  following  is 


displayed. 


THE  LEAST  SQUARES  ESTIMATED  REGRESSION 


EQUATION  IS: 


Y  -  -1 .263  +  . 060668X . 

The  next  prompt  is: 

DO  YOU  WISH  TO  ENTER  SOME  X  VALUES  TO  GET 
THE  PREDICTED  Y » S?  (Y/N). 

Enter  Y  ( If  N  is  entered,  the  program  skips 
to  hypothesis  testing  for  the  slope). 

ENTER  X  VALUES. 

Enter  61  52  46  36.  The  next  prompt  is: 


THE  PREDICTED  Y  VALUES  ARE: 


2.44  1 .89 


1.53  .92 


WOULD  YOU  LIKE  TO  RUN  SOME  MORE  X  VALUES? 
Enter  N.  The  next  prompt  is: 

WOULD  YOU  LIKE  TO  TEST  HYPOTHESIS  ON  B,  THE 
SLOPE  OF  THE  EQUATION?  (Y/N). 

Enter  Y  (If  N  is  entered,  the  program  skips 
to  confidence  interval  estimation).  The  next  prompt 


ENTER  THE  HYPOTHESIZED  SLOPE. 

Enter  .05.  The  following  is  displayed. 
SPEARMAN’S  R  EQUALS:  .5756. 


THE  P-VALUE  FOR  H0 :  B  =  .05  versus  HI :  B  > 


.05  IS: 


0079 


THE  P-VALUE  FOR  THE  TWO-SIDED  TEST  OF 


HYPOTHESIS  IS:  .0158. 


I 
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Consider  a  significance  level  of  .05.  Since 
a  P-value  of. 0079  is  less  than  .05,  we  reject  the  null 
hypothesis  that  B  »  .05  in  favor  of  B  >  .05,  and 
conclude  the  slope  of  the  regression  line  is  greater 
than  .05.  The  next  prompt  is: 

WOULD  YOU  LIKE  A  CONFIDENCE  INTERVAL  FOR 
THE  SLOPE?  (Y/N). 

Enter  Y  (If  N  is  entered,  the  mainframe 
program  ends;  or,  the  Nonparametric  Regression  test 
menu  reappears).  The  next  prompt  is: 

ENTER  THE  DESIRED  CONFIDENCE  COEFFICIENT; 
FOR  EXAMPLE:  ENTER  95,  FOR  A  95%  CONFIDENCE  INTERVAL. 

Enter  90.  The  following  is  displayed. 

A  90%  CONFIDENCE  INTERVAL  FOR  B,  THE  SLOPE 
OF  THE  ESTIMATED  REGRESSION  LINE,  IS: 

(  .05333  <.  B  <_  .07  ). 

If  the  estimated  slope,  does  not  lie  within 
the  confidence  interval,  the  following  would  be 
displayed. 

THE  LEAST  SQUARES  ESTIMATOR  OF  B  LIES 
OUTSIDE  THE  CONFIDENCE  INTERVAL.  DISCARD  THE  LEAST 
SQUARES  EQUATION  AND  USE: 

Y  =  -1 .4458  +  .060833X. 

THIS  EQUATION  IS  BASED  ON  THE  MEDIANS  OF 
THE  X  AND  Y  DATA,  AND  THE  MEDIAN  OF  THE  TWO-POINT 
SLOPES  CALCULATED  FOR  THE  CONFIDENCE  INTERVAL  ON  B. 


The  next  prompt  is: 

DO  YOU  WISH  TO  ENTER  SOME  X  VALUES  TO  GET 
THE  PREDICTED  Y’S  FROM  THE  NEW  EQUATION?  (Y/N). 

To  compare  results,  let  us  input  the 
temperatures  in  the  new  equation.  Enter  Y  (If  N  is 
entered,  the  mainframe  program  ends;  or,  the 
Nonparametric  Regression  test  menu  reappears).  The 
next  prompt  is: 

ENTER  X  VALUES. 

Enter  61  52  46  36.  The  following  is 

displayed. 

THE  PREDICTED  Y  VALUES  ARE:  2.265  1.72 


1.35  .74. 

The  next  prompt  is: 

WOULD  YOU  LIKE  TO  RUN  SOME  MORE  X  VALUES? 
Enter  N.  The  next  prompt  is: 

WOULD  YOU  LIKE  TO  TEST  HYPOTHESIS  ON  B,  THE 
SLOPE  OF  THE  EQUATION?  (Y/N). 

To  compare  results  once  again,  we  enter  Y 
(If  N  is  entered,  the  mainframe  program  ends;  or,  the 
Nonparametric  Regression  test  menu  reappears).  The 
next  prompt  is: 

ENTER  THE  HYPOTHESIZED  SLOPE. 

Enter  .05.  The  following  is  displayed. 
SPEARMAN’S  R  EQUALS:  .8287. 
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THE  P-VALUE  FOR  H0 :  B  =  .05  versus  HI :  B  > 
.05  IS:  .0000. 

THE  P-VALUE  FOR  THE  TWO-SIDED  TEST  OF 
HYPOTHESIS  IS:  .0000. 

The  mainframe  program  ends.  The  microcom¬ 
puter  program  pauses  for  input  from. the  keyboard  by 
prompting: 

PRESS  ENTER  WHEN  READY. 

Press  Enter  and  the  Nonparametric 
regression  test  menu  reappears. 
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APPENDIX  D 


MAIN  PROGRAM  LISTINGS  FOR  MICROCOMPUTER  WORKSPACE 


I 

k 

9. 

10, 

II 
12 
13 

rfr 

pi!: 

20, 

21 

22 

23] 

.2“, 

25, 

26] 

27 

28 
29, 

30 

31 

32 
.33, 
34, 

33 

36 

37, 

38 

39 
>0, 
>1. 

:51: 

;*ts; 

48 

49, 

50 

51 

52 

53, 

54, 

55 

56 

57 

:il: 

85 

:665 

66 

82 

:?o] 

'M 

[73' 


ORDER  X  IN  INCREASING  ORDER 
COMPOTE  CORRENT  RANKING  OF  I 
NOW  ORDER  7  RANKS  IN  INCREASING  ORDER 


7  KENiAiAAiBiBXiBIiCiCXtCI-.DiDDiDX-.DIiDXl-.SiPOSiNEGiXXilliNiDEN-.NNiNOM-.P ; 

B  PTBibHbhcii6N 5  cdMpbfhs !  ihf  kENDkliH '  Statistic  uaica  is  a  measure 

i  iitmxipjm  jgysggf 

a  SUBPROGRAMS  CALLED  BI  TRIS  FUNCTION  INCLUDE'.  TIES,  TIESK,  KENDAL? , 
a  INTERP,  INPUT  AND  NORMCDF . 

a  D I SPLAT  TEST  MENU  AND  INPUT  DATA 

Nl'.E+MENO  KENQSi 
*<S=1)/31 

MENU  MAINQSi 

BliR*INPOT  2 

m9^R 

^T*lQ+Z)rR  QRDER  r  IIf  INCREASING  ORDER  OF  X 
A+llt>Xl 

a3*x cax: 

n 
« 

aD*ALjl  DIES  EXIST  IN  EITHER  X  OR  T  RANKED  VECTOR  USE  MID -RANK  METHOD 
DD* 1  TIES  D 
yy+\  TTES  3 

B  FIND  ORIGINAL  RANKING  OF  I  WITH  TIES  RESOLVED 

TT+DDICl 

ftN COMPUTE  NUMBER  OF  DISTINGUISHABLE  PAIRS  ' 

NN+(NxlN-l))*2 
S*  pO 
AA*0 

a  _ 

LBxZ1xxTaA] > (AA*XX )) 

CX-KXX[AA]<UA+XX))*<  1) 

DX*BX*CX  ,  ,  ... . 

fly-*-<TJ[AA]>(AA+JT)}  . 

cr*aruA]<(AA+iT))*(  1) 

D7+B7+CI 
DX7*DX*D7  % 

P0S*(DX7> 0)  , 

WEOU}XX<0)*(  1) 

S-S.POS  ,w?<; 

b-(AA<(W-1))/C1  vector  to  determine  s 

C+a, /C 

B  OBTAIN  THE  NUMBER  OF  TIES  IN  EACH  VECTOR  US INC  THE  TIESK  FUNCTION 
U*TIESK  B 
V*TIESK  D 
SU++/l2’.U\ 

SV++ / ( 9  i  V ) 

*  CALCULATE  TBR  B  STATISTIC  INCLUDING  THE  CORRECTION  FOR  TIES 

t+s*1<.nn-su)*Xnn-sv))*o.s 

AT*\T 

aN>13)/NCALL  KENDALP  TO  CALCULATE  THE  RIGHT  TAIL  OF  THE  CDF  OF  B 

3P*ALID*'NTER?  TO  CALCULATE  P -VALUE  31  INTERPOLATION 
PVAL+AT  INTERP  ? 

~<PVAL*  U/L 3 

PVAL* 0 . 5 

n* CALCULATE  P  VALUE  USING  NORMAL  APPROX. 

NORMlNUM+T.3*AT)*{  {2*NN)*Q .  5) 

MM2X^  +  5”*°-5 

PVAl+1- {NORMCDF  Z) 

a  IF  B  IS  POSITIVE  PRINT  OUT  DIRECT  ASSOCIATION. 

L3i*(.T>0l/Li 
CB A* 'INDIRECT' 


Positive  ones  come  from  a  runs  up  condition s  necattve  i  from  runs  down 

ZERO  IS  SCORED  FOR  TIES  .  MULTIPLE  THE  RESULTS  FOR  EVERT  ELEMENT  AND  SUM 


78 


|7ui  *n 

75  Li :CBA*' DIRECT' 

76  LT.PV+2*PVAL 
77  +(P7Sl)/£8 

78  PV+1 

79 !  LB:' KENDALL" S  B  EQUALS:  ' . (u«r) .QTCNL 

80  'TBS  P-VALUE  POR  SO:  NO  A iSOCIAtlON  EXISTS  VERSUS ' 

81  1  Hi:  '.(*CBA),'  ASSOCIATION  EXISTS  IS: 

82  'TBS  P-7ALUS  POR  TBE  Ti/O-SIDBD  TEST  OP  BIPOTBESIS  IS: 
83  ' PRESS  ENTER  WBEN  READS. ' 

34,  ww+a 
[853  *N 1 
7 


'  .(.UvPVAD.OTCNL 
, l*9PV) .UTCNL 


7  KRUL :NUM: DENOM: A:C;B:D:K;AA:BB:DDiE:F:N:OP:P:PVAL:R;SOFR:SR; TSORxCBA : B 
A  TBlS  PbNCTIofi  COMPfjlkS  TBE  XUSXAL-W ALLIS  TEST  STATISTIC  3  WHICH  IS 
a  A  MEASURE  OP  TBE  EQUALITX  OP  X  INDEPENDENT  SAMPLES, 
a  SUBPROCRAMS  CALLED  BI  THIS  FUNCTION  INCLUDE:  TIES,  TIESK,  INDEXPLS , 
a  FDISTN.  INTER P  AND  THE  VARIABLES  PMATKW20 .  PMATKWS 1.  PMATXW 33,  PMATKN 34 
a  PMATXW'iX.  PMAT.KM 42.  AND  PMATXW 43. 

a  MENU  CHOICES  AND  ROUTE  TO  PROPER  STATEMET  POR  ACTION. 

Nl:B*MENU  KRULQH 
♦(S=l )/Hl 
MENU  MAINQBi 

•*0 

3PP*u 

31 : ' ENTER  TBE  NUMBER  OP  POPULATIONS  TO  BE  COMPARED  ( MUST  BE  CHEATER  THAN  Z 
WO).  ' 

♦Tor<3  )vaitJO*0))/Hl 
*■(  lpJT)>l  )/E  1 

A  INITIALIZE  VECTORS  E  AND  P  AND  VARIABLE  C 

B+P+SOFH+pO 
OO 

A THIS  LOOP  FACILITATES  ENTERING  TBE  SAMPLE  VECTORS  AND  STORING  THEM 
ISA*' FIRST' 

31:001 

’ENTER  JOUR  '  ,  (.nCBA  ) , '  SAMPLE. ' 

♦7? Oo 3  )*0 )/NEXT 
OlpH 

a  CONCATENATE  SAMPLES  AS  TEST  ARE  ENTERED  AND  STORE  THEM  IN  VECTOR  E 
NEXT  *  B+E  D 

a  '  ’RECORD  TBE  LENGTHS  OP  THE  SAMPLES  AS  TBEI  ARE  ENTERED 
P+F.pD 
CBA*'NEXT' 

*\c< (H-ljO/il 

CBA*'LAST' 

~lC<K)/Ll 

a  RECORD  SIZE  OP  ALL  SAMPLES  WBEN  COMBINED 

N*+/P 

A  ORDER  SAMPLE  SIZES  LARGEST  TO  SMALLEST 
0F*PZ9P 3 

a  ORDER  COMBINED  SAMPLE  VECTOR  TO  BE  USED  BI  TIES  FUNCTION 

* INDEXPLS  TO  INCREMENT  INDEXES  WBEN  TIES  OCCUR  WITHIN  ONE  SAMPLE 
AA*P  INDEXPLS  E 

A  CALL  TIES  TO  BREAK  TIES  BI  MIDRANK  METBOD 

BB*  1  TIES  D 
OO 

a  THIS  LOOP  CALCULATES  TBE  B  STATISTIC 

£2  :OC+l 

l>SR*+/BBi(PiCl1t(AAiC^l^-]F0R  ZACB  SAHFLS  IS  CALCULATED 

A  CALCULATE  SUM  OP  kANKS  SQUARED  DIVIDED  BI  TBE  INDIVIDUAL  SAMPLE  SIZE 
SR+(SR*2)*PZC\ _ 


SOFR+^OPR i SR 


TORE  EACH  CALCULATION 


A  SUM  ACROSS  ALL  SAMPLES 

TS0R++/S0FR 

a  CALCULATE  FINAL  3  STATISTIC 

3~<TSORx(12*(Nx(N+l  )  )  )  -  ( 3  *  GV+l ;  ) 
a  RECALCULATE  B  WITH  CORRECTION  FOR  TIES 
A*TIESK  E 

NUM+(  +  /(.A*3)')-(.+/A) 


DEN0M*Nx((N* 2)-l) 

B*H* (1- (NUM* DENOM)) 

, SI STEM  QP  LOGICAL  STATEMENTS  ENSURE  PROPER  PROB .  IS  ACCESSED 


)/P APPROX 


+(0Fll]<2}/0UTPUT 
*(K>*)/P APPROX 
*\K=2)/IF 

3  3  3  5))''(OPC13>3) 
*U/(0P*  2  1  1  1))/00TPUT 
♦(  (H=4  )a«7PC13  =  2  )  J/P42 


C7i] 

[72 

73 

74 

75 

76 

77 

78 

79 

80 
,31 
82 

83 

84 

85 

86 

87 

88 

89 

90 
,91 
92 
33 
,94 

95 

96 

97 

98 

99., 

;  1003 
101 
[102] 
.  103  J 

[104] 

[105 

[106 

[l07 

[108] 


*}pui(°FS  31  1  i;)v(A/(<7fs  321  1>)V<*/«?P=  3  3 
T!2^C1] >4 )/FAPPROX 

v/(OP=  2  1  i})v(a/(0P=  3  1  1))) /OUTPUT 
)^C11=4)a(0^[?]=1)5/?31  „  ,  , 

lOF=  3  3  ' 

[l]=4)>/ 

PROPRIA! 


.  1  1))v(a /(OF-  3321)) 

..P41 

IF : * ( ol { 1] >4 ) /FARRROX 

♦  (  (  A  / 

fQ'pf  _ 

*SS' /)0>='3  '2'l})vU/l'OF=~3  3  1))v(a /{OF*  3  3  2))v(a  /(OP=  3  3  3)))/P33 
♦<(a/(OP=  3  2  2)jv<0pm=4))/P34 
#  CAIX  APpSOPSXATS  7>LRJASI2  ACCESS  CDF 

P23:P+P*ArXV20f  • 

-PH 

P3 1 : ?*PMATKW3 1 C 08- 5  ) ; ; ] 

^  p  ||^ 

P33:P+PMATKV33ZIN-S): ;] 

*PH 

P3'i:P+PMATKW2'*ZZN-S)i ;] 

♦PM 

Pm  lP*PHATKWm  [  08-5  ): ;  ] 

+PIH 

P42 :P*PMATKWm  C  08-5  ) : :  ] 

P43  :  P+PMATKW* 3  C  (38- 7  ) :  :  ] 

fl  CALL  INTERP  TO  bXLCULATE  P -VALUE  ST  INTERPOLATION 
?M:PVAL-*-3  INTERP  ? 

♦(P7A£=  D/OUTPUT 

^^ALCULATR  .P-VALUE  .USING  IRE  F .DIRT  U/ONE  LESS  D.F.IN  DENOM  APPROX. 

P7AW' 

*C  5 

OUTPUT \PVAL*' GREATER  T3AN  .25' 

Li: 'TBS  3  STATISTIC  EQUALS: 

'TEE  P -VALUE  FC~  - 

'  31:  AT 

0 TCNL 

' PRESS  ENTER  WREN  READ T . ' 

AA*0 

♦/VI 

Si: 'ERROR:  ENTER  A  SINGLE  INTECER  VALUE  GREATER  TRAN  2;  TRT  AGAIN. 1 ,CTCNL 
♦SI 
7 


CALCULATE  P-VALUE  USING  TEE  F  DIRT 
?PROX:F~UN-K)xB)*(CK-l)xJjN-l)-B)l 
VAL*l-(.  (  vS-1 )  ,  (.(.N-k)-!} )  FDISTN  F) 

7REATER  TRAN  .25' 
riSIIC  EQUALS:  ',(4 *R).2TCNl 

FOR  30 :  T3E  POPULATION  MEDIANS  ARE  EQUAL  VERSUS  >  , TTCNL 
?  LEAST  TWO  POPULATION  MEDIANS  ARE  NOT  EQUAL  -I*  ' .X^vPVAL), 


NCI 


W’N:M;PV2:A;B:C;G:MM}NN;RX;U: NM1 ; P:NU ; PVAL :NM : NUMZ ;  J8CWZ 1 ; DEN: DENC ; DE 
•JciTCUNOM-.  Z ,  21  dLPkA ;  CDF i INDEX ilPx :CX :  UALprX  ;BB  ;CC;U1  ;l/2  -.Pp-.NNl  ;NN 


:io 

13 

,14 

15 

16 

17 

18 
.19 
20 
,21 
,22 
.23 
,24 

25 

26 
,27 

! 

in 

L  3  2 
,33 
,34 
,35 
36 
,37 
.38 

39 

40 
.41 
.42 
>3 
44 


C.I.  FOR  {NT -MX ) ,  THE  SR I FT  IN  LOCATION  __ 
LIED  BY  TRIE  FUN bT ION  INCLUDE :  TIES.  TIES 2 


B  WBITNl.  _  _ 

B  OF  LOCATION  AND  SCALE.  TRE 

B  ALSO  COMPUTED .  SUBPROGRAMS  CA _ _ _  _ 

tk^INDEXPLS,  VARMN ,  HANWP ,  INPUT.  CONFMW,  NORMCDF,  AND  NORMPTR. 

N1 : MENU  MANWRELP 

B  MENU  CBOICES)  AND  ROUTE  TO  PROPER  STATEMENTS  FOR  ACTION. 

S3 : D+CBQICEM  PACEDMENU  MANWQRi 


CALCULATE  TRE  MANN- 
~  VALUE  FOR  TRE  TEST 
- - JS 


DIFF* 0 
♦(0= l)/JVl 


MENU  MAINQH 

S3: ’ENTER  TRE  DIFFERENCE  OF  TEE  MEANS  OR  MEDIANS  (W£  -Ml). 
DIFF* 0 

♦<(pCXPP)>l)/S3 

fl  _ _  _  „  ,  ENTER  DATA  VECTORS 

SI  '.RENTER  X  DATA  ( MORE  TRAN  ONE  OBSERVATION  IS  REQUIRED).1 


N*D 

1  ENTER  T  DATA.' 


H+h 

* Nlf -CALCULATIONS  INVOLVE  VARIANCES  ADJUST  X  BT  TRE  DIFFERENCE  IN  MEANS 

fl  CONCATENATE  X  AND  I  SAMPLE  VECTORS 

A+N  ,M 

o  DETERMINE  SIZE  OF  X  AND  I  VECTORS  AND  ASSIGN  TO  NN  AND  MM 

NN+oN 
HM*oM 

*  COMPUTE  SIZE  LIMIT  OF  LEFT  TAIL  OF  NULL  DISTRIBUTION 

NM+  \rrrtxMM  )  ♦  2 
NM1+INM 

fl„  ORDER  A  AND  ASSIGN  TO  B 

B+A  [  aA  ] 

C+(NN,MM)  INDEXPLS  A 

«  „  mrQALL  TIES  FUNCTION  TO  BREAK  TIES  USING  MI DRANK  METRO D 

G* 1  TIES  B 

fl  IF  FALSE  CALCULATE  TEST  FOR  VARIANCES 
♦(£=2,3  .4,5  )/S5 

« „  -A^pL  VARm/  TO  GENERATE  RANKS  REQUIRED  FOR  VARIANCE  TEST 
GG+VAJiMW  ("ff+Wf ) 

fl  CALL  TIES  TO  RECORD  TIES  IN  TRE  DATA  AND  BREAK  TIES  IN  GO 
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SH 


$ 


Si 


a0*01 7  TJ?5  5  CA0CUMIP  SUM  OF  X  RANKS 

ss.fx*+/(oC(cci,  a ju/nubit  a  statistic 

0*RX-((NNx(NN+l))+2) 

ui*u 

a  IF  SIZE  OF  X  TIMES  SIZE  OF  I  >  80 ;  GO  TO  NORMAL  APPROX 

♦  <OVMx2>80)/£2 
NN1+L /NN , MM 
NN2*\ /NN, MM 

aMANWP  FUNCTION  CALCULATES  LEFT  TAIL  CUMULATIVE  PROBS .  OF  U  STATISTIC 
P*NN2  MANWP  NN1 
♦(0  =  5  )/L 10 

a  LOGICAL  STATEMENT  ENSURES  ONLI  LEFT  SIDE  OF  NULL  DIST  IS  USED 

♦  (US/VM1 )/03 

a  CONVERT  U  STAT  WEEN  CREATER  THAN  LEFT  TAIL  VALUES 

U*(NNxMM)-U 

a  IF  C/1  IS  A  FRACTIONAL,  INTERPOLATE  P  VALUE 
L3 :♦( (1 1 U)=0)  /NON 
U2*‘.U 
♦<U2>0 )/Pl 
PV^l-  ( (PCU2+1]  )*2) 

♦P3 

?l:?T/*l-(  (PCU2]  T?rU2-Ll2  )  +  2  ) 

P3  :PVI*(.PiU2+l}  +PZU2+22  )  +  2 
+CBECK  % 

NON:*lU>Q)/GO 
PV*  1 
♦P2 

GO:PV*l-prui 
P2:PVI*?r>u*\ n 
CHECK: *(UlZNMl) / 04 

OT/O^PV 
?V—?VI 
PVI*?V 2 
♦04 

a  IF  CONFIDENCE  INTERVAL  DESIRED  GO  TO  010 
L2:*(D=z )/0lO 

a  COMPUTE  TEE  NORMAL  APPROXIMATION  W/CORRECTION  FACTOR 
NUMZ*(U+0 . 5 )-NM 
NUMZl*(U-0.3)-NM 
DEN~i(MM*NNx  (MM+NN+ 1 )  )  +  12  )*0 . 3 
Z*NUMZ+DEN 
Z 1 *NUMZ 1 ♦ DEN 
NUM* I (U-NM ) 

nnii/'  ^  /  t  f  it  if  -  -  * *  •  -  .  .  •  * 


DENC*(  ( ( (NN+MM-l  )x  (DEN *2  )  )♦  (NN+MM-2  ) )-  (  (  (MUM-0 . 5  )*2  )+  (NN+MM-2  ) )  )*0 . 5 
DENCI+UUnN+MM-I  )x  (VEN*  2  )  )♦  (NN+MM-2  )  )-  (  ( (NUM+ 0 . 5  )*2  )♦  (NN+MM-2  )  )  )*0 . 5 
TC+\NUn“vj  •  5  j+DENC 
TCl*(NUM+0 . 5 )+DENC 1 
*{USNM) /SECOND 

Pvl*( (NORMCDF  Zl+( (NN+MM-2)  TDISTN  TC))+2 
PV* (11- (NORMCDF  Zl))+(1- ( (NN+MM-2 )  TDISTN 

♦04 


SIN  rcD))  +  2 


SECOND : PVI*( ( NORMCDF  Z )+<l- ( (NN+MM-2 )  TDISTN  TC)))+ 2 
,  PV*l(l- (NORMCDF  Zl))+( (NN+MM-2)  TDISTN  TCl))+2 
L*:PV3*2x(\/(PV ,PVI)) 

♦(P73S1 )/W5 


M5  :PVM*(3  ,l)p  (P7I,P7.P73  ) 

■rPP  SUM  OP  TEE  X  RANKS  IS:  1 , (»FX) , ' .  TUP  U  STAIISHC  EQUALS:  ' ,(®U1), 
QTCWS 

a  LOGICAL  STATEMENT  FOR  VARIANCE  OUTPUT 

*(D=B,7,8}/VAR 

(6  4  iPPM[0?Vi])ffarc^  =  H1  255525  S1:  '’^LOCICZD- lsl]).'  MJ  IS: 

PAP?PPM*(3,1 )p (PP,PVT,P73 ) 

(6H?  SpVM[0-5-lf  =  71  ffl:  ^2  '  ,  (vLOCIClD-Sill  ),  '  7Z  IS: 

'PPPSS  FWIPP  &3P/V’pP>U)r.  ' 

BB+a 

♦M3 

08: ^WOUIU^IOU  £7FP^  CONFIDENCE  INTERVAL  FOR  THE  SHIFT  IN  LOCATION  (MI  - 
SS*0 

-<33=  'V  VZV3 
013 : CC* INPUT  z 
ALPHA*1 10C-'U:+23C 

a  ROUTE  TO  NORMAL  APPROX.  FOR  CONF .  INT .  OF  LARGER  SAMPLE  SIZES 
♦( (NMx2 )>8O}/05 

a  COMPUTING  CONFIDENCE  INTERVALS  BI  EXACT  P -VALUE 

CDF+P 

a  INDEX  POSITION  OF  VALUE  IN  CDF  S  ALPHA 

INDEX* (+/( CDFSALPHA ) ) 

*(INDEX>0 )/L6 
INDEX* 1 
♦06 

a  COMPUTING  CONFIDENCE  INTERVALS  USING  NORMAL  APPROX.  W/C.F. 

05 :DEN*( (MMxNNx (MM+NN+1 ))  +  12  )*0 . 5 
UALPHA*(DENx (NORMPTH  ALPHA)  )+NM-0. 5 
a  ROUND  UALPHA  DOWN  AND  INCREMENT  SI  O/VP 


.IWI 


v’v'v'v' 
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INDEX*iUALPHA*l 
16 : IPX*NN. INDEX 

CI^PX^CONFMW  A  confidencb  interval  for  the  shift  IN  ' 

'LOCATION  BETWEEN  POPULATIONS  X  AND  I  IS'.'.QTCNL 
<  (  ' .(JCIC1] S  HI  •  Hi  S  '.(•CXC23),*  )',Q ICNL 

'PRESS  ENTER  WHEN  READ Y. ' 

Bfl-Q  0 

*N  3 

El '.'ERROR:  SAMPLE  CONTAINS  LESS  THAN  TWO  ENTRIES’,  TRI  AGAIN .' .UICNL 
* 31 

E3:' ERROR:  YOU  HAVE  ENTERED  MORE  THAN  ONE  VALUE.  TRY  AGAIN .' .UTCNL 
*B  3 
7 


7  NPLR: N :  SUMX :SUMY  :XBAR  :YBAR :  SUMX2  :SUMXY:B ;  A  :WW  iXX-.BB  ;U :D ;  ALPHA  :P:CC\  CDF ; 
TALPHA  :  NN :  cl :  SLOPES ;  P.R :  SP.  :YY :  DENOM : INDEX  :FF:X:Y;q;R;  CHA  ;E-tPV 
a  PROGRAM  CONDUCTS  NONPARAMETP.IC  LINEAR  REGRESSION.  IRE  LEAST  SQUARES 

a  ESTIMATED  REGRESSION  LINE  IS  COMPUTED  WITH  HYPOTHESIS  TESTINC  AND 
a  CONFIDENCE  INTERVAL  AVAILABLE  FOR  THE  SLOPE  B .  IF  B  DOES  NOT  LIE  IN  THE 
a  C.I.  AN  ALTERNATE  REGRESSION  LINE  IS  PROPOSED.  SUBPROGRAMS  CALLED  ARE 
a  SPMANP,  ZEND ALP,  NORMPTH ,  INPUT,  AND  CONFLR. 

QPP*5 

a  DISPLAY  MENU  AND  INPUT  DATA. 

Nl:E*MENU  NPLRQSi 
*(E-1] /N2 
MENU  MAI NO M 

N2:H*INPUT  2 
Q*l-R 

I-KQ+D+H 

a  ASSIGN  THE  SIZE  OF  X  ( AND  Y)  TO  N 

N*oX 

a  COMPOTE  THE  SUM  OF  X'S  AND  Y'S 

SUMX**/X 
SUMY++/Y 

a  COMPUTE  THE  MEAN  OP  X  AND  Y 

XBAR*SUMX*N 
YBAR*SUMY*N 

a  COMPUTE  THE  SUM  OF  THE  X'S  SQUARED 

SUMX2++/ (X*2 ) 

a  COMPUTE  THE  SUM  OF  X  TIMES  Y 

SUMXY*+/(X*Y ) 

a  COMPUTE  <B'.  THE  SLOPE  OF  THE  ESTIMATED  LEAST  SQUARES  REGRESSION  LINE 
B*(  (.NxSUMXD-lsUMXxSUMY))*  (.(N*S0MX2  )-  (SUMX* 2  )  ) 
a  COMPUTE  'A',  THE  Y -INTERCEPT 

A*YBAR-(.BxXBAR) 

F<THE) LEAST  SQUARES  ESTIMATED  REGRESSION  EQUATION  IS'.'.UTCNL 
'Y  =  «.(5»A),'  +  ' AS*B)'X.',VTCNL 

'DO  Y0b  WISH  TO  ENTER  SOME  X  VALUES  TO  GET  THE  PREDICTED  Y"S?  (Y/N) .  < 
WW*0 

*(WW='N<  )/Ll 
L2:' ENTER  X  VALUES.' 

XX*Q 

a  CALCULATE  PREDICTED  Y'S 

YY*A +B*XX 

'THE  PREDICTED  Y  VALUES  ARE:  ' .(vYY) .DTCNL 

'  WOULD  YOU  LIKE  TO  RUN  SOME  MORE  X  VALUES?  (.Y/N).  ' 

ww*a 

-■'WW='Y'  )/L2 

Ll: 'WOULD  YOU  LIKE  TO  TEST  HYPOTHESIS  ON  B,  THE  SLOPE  OF  THE  EQUATION?  07 
N).' 
ww*a 

*(WW='N>  )/L3 

'ENTER  THE  HYPOTHESIZED  SLOPE.' 

3B*  □ 

a  ,  .  COMPUTE  Ul'S 

U-Y-(BBxX) 

CHA* 1 >  1 

a  CALL  SPMANP  TO  COMPUTE  R  AND  ASSOCIATED  P-VALUES. 

D*X  SPMANP  7 
*(Dri:>0 )/L 11 
CHA* 1 <  1 
£11 :PV*2*DZ2l 
*(PVS1)/L 19 

PV*  1  _  _  .  . .  . 


£19 : 'SPEARMAN' ' S  R  EQUALS:  ' . (4»C C13 ) ,  1TCNL 
'THE  P-VALUE  FOR  HO:  B  =  '  .  ( *BB  ) .  1  VERSUS' 


Hi:  B  '.(.9CBA  [1  2]77£*BBV  IS:  <  .  (4»C[2]  )  ,Or£W£  ^  „ 

THE  TWO-SIDED  TEST  OF  HYPOTHESIS  IS:  '  ,('*9(2xDl2l  ))  ,DTCN 


'THE  P-VALUE  FOR  THE  TWO-SIDED  TEST  OF  HYPOTHESIS  IS:  ' , (u«(2*£>[2]  )  ), 

aL  IF  USING  THE  NEW  REGRESSION  EQUATION  BASED  ON  MEDIANS,  EXIT  HERE. 
L3:*(FF='N' )/L 18 
'PRESS  ENTER  WHEN  READY. ' 

ww*a 
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a  COMPUTE  CONFIDENCE  INTERVALS  ON  B 

£18 : 'WOULD  TOO  LIKE  A  CONFIDENCE  INTERVAL  FOR  THE- SLOPE?  (Y/N). ' 

WW*  0 

*0 W='AM  5/W1 
£10  :CC*INPUT  5 

a  CHANCE  ENTERED  VALUE  TO  ALPHA 

A£PBJl*(l00-CC)  +  200 

a  B0WB  TO  NORMAL  APPROX .  FOR  CONF .  IBZ\  OF  LARGER  SAMPLE  SIZES 
♦(W>13)/£5 

a  COMPUTING  CONFIDENCE  INTERVALS  BE  EXACT  P-VALUE 

P+RENDAL?  N 
CDF*P [2;: 

A  POSITION  OF  VALUE  IN  CDF  S  ALPHA 

INDEX*  <,*/ (CDFSALP3A ) ) 

*(INDEX>0 )/£6 
INDEX* 1 
+■£  6 

a  COMPUTING  CONFIDENCE  INTERVALS  USING  NORMAL  APPROX.  W/C.F. 

LS:DENOM*((Nx(.N-l)x(  ( 2*W 5+5 5 5 +19 5*0 . 5 
TALPHA-DENOMx ( | ( NORMPTH  ALPHA  )  ) 

-£3 

£6  :  rX£?ff  A+-P  C  3  ;  IJVBBX] 

TALPHA 

L3:CI*X  CONFLR  I 
NN*1 *CI 
SLOPES*\*CI 
RR* L ( iNN-T ALPHA 5+2) 

*^RR*0 5/£20 
RR*1 

1  £20  Cl+(  (NN+TALPBA)*2  ) ) 

SR 

->SRS(s SLOPES' )/£21 
SB+-oS£JPBS 

£21:  M  ’.(•CC),'  CONFIDENCE  INTERVAL  FOR  B.  THE  SLOPE  OP  ' 

'  TEE  ESTIMATED  REGRESSION  LINE,  IS:',QTCNt 

<  <  •  (5«5£0P£S;PP]  5 ,  '  S  B  S  '  ,  (S9SLOPESISR1  ),  '  ).',0T<;.V£ 

' PRESS  ENTER  WHEN  READY . 1 
WW*$ 

a  IP  5  OUTSIDE  THE  C.I.  CALCULATE  NEW  EQUATION  BASED  ON  MEDIANS 
-U32SLOPESZRR2  )a (3SSL0PESCSP.:  ))/N  1 
a  ORDER  X  AND  2 

X-+XTAX] 

A  ,  CBBCX  TO  SEE  IF  THE  SIZE  OF  SS  IS  EVEN  OR  ODD  FOR  FINDING  MEDIANS 
+  <(2|JW)  =  0)/S1 

a  ,  COMPUTE  MEDIAN  FOR  ODD  CASE 

B*SLOPESZ  (  (MV+1 5+2  5] 


DO  THE  SA _  _  _ 

S2  :♦(  (2lJV)  =  0  )/S3 
JB4B-+JL  ( (N+15+23] 

Xfl4H+-X[((»+l5  +  2)] 

*OUT 

53:XSAP-+(TC(P+2)]+rC((P+2)*2)3  )*2 
XBAH+-CXC  (W+2 )] +XC  (  (N+2  5  +  2  )]  5  +  2 
»  ,  _D,D  COMPUTE  NEW  INTERCEPT  < A' 

OUT:A*IBAR-(BxXBAR) 

'THE  LEAST  SQUARES  ESTIMATOR  OF  B  LIES  OUTSIDE  THE  CONFIDENCE  INTERVAL 
'DISCARD  THE  LEAST  SQUARES  EQUATION  AND  USE: ' ,QTCNL 
I  =  ',(•<0.'  +  'J-*Ul,'X  1  ,'OTCNL 

'  THIS  EQUATION  IS  BASED  ON  THE  MEDIANS  OF  THE  X  AND  I  DATA. AND  TH 

'MEDIAN  OF  THE  TWO-POINT  SLOPES  CALCULATED  FOR  THE  CONFIDENCE  INTERVAL  0 
N  fl.  ' .0 TCNL 

A  ALLOW  USER  TO  DO  SOME  ANALYSIS  ON  NEW  EQUATION 
’  DO  YOU  WISH  TO  ENTER  SOME  X  VALUES  TO  GET  PREDICTED  Y"S  ' 

'FROM  THE  NEW  EQUATION?  (.Y/N).  ' 

FF*  0 

♦( PF-'Y' )/£2 
FF*'Y' 

♦£1 

7 


V  SICN ;A:C:B :D :PVAL ;X:MO:N :CDF: ALPHA ;CI:Y :AA  ;BB;CC:DD:PV', PVI ; NNN ; EPOS : ORD 
D ; 27 lORDX : RALPH A :Z:Zl:QUA ;WW;PVM : pfo ; A : 6 

*  function  uses  the  ordinary  Sign  test  to  calculate  the  k 

a  swrisrrc.  p-value.  and  confidence  interval  as  a  test  for  medians. 

a  THE  LAST  OPTION  WILL  DISPLAY  A  TABLE  OF  CONFIDENCE  INTERVALS  OF  ORDERED 
a  STATISTICS  WITH  CONFIDENCE  COEFFICIENTS. 

a  57  FUNCTION  INCLUDE:  BINOM,  NORMCDF,  NORMPTH. 

*N<*  ’ 

NX -.MENU  SICNHELP 


103  N*:C*CBOICSS 


asm 

ICES  PAC, 


WICES  AND  ROUTE  FOR  PROPER  ACTIONS 
'ACEDMENU  SIGNQS.J. 
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♦(0=2,3 ,4,5 )/£8 
♦(0=7.8,3.10)/£9 

*(0= ll)/£20 

MPJVO  MAINQM 

u*°  INPOX  DATA  FOR  SINGLE  SAMPLE  CASE 

Li:AA*l 
X*INPUT  1 
NNN*oX 
♦  (0= 5 5/116 
MO+INPOT  3 
D*X-MO 
*£11 

p  PAIRED  SAMPLE  CASE 

L<3:AA*2 
R* INPUT  2 

i*itfo+i)+p 

y*(q+i)+P 

dd*x-i 

NNN+oDD 
*(C= 10 5/£l6 
M0*INPrJT  a 
0*«-X5-M0 

a  COMPRESS  D  TO  REMOVE  ZEROS 

L11:A+(D*0 )/D 

p  RECORD  LENGTH  OF  A  AND  ASSICN  TO  N 

N+pA 

p  KEEPING  TRACK  OF  POSITIVE  SIGNS 

ZPCS++/(A>0 5 
-OV>  2  5  5  /NORM 
PVAL*3INOM  .V 
KPOS>0  J/Pl 

pT/7*i 

pT?P7X*i -pimxcppos;; 

P2 :  PV*PVAL  L  (.KP0S+ 1 5  3 
*£6 

p  x?  .v  is  opp-irsp  th>iw  30  osz  normal  apprqx  w/  continoitt  correction 

NORM:  31-<  (XFOS-O  .  S 5- ( 0 . 5*AJ  ))♦  (C  .  S *  (tf*0 . 5  5 ) 

2*L  (KPOS+O  . 5 5- CO .3*^5  5+ (0  .  =*(.iV*0 .5)5 

pv+normcdp  z 

PVI+l - (NORMCDF  215 

p  IP  PAIRED  SAMPLE  TEST  GO  TO  L 17  POP  OOTPOT  STATEMENT 

£6:P73+2*(L/(Pl\P7X)) 

♦(P73S1)/M5 
P73*l  ,  %  , 

J?5 :P7M*( 3.1  )o  (P7.P7X.P73  ) 

' COMPUTATIONS  AkE  BASED  ON  A  SAMPLE  SIZE  OP:  ' .  («tf) ,  ’ • ' .UTCNL 
'J^BEfOTJ^  WVLBER  OF  POSITIVE  SIGNS  IS:  ' .  ( *KPOS )  ,  UTCNL 

*TBE? P-tALUE^FOR  HO:  M  =  '.(*M0).'  VERSUS  Hi:  M  ' , (»£OCXO[ (0-1 ) ; 1] ) , '  ' 
•MO ) , '  IS:  '  ,(5*P7MCC-l5l]5,Q2,OP£ 

♦£  18 

£17: 'I77F  P-VALUE  FOR  COMPARING  THE  MEDIAN  OF  THE  POPULATION  OF  ' 
'DIFFERENCES  TO  THE  HYPOTHESIZED  MEDIAN.  '.UTCNL 

'HO:  MU-1)  =  '.(•MO),'  VERSUS  HI:  M(X-I)  ' , (•XOCXOC (0-6 ) ; 1] ) , '  ',(»M0) 
,  XS:  '  (5»P7MC(C-6 );1] )TuToro 

£18: 'WOULD  IOU  LIKE  A  CONFIDENCE  INTERVAL  FOR  THE  MEDIAN?  CI/N) . ' ,UTCNL 


•  .  V  7 

£18: 'WOULD  IOU 
BB*  g 

*(BB='I' )/£ 16 
+QUANT 

£16 :CC*INPUT  5 

aZppwiqo-oo 


XPPOT  SXZP  OP  CONFIDENCE  INTERVAL 


ALPHA+\100-CC)*200 

+INNN>2S)/N0RM1 

a  COMPUTING  CONFIDENCE  INTERVALS  BI  EXACT  P-VALUE 

CDP*BINOM  NNN 

p  ,  INDEX  POSITION  OF  CDF  FOR  ALPHA  *  2 

B++ /  (CDPSALPBA') 

*(.B>0>  /  SKIP 
3*  1 
■+SKIP 

o  COMPUTING  CONFIDENCE  INTERVALS  BI  NORMAL  APPROX. 

NORM!  :  RALPH  A*  <  (0  .  *s*<NNN*0 .515*  ( NORMPTR  ALPHA)  WO  .  ‘xWlf'-O.  5 
a  ROUND  ’’.ALPHA  DOWN  TO  NEAREST  INTEGER  AND  INCREMENT  31  ONE 
3—\ EALPBA *  1 

a  IF  SINGLE  SAMPLE  CASE  GO  TO  27 

SKIP:*(.C=2 ,3 ,4 , 5  )/£7 

p  CALCULATE  AND  PRINT  OUT  CONF.  INT .  FOR  PAIRED  SAMPLE  CASE 

LS:ORDD*DDliDDl 
II*0ORDD 

CI+ORDDiB] ,0P00C(XX-(S-1))] 

1  A  ' .(•OO).'  CONFIDENCE  INTERVAL  FOR  THE  MEDIAN  OF  THE  ' 

'POPULATION  OF  DlPPERENCES  IS: '.UTCNL 

<  X  ',(*cxtl])f'  S  MEDIAN (X-X)  £  ' , (*0X C2]).'  )' .UTCNL 


+QUANT 

LT.ORDX+XltXl 

II+oQRDX 


UTCNL 


CALCULATE  AND  PRINT  OUT  CONF.  INT.  FOR  ONE  SAMPLE  CASE 


H~?R,f*'cc]  \°'RD  confidence  )  interval  for  the  median  of  TEE  POPULATION  IS:  ' 

f02^  (  '  (•CJClD),’  £  MEDIAN  &  <  AvCIZ 2]).'  )'.OrCW£  _ 

QUANT: 'MOULD  TOO  LIKE  cdNFIDENCE  INTERVALS  FOR  A  SPECIFIED  QUANTILE?  C I/N 
)  .  ' . Q TCNL 

ww*is 

♦AJ4 

120: 'ENTER  TEE  SIZE  OF  TEE  SAMPLE 
MW*-0 

El:' ENTER  DESIRED  QUANTILE :  FOP  EXAMPLE:  ENTER  20,  FOP  TEE  20 TE  QUANTILE. 

t 

♦  ( ( CCMSO  )  v  f  OCM >  1 0 0  )  ) /Fl 
mw  ccmjvc  <jm 

I  *****  TJJJS  TABLE  GIVES  CONFIDENCE  COEFFICIENTS  FOR  VARIOUS  INTERVALS  ' 
'WITS  ORDER  STATISTICS  AS  END  POINTS  FOR  TEE  ' , (vQUA ) ,'IU  QUANTILE.  < .QIC 
NL 

' PRESS  ENTER  VEEN  READS . ' 
bb*  a 

*N  4 

El -.'ERROR:  TEE  QUANTILE  VALUE  MUST  LIE  BETWEEN  0  AND  100;  TR1  AGAIN.' 

*B  1 
V 


7  S?MAW:X:I:4:!J:P:CPjI:3B:3:.P7 

a  THIS  FUNCTION  COMPOTES  TEE  SPEARMAN  R  STATISTIC  NR ICE  MEASURES 

a  TEE  DEGREE  OP  CORRESPONDENCE  BETWEEN  RANKINCS  OF  TWO  SAMPLES.  TEE  ?- 
a  VALUE  IS  GIVEN  FOR  TESTING  ONE  AND  TWO-SIDED  SYPOTEESIS  OF  ASSOCIATION 
a  SUBPROGRAMS  CALLED  BY  TEIS  FUNCTION  INCLUDE:  TIES,  IIESK ,  SPEARP , 

A  INPUT,  SPAPROX,  INTERP,  AND  TEE  VARIABLE  PMATSP. 

a  DISPLAY  MENU  AND  INPUT  DATA. 

N1 : B*MENU  SPMANOBJ 
*{b=1)/B1 
MENU  MAINQSi 

Bl:R*INPUT  2 
Q*l+R 

X*1*(Q+1)+R 

Y*(Q+1)*R 

A  CALL  SPMANP  TO  CALCULATE  TEE  STATISTIC  AND  ASSOCIATED  P -VALUES 
A*X  SPMANP  Y 
*UZ12>01/L1 
CEA*1 INDIRECT1 
*L2 

LI : CEA* 'DIRECT' 

L2:PV*2*Al21 

*gvsi)/L3 

L3:*SPEARMAN<  ' S  R  EQUALS:  ' , <4«*ri] ) ,QTCNL 
'TEE  P-VALUE  FOR  NO:  NO  ASSOCIATION  EXISTS  VERSUS ’  ,  ,  , 

'  El:  9CEAI ASSOCIATION  EXISTS  IS:  1  ,  (4»4  [21 )  .OTCJVC 

'TEE  P-VALUE  FOR  TEE  TWO-SIDED  TEST  OF  EYPOTBESIS  IS:  '  ,  (4*  (.2*A12  j T )  ,02 


SIS:  '  ,(4»4[2l).0rCJV£ 
IS:  »,t4*(2*At2jT),arc 


'PRESS  ENTER  WEEN  READY. ' 
BB*  Q 
*N  1 
V 


7  WISIC ;A:B:D:E: F:PV 2 ; Z 1 ; Z : DEN ; NUMZ :NUMZ 1 : PVAL ;X:MO;N; TPLUS :CDF :TALPEA : AL 
PBAiB:C}:Y;AA;B8:CC:NN:Db;PV;POS;TPOS;NM;PVI;TPOSUf/NN;C;PVM;PV3;B;Q;NU 
M:TC  :TC  1  -.TRAP  :DENT:DENT1 

a  TEIS  FUNCTION  USES  tEE  WILCOXON  SIGNED  RANK  TEST  TO  -CALCULATE  TEE  TPLUS 
a  STATISTIC.  P-VALUE,  AND  CONFIDENCE  INTERVAL  AS  A  TEST  FOR  MEDIANS. 
a  SUBPROGRAMS  CALLED  BY  IRIS  FUNCTION  INCLUDE:  TIES,  WILP ,  NORMCDF , 
a  NORMPTB ,  CONPW,  INPUT  AND  TEE  VARIABLE  PMATRIX . 

*<V  4 

.Vl  -.MENU  WILEELP 

a  "MENU  CHOICES  AND  ROUTS  FOR  PROPER  ACTIONS. 

N'i:C*CBOICEW  PAGEDMENU  WILQSi 
♦CC=2,3,4,5)/£8 
♦  (C=6 , 7 , 8 , 9 )/L9 
*{C=1  )/fll 
MENU  MAINQH 

9*  INPUT  DATA  FOR  SINGLE  SAMPLE  CASE 

L3:AA*1 
X*INPUT  1 
NNN* oX 
♦(C=5}/£16 
MO*INPUT  3 


PAIRED  SAMPLE  CASE 


20 
21 
22 
23 
l2“ 
2  5 
26 

27 

28 
29 

[30 

31 

32 

33 
3“ 

35 

36 

37 
.38 
139 
?40 
,41 

,'-•4 

>3 

>4 

>5 

46 

47 

48 
[ug 

50 

[51 

:?2. 

33 

54 

55 

56 

57 
53 

.59 

[60 

,61 

62 

,63 

,64 

65 

66 

67 

68 
69 
70„ 
7lj 
72*1 

73 

74 

75 

76 
,77 

78 

79 
,80 
,81 
,82 
,83 
.84 

85 

86 
,87 
88 
89 
.90 

91 

92 
.93 

IS 


D*X-MO 

♦Ill 


£9  :JUt*2 
R*INPUT  2 


!♦( 

DD+X-Y 
NNN*oDD 
♦  ( C=9  )/£16 
MO*INPUT  4 
D*(X-I)-MO 


Lll:A~(D*0)/D 

a 

N+oA 

*POS+(A> 0) 


COMPRESS  D  TO  REMOVE  ZEROS 
RECORD  LENCTB  OP  A  AND  ASSIGN  TO  N 
KEEP INC  TRACK  OP  POSITIVE  SIGNS 


a  TAKE  T3E  ABSOLUTE  VALUE  OP  A;  ASSIGN  TO  3  AND  ORDER  3 

B+IA 
B*B Z&Bl 

a  Reorder  positive  signs  to  coincide  wits  proper  positions  in  b 
pos*posu(  u): 

A  CALL  FUNCTION  TO  BREAK  TIES 

£♦1  TIES  B 

a  CALCULATE  TPOS  B1  ADDING  ACROSS  ALL  POSITIVE  VALUES  OP  E 

TPOS**/ (POSxE) 

TPOSl*TPOS 

a  GIVES  SIZE  OP  LEFT  TAIL  OP  PROBABILITY  DISTRIBUTION 

NM*(i((*/\N)*2))+l 

a  GO  TO  STATEMENTS  BASED  ON  LENCTB  OP  VECTOR  E 

*(N>9)/L3 

a  GENERATE  NULL  DISTRIBUTION  FOR  TPLUS 

F+WILP  N 

A  IP  TPOS  PALLS  IN  LEFT  BALP  OP  PROB  DIST  CALCULATE  RVALUE  AS  NORMAL 
L1:*(TP0SS  (NM-lj  )/TREC 

a  OTHERWISE  USE  TBE  NEGATIVE  T  STATISTIC 

TPOS*(*/\N ) -TPOS 

i  IP  TPOS  IS  FRACTIONAL  USE  BOTB  TBE  INTEGER  ABOVE  AND  BELOW 


TP.EG :♦<  (1 1  TPOS ) = 0  J /NON 
TPLUS*LlPOS 


AS  TPLUS 


Pl:PV*l-((PlTPLUSl*P£TPLUS*%)  )*2) 
PVI+(F£TPLUS+ll +PLTPLUS+21 )+2 

NON:*(TPOS>0  )/CO 


PV* 1 
♦P  2 

GO:PV*2-Pt(TPOS)\ . 

P 2 :PVI+F£lTPOS+l)l 
CHECK :  ♦  (  TPOS  IS  (NM- 1 )  )/£6 
PV2*PV 
PV*PVI 
PVI*PV 2 
♦£6 

a  COMPUTE  NORMAL  APPROX.  W /CONTINUITY  CORRECTION  FACTOR 
£3  :TRAP*(Nx  (N*l  5  3  +  4 
NUMZ*(TP0S*0 . 5  3 -TRAP 
NUMZ1*(TPOS-0 . 5 ) -TRAP 
DEN*((Nx(N*\)x((2xN)*l))*2'*)*0.S 
Z*NUMZ*DEN 
Z1+NUMZl*DEN 

A  COMPUTE  STUDENT  T  APPROXIMATION  WITH  CONTINUITY  CORRECTION  FACTOR 
NUM*I (TPOS-TRAP) , 

DENT*(X  (Nx  (DEN* 2  )  H (N- 1 ) )-  ( < (NUM- 0 . 5  )*2  )  +  (N- 1  )  )  )*0 . 5 
DENT1* ( ( (Nx (DEN* 2) )+ (N- 1 3 )-  (  ( (NUM*0 . 5 )*2 )+ (N-l 3  3  )*0 . 5 
TC* (NUM-0 . 5 ) *DENT 
TCl*(NUM*Q . 5 )*DENT1 
A  COMPUTE  AVERAGE  QP  TC  AND  ZC 
*\TPOSS,l  (+/\N  )*2  )) /SECOND 

PV*((\-(NORMCDF  il)l+(l-((N-l)  TDISTN  rci))>*2 
PVI*  (  (NORMCDP  Z)*((N- 1)  TDISTN  TC))*  2 
♦£  6 

SECOND :PV*((1- (NORMCDP  Z1))*((N-1)  TDISTN  TC1))*2 
PVI*(  (NORMCDP  Z)*(.\-((N-\)  TDISTN  TC)))*  2 
La :PV2*2x (i / (PV .PVI ) ) 

♦<  ?T/3£i  )/‘N5 
PV  3*1 

N 5 :PVM*( 3.1 )p (PVI.PV.PV3  ) 

' COMPUTATIONS  ARE  BASED  ON  A  SAMPLE  SIZE  OF  ' ,(wN ) , ' . ' .UTCNL 
'TBE  TOTAL  SUM  OP  POSITIVE  RANKS  IS:  ' .(*TPOh\) .uicNL 
"  ...  IP  PAIRED  SAMPLE  TEST  GO  TO  L 17  FOR  OUTPUT  STATEMENT 

♦<AA=2  )/£17 

iicVy^un^uih'.tiisv-'  k!sbs,,! » ■.  <•«««<:- !,»>.•  ,,<. 

♦£18 

Ki'fekc-iSWfireM*585  S1’  1  A’OOOICU-S.I l),'  '.(MO).*. 


86 


3 


•a 


m  ^ 


£18 : 'WOULD  ZOU  LIKE  A  CONFIDENCE  INTERVAL  FOR  IRE  RED I AN?  ( Z/N ).' 

ss+a 

*(BB='Z')/L 16 

*N  u  ° 

£16 : CC*INPUT  5 
ALPHA*  (1QQ-CO+200 

A  ROUTS  TO  NORMAL  APPROX .  FOR  CONF  INT  OF  LARCE  SAMPLE  SIZE 

♦  <AWJV>9)/£4 
CDF+WILP  NNN 

A  INDEX  POSITION  OF  CDF  FOR  ALPHA  *  2 

CDF*(CDF*0)/CDP 
TALPHA*( f / (CDFSALPBA ) ) 

*(TALPBA*Q j/JUMP 
TALPHA* 1 
*JUMP 

a  COMPUTING  CONFIDENCE  INTERVALS  BE  NORMAL  APPROX.  W/C.F. 

£4 : TRAP* ( NNN * ( NNN* 1 )  )  ♦  4 
DEN*( (NNN* (NNN+1 )x ( (2*NNN)*1 ) )  +  24 )*0 . 5 
TALPHA* (DEN* (NORMPTS  ALPBA ) ) *TRAP- 0 . 5 
TALPHA 

A  ROUND  TALPHA  DOWN  TO  INTEGER  VALUE  AND  INCREMENT  31  ONE 

TALPBA*ITALPBA*1 

a  IF  ONE  SAMPLE  CASE  GO  TO  £7 

IUMP:*(C=2 ,3,4,5 )/£7 

a  CALCULATE  AND  PRINT  OUT  CONF.  INT.  FOR  PAIRED  SAMPLE  CASE 

£5 :CI*TALPBA  CONFW  DD 

'  A  1 ,(9CC).<  CONFIDENCE  INTERVAL  FOR  TBE  MEDIAN  OF  THE  ' 
'POPULATION  OF  DIFFERENCES  ISi',DTCNL 

'  (  (sCICl]),'  i  MEDIAN (X- Z)  i  ' , («CJC2] ) , '  )',0 TCNL 

'  PRESS  ENTER  WHEN  P.EADZ .  ' 

3B*a 

*(V4 

a  CALCULATE  AND  PRINT  OUT  CONF.  INT.  FOR  ONE  SAMPLE  CASE 

L?:CI*T ALPHA  CONFV  X 

0 '  C0NFIDENCE  INTERVAL  FOR  TBE  MEDIAN  OP  THE  POPULATION  IS'.'-  . 

'  (  ' ,(«crci: ) . '  i  MEDIAN  Z  ' , (9CIZ21 ) , 1  )',0 TCNL 

' PRESS  ENTER  WHEN  &EADZ. 1 

BB*Q 

*N  4 


IjWj  »-  «.«  * 


«'i  •»!  4  M  t.l  M  •.(  i  '%  llj  «J  *-  ^  rtti 


APPENDIX  E 

RAIN  PROGRAM  LISTINGS  FOR  MAINFRAME  COMPUTER  WORKSPACE 


7  JUWM££:A:ilil;fl;BX;BrjCiCX:Cr;Z>;Z>D;:>X;I>ri0Xy;S;POS;/reCsXXiX7;iV;i>£WsAW;JVa 
H : P  :PV:PviL-Ja;SV:T:b:1r}Ai'iZiX;XiCBAi<)ik 
A  MS  FUNCTION  COMPUTES  THE  KENDALL  B  STATISTIC  WHICH  IS  A  MEASURE 
a  OF  ASSOCIATION  BETWEEN  SAMPLES.  P-VALUES  ARE  GIVEN  FOR  TESTING  ONE 
a  AND  TWO-SIDED  HYPOTHESIS  FOR  NO  ASSOCIATION  VERSUS  ASSOCIATION . 

A  SUBPROGRAMS  CALLED  31  THIS  FUNCTION  INCLUDE :  TIES,  TIESK,  KENDAL? . 
a  INTER P  INPUT  AND  NORMCDF . 

R-INPUT  2 


$*l+fa-H)+R 

Y+(Q+1)*R 

"a-^tcax: 

A  C 

3*XZ&X1 

"c+iUA 


ORDER  I  IN  INCREASING  ORDER  OF  X 
ORDER  X  IN  INCREASING  ORDER 
COMPUTE  CURRENT  RANKING  OF  I 


a  NOW  ORDER  I  RANKS  IN  INCREASING  ORDER 

D+AZ&A 1 

A  IF  TIES  EXIST  IN  EITHER  X  OR  I  RANKED  VECTOR  USE  MID -RANK  METHOD 
DD+ 1  TIES  D 
XX+i  mISS  3 

a  *  FIND  ORIGINAL  RANKING  OF  I  WITH  TIES  RESOLVED 
II+DDZCl 
N+oX 

a  COMPUTE  NUMBER  OF  DISTINGUISHABLE  PAIRS 
NN*(N*  (N-l))*2 
S+p  0 
AA+0 

a  POSITIVE  ONES  COME  FROM  A  RUNS  UP  CONDITION :  NEGATIVE  1  FROM  RUNS  DOWN 
A  ZERO  IS  SCORED  FOR  TIES.  MULTIPLY  THE  RESULTS  FOR  EVERY  ELEMENT  AND  SUM 
Ll:AA*AA+\  . 

3  <  [aa +XX  ] }  *  C 1 ) 

DXY*DX*DY 

P0S*(DXY>0) 

NEG*(DXY<0 )*<1) 

S+S.POS.NEC 

a  SUM  FINAL  VECTOR  TO  DETERMINE  S 

S**/S 

A  OBTAIN  THE  NUMBER  OF  TIES  IN  EACH  VECTOR  USING  THE  TIESK  FUNCTION 
U+TIESK  B 
V*TIESK  D 
SU++/(2\U) 

S7*+/(25K) 


S7*+/(2!7 
A  CALCULATi 
T*S*X(NN~. 


CULATE  THE 
UNN-SU)xXl 


AT*\t 

*(N>12  )  /NOR, 
A  CALL  . 

P+KENDALP  N 


'  8  STATISTIC  INCLUDING  THE  CORRECTION  FOR  TIES 
NN-SV))*Q.S 


N>l2)/NORM 

CALL  KENDALP  TO  CALCULATE  THE  RIGHT  TAIL  OF  THE  CDF  OF  B 


a  CALL  INTERP  TO  CALCULATE  P -VALUE  BY  INTERPOLATION 
PVAL*AT  INTERP  P 
+IPVAL*  1 )/L3 
PVAL*0 . 5 

-C  2 

a  CALCULATE  ?  VALUE  USING  NORMAL  APPROX. 

NORMiNUM*  1*AT)*<.  . 2*MV)*0.  i  j 
DEN*(2*(<.2*N)  +  S  )  )*0. 5 
Z*NUM*DEN 

PVAL*i- (NORMCDF  Z) 

A  IF  B  IS  POSITIVE  PRINT  OUT  DIRECT  ASSOCIATION . 
L3:*(T>0}/LS 
CHA*' INDIRECT' 

*LV 

LS-.  CHA*' DIRECT' 

LV :PV+2*PVAL 

L8 : ’KENDALL’ 'S  8  EQUALS  ',<4*r) 


'TBE  P- VALUE  FOR  BO:  SO  ASSOCIATES  EXISTS  VERSOS' 

|  Si:  1  .  (tea A) ,  'ASSOCIATION  EXISTS  IS:  '  .  (mPVAL) 

<TBE  P-VALUE  FOR  THE  TWO-SIDED  TEST  OP  EIPOTBESIS  IS: < , (4«P7) 


7  KRUSKAL :NUM:DENOM :A:C:B :D : K:AA :Bb:DD:E:P:NiOF:P; PVAL :R:SOFR : SR :TSOR ;CBA 
A  THIS  FUNCTION  COMPUTES  TBp  K'USKAL-U ALLIS  TEST  STATISTIC  B  WBICB  IS 
a  A  MEASURE  OF  TBE  EQUALITY  OF  K  INDEPENDENT  SAMPLES. 
a  SUBPROGRAMS  CALLED  BY  TBIS  FUNCTION  INCLUDE:  TIES.  TIESK,  INDEXPLS . 

A  PDISTN.  INTERP  AND  TBE  VARIABLES  PMATKW 20.  PMATKW 31,  PMATKW 3 3 , 

A  PMATKW3 4,  PMATKW'*!,  PMATKW*2 ,  AND  PMATKW*  § . 

QPP+S 

Si:' ENTER  TBE  NUMBER  OF  POPULATIONS  TO  BE  COMPARED  (.MUST  BE  GREATER  TBAN  T 
WO).  > 

*7(£<3)v((l!*>*0))/El 

♦((pS)>l)/Sl 

a  INITIALIZE  VECTORS  E  AND  F  AND  VARIABLE  C 

E*F*SOFR*o 0 

a THIS  LOOP  FACILITATES  ENTERING  TBE  SAMPLE  VECTORS  AND  STORINC  THEM 
CBA+' FIRST' 

£l:C+C+l 

' ENTER  YOUR  <  ,(*CB A),'  SAMPLE . ' 

D*-Q 

*<(oaO)*0)/NEXT 

D+laD 

a  CONCATENATE  SAMPLES  AS  TBEY  ARE  ENTERED  AND  STORE  TBEM  IN  VECTOR  E 
NEXT : 3*E , D 

a  RECORD  TBE  LENGTHS  OF  TBE  SAMPLES  AS  TBEY  ARE  ENTERED 

F*F,0D 
CBA*' NEXT' 

*(C<(K- 1 ) 5/£l 
CBA*-' LAST' 

-r(C<K)/L  1 

A  RECORD  SIZE  OF  ALL  SAMPLES  WBEN  COMBINED 

N-^/F 

a  ORDER  SAMPLE  SIZES  LARGEST  TO  SMALLEST 
OF+PlWPZ 

A  ORDER  COMBINED  SAMPLE  VECTOR  TO  BE  USED  BY  TIES  FUNCTION 

D*EZtE) 

a  CALL  INDEXPLS  TO  INCREMENT  INDEXES  WBEN  TIES  OCCUR  WITHIN  ONE  SAMPLE 
AA*P  INDEXPLS  E 

a  CALL  TIES  TO  BREAK  TIES  BY  MIDRANK  METHOD 

BB+ 1  TIES  D 
C*  0 

a  TBIS  LOOP  CALCULATES  TBE  B  STATISTIC 

£2:OC+l 

RSS^+/BBC(F[C]r+  AA  C™*^  BACB  SAHPLZ  IS  CALCULATED 
*  CALCULATeIu^O^PInkI  SQUARED  DIVIDED  BY  TBE  INDIVIDUAL  SAMPLE  SIZE 
SR*(SR*2)*FIC\ 


U+SCC^ 


SOFR*SOFR,SR 

*(C<K)/L2 

TSOR++/SOFR 


TORE  EACH  CALCULATION 


SUM  ACROSS  ALL  SAMPLES 


+  (TSOR*  (12* (Nx  (N*l))) )  -  (3x  (JV+1  )  ) 
RECALCULATE  B  WITS  CORRECTION  Ft 


A*TIESK  E 

NUM*(+/(A*3))-(*/A) 
DENOM+Nx((N*2)-l) 
H*B * ( 1 - (NUM*DENOM ) ) 
A  SYSTEM  OF  LOGICAL 


'OR  TIES 


STATEMENTS  ENSURE  PROPER  PROB .  IS  ACCESSED 


♦((a /(OF=  333  3))v(OBri]>3) 

(OF-  3  11  1;;v(.a /(OP-  3 


:/-*»“  1 

■‘kRxu.  )/?u3 

IF:* (OP [!]>“) IF  APPROX 
♦((a /{OF-  2  1  l))»(A/j 
♦((OS£i1=u)a?0W3]=1  )) 

*{(A/rof=  3  2  lJ)vfA/( 

*((A /(OF*  3  2  2  )M0Fr 
a  CALL  API 

P23 :P+PMATKW20l(N-* ) ;  j  j 
*PM 

P3 1 : P*PMATKW3 1 C (N- 5  ) ; ; ] 
P3 3 : P*PHATKW3 3 C (N- 5 ) | ; ] 


)/F APPROX 

2  i  i:)v(a/(ob= 

1))) /OUTPUT 
1))v(a /(0F=  3  3 


: ) v  c  a / (ob= 


3  3  1  ;  ) 


2))v(a  nop= 


3  )  )  )/P33 


VARIABLE  ACCESS  CDF 


P34  :P^PMATEW3'i  C  (S-6  ) : ;  3 
+PM 

P*1:P+PMATKW*1Z (S-5 ) ; ; 3 

P4 2 : P+PMATKW* 2[(S-5); ;] 

PU3:P*PMATKWU 3  C (S- 7 ) : : 3 
a  CXXS  INTERP  TO  tit 
PM: PVAL* g  INTERP  P 
* (PVAL- 1 ) /OUTPUT 


CULATE  P-VALUB  BI  INTERPOLATE OR 


W/ONE  LESS  D .P .IN  DENOH  APPROX. 


OUTPUT: PVAL*1 GREATER  TRAN  .25' 

L5: ;rss  s  srxrxsr/c  equals:  •.(mb) 

’TBS  P-7ALUE  POR  SO:  TBE  POPULATION  MEDIANS  ARB  EQUAL  VERSUS 

i  i 

'  Si:  AT  LEAST  TWO  POPULATION  MEDIANS  ARE  NOT  EQUAL  IS:  ' , 


. ( 9 PVAL ) 


El: 'ERROR:  10U  MUST  ENTER  A  SINGLE  INTEGER  VALUE  GREATER  TRAN  2:  TRI  AGAI 
N.  ' 

■  t 

♦SI 

7 


7  MANNWRIT : N ; M:PV 2 :A;B  :C  :G;MM-,NN  :RX;U : NM1 :P:NU: PVAL : NM ; NUMZ : NUMZ 1 DEN Cl 
;  2 1 :  ALPS  A : CD * :  INbEA  :tpX ;  Cl ;  UALpHA  :3B  :CC  ‘,Ul‘,U2iPV ;NN1  iNN2  :PVI ;  DIP F ;  ; 

;PVM :PV 3 :D:Q:R: DEN:DENC:TC: TCI :NUM 
A  THIS  FUNCTION  USES  TBE  SUM  OF  RANKS  PROCEDURE  TO  CALCULATE  TBE 
A  MANN-WHITNEI  U  STATISTIC  USED  IN  COMPUTING  TBE  P-VALUE  FOR  THE  TEST 
a  OF  LOCATION  AND  SCALE.  THE  C.I.  FOR  M(I)-M(X),  TBE  SHIFT  IN 
a  LOCATION.  IS  ALSO  COMPUTED.  SUBPROCRAMS  CALLED  INCLUDE:  TIES,  TIES 2 
A  INDEXPLS,  VARMW,  MANWP ,  INPUT,  CONPMU ,  NORMCDF ,  AND  NORMPTH . 

DIPF*0 

I  r 

'DO  10 U  WISH  TO  COMPARE  TBE  MEDIANS  OR  VARIANCES  OP  TBE  POPULATIONS?' 

i  i 

B2:'_  ENTER:  1  TO  COMPARE  MEDIANS i  2  TO  COMPARE  VARIANCES. ' 


1  TO  COMPARE  MEDIANS : 


2  TO  COMPARE  VARIANCES . ' 


X*1)aOA*2))/Z2 


S3:'  „  ENTER  TBE  DIPPERENCE  OP  MEDIANS  Of(X)  -  Mil))  OR  900  TO  QUIT).' 

,>/« 

♦(BISSs 900 )/0 

'THE  NULL  BIPOTBESIS  STATES  -  TBE  POPULATION  VARIANCES  ARE  EQUAL;  V(X)  = 
711) . 1 

i  i 

'  WBICB  ALTERNATIVE  DO  IOU  WISB  TO  TEST?' 

i  i 

J72: 'ENTER:  1  POR  SI:  7(X)  <  7(X);  2  POR  SI:  7(X)  >  7(1);  3  SOS  Bl :  7(X) 

*  7(1). ' 

)a(S*2 )a(0*3 ) )/S4 

Slf'  TSS  NULL  BIPOTBESIS  STATES  -  TBE  MEDIANS  OP  X  AND  I  ARE  EQUAL ;  M(X) 

7 

'  WBICB  ALTERNATIVE  DO  IOU  WISB  TO  TEST?' 

i  i 

S6:'XJVr^S:  1  POR  Si:  M(X)  <  M(I);  2  SOS  Si:  M(X)  >  M(I);  3  SOS  SI:  S(X) 

2T?s*i  )  a(S*2)a(D*3 ))/SU 

a  ^vrss  oxrx  vectors 

31 : ' ENTER  X  DATA  ( MORE  THAN  ONE  OBSERVATION  IS  REQUIRED).' 


'*!  ro3S)=o  )/si 

' ENTER  I  DATA. 


'Enter  i  data.' 

M*Q 

a  IS  CALCULATIONS  INVOLVE  VARIANCES  ADJUST  X  BI  TBE  DIPPERENCE  IN  MEANS 
N*N-DIPP 

A  CONCATENATE  X  AND  I  SAMPLE  VECTORS 

A*N,M 

A  DETERMINE  SIZE  OP  X  AND  I  VECTORS  AND  ASSIGN  TO  NN  AND  MM 

NN+pN 
MM*pM 

A  COMPUTE  SIZE  LIMIT  OF  LEFT  TAIL  OP  NULL  DISTRIBUTION 

NM*(NN*MM)*2 
NM1+INM 


ClU 


a  ORDER  A  AND  ASSIGN  TO  B 

cZil)NA,]tM)  INDEXPLS  A 

A  CALL  TIES  FUNCTION  TO  BREAK  TIES  USING  MIDRANK  METHOD 

G*1  TIES  B 

a  IF  FALSE  CALCULATE  TEST  FOR  VARIANCES 
*UjI  =  1)/35 

a  CALL  VARMH  TO  GENERATE  RANKS  REQUIRED  FOR  VARIANCE  TEST 
GG*VARM»~NN+MM) 

A  CALL  TIES  TO  RECORD  TIBS  IN  THE  DATA  AND  BREAK  TIES  IN  GG 
G*GG  TIES  B 

a  CALCULATE  SUM  OF  X  RANKS 

SS :RX*+/ (GC (C[l ; ( \NN )? ) 

a  CONVERT  TO  MANNWHIT  U  STATISTIC 

U*RX- ( (NNx (NN+1 ))*2) 

01*0 

a  IF  SIZE  OF  X  TIMES  SIZE  OF  1  >  80;  GO  TO  NORMAL  APPROX 

*( (WM*2 )>80 )/L2 
NNl*l/NN,MM 
NN2*\  /NN.MM 

aMANWP  FUNCTION  CALCULATES  LEFT  TAIL  CUMULATIVE  PR03S.  OF  U  STATISTIC 

a  LOGICAL  STATEMENT  ENSURES  ONL1  LEFT  SIDE  OF  NULL  DIST  IS  USED 
*(0ZNM1 )/L 3 

a  CONVERT  U  STAT  WREN  GREATER  TRAN  LEFT  TAIL  VALUES 

0*(NN*MM)-U 

A  IF  U 1  IS  A  FRACTIONAL,  INTERPOLATE  P  VALUE 


A  IF  U 1  IS  A  FI i 
L3:*dl\U)=0)/NON 

*\u2>0)/Pl 


82]  *(U2>0)/P\ 

831  PV*l-((Pr.U2+U  )*2) 

Is]  Pl:ly*l-((PC02]+i>Cff2+l]  )  +  2) 

86  ?3:PVI*(PlU2+ll*PlU2+2l)+2 

87  *CBECK 

88 J  NON:*(U>0)/GO 

89  PV* 1 

90  *P2 

91  CO:PV*l-prui 

92  P2:f>VI-*-P[;  {0+1)1 

93  CHECK  :*(U1SNM1)/L'+ 

94  PV2*?V 

95  PV*PVI 

96  PVI*PV 2 

97  +£u 

98  A  COMPUTE  THE  NORMAL  APPROXIMATION  W /CORRECTION  FACTOR 

99  L2:NUMZ*(.U+0.5)-NM 

100]  NUMZl*(U-0 . 5 )-NM 

101]  DEN*UMMxNNx  (HM+NN+1  )  )»12  )*0  •  5 


Z1*NUMZ1+DEN 
NUM*\(U-NM ) 

DENC*( ( ( (NN+MM-1 )> 


x (DEN* 2 ) ] ♦ (NN+MM-2 )  j  -  (  £ [NOM- 0 . 5 )*2 ) ♦ (NN+MM-2 ) ] ]*0 . 5 


DENC !♦(< ( (NN+MM-1 )x (DEN* 2 ) }♦ (NN+MM-2 ) 
TC*(NOM-0.S)+DENC 
TCl*(NUM+0 . 5 )+DENCl 
+IUSNM) / SECOND 

PVI+l (NORMCDF  Z)+((NN+MM-2)  TDISTN  TC))+2 
PV* ( < 1 - (NORMCDF  Z1 ) )+(l- ( (NN+MM-2 )  TDISTN  TC l)))*2 
♦£u 

SECOND : PVI+ ( (NORMCDF  Z)+(l-((NN+MM-2)  TDISTN  TC )  )  )♦ 
PV+((1- (NORMCDF  Zl}}+( (NN+MM-2)  TDISTN  TCl))+2 
M:PV2*2x ( L / (PV ,PVI ) ) 


(NUM+Q . 5 )*2 )♦ (NN+MM 


2  J  lift 


115  m:PV3+2x(\/  (PV  ,PVI )) 

116,  ♦(P73Sl)/ff5 

117  PV3+1  J 

118  NS:PVM*(3,l)p(PVI,PV.PV3 

119  • THE  SUM  OF  THE  X  RANKS 


TDISTN  TC)))*2 
N  TCl))+2 


]ll9]  < THE  SUM  OF  THE  X  RANKS  IS:  ' , (*RX), ' .  THE  0  STATISTIC  EQUALS:  ' . (vUl) 

120  '  ' 

:i2i:  A  LOGICAL  STATEMENT  FOR  VARIANCE  OUTPUT 

122  *  (  AA  -2) /VAR. 

]l23J  '  THE  P-VALUE  FOR  BO:  M(X)  =  M(7)  VERSUS  El:  M(X )  ' , (vLOGICtD ; 1] ) , '  M(I) 

IS:  ' ,(4*P7H[0;i]l 

ill, 

]l27]  7*TBEV?*kll,UEP^OR'HC* ’7«!  -  7(1)  VERSUS  31:  V(X)  '  ,  (iLOCICZD:  i:  ) . '  7(1) 
IS:  '  .  (‘**PVMZD:i:  J 
;i29]  '  * 

l129J  *0 

tl30 j  £8: 'WOULD  IOU  LIKE  A  CONFIDENCE  INTERVAL  FOR  THE  SHIFT  IN  LOCATION (M(7 )  - 
M(X))7  (Z/N).' 

131  BB+ Q 

132  *(BB-'N')/0 

133  L10:CC+INPUT  5 

134,  ALPBA*(100-CC)+200 

135  A  ROUTE  TO  NORMAL  APPROX.  FOR  CONF.  TNT.  OF  LARGER  SAMPLE  SIZES 
136,  *((NM*2)>iO)/LS 


CDF+P 

INDEX* ( 
*  (INDEX 
INDEX* 1 


COMPUTING  CONFIDENCE  INTERVALS  BT  EXACT  P-VALUE 

INDEX  POSITION  OF  VALUE  IN  CDF  S  ALPHA 
(  +  / (CDF&ALPHA)  ) 

X>0}/£6 


a  1  COMPUTING  CONFIDENCE  INTERVALS  USING  NORMAL  APPROX.  W/C.F. 

LS : U ALPS A* (DENOMZ* (NORMPTB  ALPBAl  )-NM-Q . 5  _  _ 

A  ROUND  UALPBA  DOWN  AND  INCREMENT  BT  ONE 

INDEX+lOALPBA+l 
L6-.IPX-NN  .INDEX 

<*IAI*X.(jicc¥!!'  PERCENT  CONFIDENCE  INTERVAL  FOR  TEE  SB  I  FT  IN  ' 

•  LOCATION  BETWEEN  POPULATIONS  X  AND  I  IS:' 

;  (  <  ' , (»crci] > . '  S  M(I)  -  Mix )  S  '.(«CIC23).'  )• 

■*■0 

El ‘.'ERROR:  TEE  SIZE  OF  IOUR  SAMPLE  IS  LESS  TRAN  TWO ;  TRI  AGAIN.' 

■  i 

♦SI 

E 3: 'ERROR:  IOU  RAVE  ENTERED  MORE  TRAN  ONE  VALUE.  TRI  AGAIN.  > 

i  i 

♦S3 

E2 : ' ERROR :  IOU  RAVE  NOT  ENTERED  A  VALUE  OF  1  OR  2;  TRI  AGAIN.' 

i  i 

-32 

EM :< ERROR:  IOU  RAVE  NOT  ENTERED  A  VALUE  OF  1,  2,  OR  3:  TRI  AGAIN.' 

!  I 

*(AA=2)/N2 

♦86 


7  NPSLR-.N :  SUMX  '.SUM!  ;XBAR :  13  AR :  SUMX2  :  SUMXI  :3-,A  :WW :  XX :  3B :  U :  D :  ALPS  A  \P:CC ;  CDF 
■.TALPiiA ;  NN  -.d:  SLOPES :  3R  :SR  -’ll  iDENOM-.. INDEX :  FF :  X  :I:Q  ;R:CBA:?V  „„„ 

n  PROGRAM  cbNDbCTS  NONPAkAMkTRi C  LINEAR  REGRESSION.  TEE  LEAST  SQUARES 
a  ESTIMATED  REGRESSION  LINE  IS  COMPUTED  WITH  RIPOTEESIS  TESTING  AND 
A  CONFIDENCE  INTERVAL  AVAILABLE  FOR  TEE  SLOPE  3.  IF  B  DOES  NOT  LIE  IN 
a  TBE  C.I.  AN  ALTERNATE  REGRESSION  LINE  IS  PROPOSED.  SUBPROGRAMS  CALLED 
a  ARE :  SPMANP,  KEN D ALP .  NORMPTB,  INPUT,  AND  CONFLR. 

QPP-5 

a  INPUT  DATA 

R-INPUT  2 
Q—l+R 

hsliM'"  .  . . 


SUMX-+/X 

SUMI-+/I 

*XBAR-SUMX*N 
IBAR—SUMI *N 

"SUMX2*+/(X*2) 


ASSIGN  TEE  SIZE  OF  X  {AND  I)  TO  N 
COMPUTE  THE  SUM  OF  X'S  AND  T'S 

COMPOTE  TBE  MEAN  OF  X  AND  I 


COMPUTE  TBE  SUM  OF  TBE  X'S  SQUARED 
COMPUTE  TBE  SUM  OF  X  TIMES  I 

,  TBE  SLOPE  OF  TBE  ESTIMATED  LEAST  SQUARES  REGRESSION  LINE 


aS  COMPUTE1 'S',  TBE  SLOPE  OF  TBE  ESTIMATED  LEAST  SQUARES  REGRESSION  LINE 
B+i  iNxSUMXD-iSUMXxSUMI}  )*((N*SUMX2  )  -  (SUMX* 2  )  ) 

A  COMPUTE  'A1 ,  TBE  1 -INTERCEPT 

A+IBAR- C BxXBAR) 

FF ♦ 1 N ' 

|  TBE  LEAST  SQUARES  ESTIMATED  REGRESSION  EQUATION  IS:  ' 

J  I  =  ' ,{9A),<  +  ' ,{vB),'X.' 

'DO  TOO  WISH  TO  ENTER  SOME  X  VALUES  TO  GET  TBE  PREDICTED  I"S7  il/N).' 
WW+Q 

-CWW='N' )/L\ 

L2:' ENTER  X  VALUES.' 

a*  CALCULATE  PREDICTED  I'S 

I7+A+B*XX 

'TBE  PREDICTED  I  VALUES  ARE:  ' , (All) 

i  i 

'  WOULD  IOU  LIXE  TO  RUN  SOME  MORE  X  VALUES?  il/N). ' 

WW*  Q 

)/L2 

LU 'WOULD  IOU  LIKE  TO  TEST  3TP0TBESIS  ON  3.  TBE  SLOPE  OF  TBE  EQUATION?  (1/ 
N).' 

ww*a 

*(WW= ' N ' )/L 3 

' ENTER  TBE  BIPOTBESIZED  SLOPE.' 

B3-Q 

A  COMPUTE  U£'S 

U*I-(BB*X) 

CBA* 1 >  ' 

a  CALL  SPMANP  TO  COMPUTE  RBO  AND  ASSOCIATED  P-VALUES. 

D-X  SPMANP  U 
♦  U>C1]>0)/I11 
CBA*'<  < 

Lll:PV+2*Dl2l 


♦<i>VSl)/£19 

PV*1 

£19:'  SPZARMAN"S  R  EQUALS:  ',(4*CC1]> 

1  TBS  P-VALUE  FOR  30;  B  =  *,<«3S)#'  VZRSOS  Bit  B  ' , (fCBA  Cl  2  IS 

t‘ i{*9DC21) 

'tTBZ  P-VALUE  FOR  TBZ  TWO  SIDZD  TEST  OP  BIPOTBZSIS  IS:  ',(4 vPV) 

A  IP  USING  TBZ  NEW  RECRZSSION  EQUATION  BASED  ON  MEDIANS,  EXIT  BZRZ . 

£3 :*(FF- ' I' )/0 

A  COMPOTE  CONFIDENCE  INTERVALS  ON  B 

' WOULD  100  LIKE  A  CONFIDENCE  INTERVAL  FOR  TBE  SLOPE?  (1/N).< 

WW*  a 

*(WW='N' )/0 
£10 : CC*INPUT  5 

A  CHANGE  ENTERED  VALUE  TO  ALPBA 

ALPBA+ (10  0-CC )  +  20  0 

A  ROUTE  TO  NORMAL  APPROX.  FOR  CONF .  INT.  OF  LARGER  SAMPLE  SIZES 
*(.N>12WLS 

A  COMPUTING  CONFIDENCE  INTERVALS  31  EXACT  P-VALUE 

P*KENDALP  N 
CDF*Pl 2:2 

A  INDEX  POSITION  OF  VALUE  IN  CDF  &  ALPBA 

INDEX* ( +/ <  CDFSALPBA) ) 

*(INDEX>0)/L6 
INDEX* 1 
*L6 

a  COMPUTING  CONFIDENCE  INTERVALS  USING  NORMAL  APPROX.  W/C.F. 

£3  :DENOM*(  CiVx  (N-\  )x  ( ( 2*iV  )  +  5  )  )*18  ) *0 . 5 
TALPBA*DENOMx  (  |  (NORMPTB  ALPBA  ) ) 

*L5 

Do :  TALP3A*?Z3  -.INDEXl 
£9 : CI*X  CQNFLk  2 
NN*1*CI 
SLOPES*! *C1 
RR* L ( ( NN - T ALPBA )*2  ) 

*(&ft*0)/£20 
RR*  1 

£20 : SR*r (!*( (NN*TALPBA )  +  2 ) ) 

~<SRZ(s SLOPES ))/L 21 
SR*pSLOPES 

£21:'  A  ' , (*CC)  '  PERCENT  CONFIDENCE  INTERVAL  FOR  B,  TBE  SLOPE  OF  ' 

'  TBE  ESTIMATED  REGRESSION  LINE,  IS:' 

I  I 

;  (  <  >, (•SLOPES CRR3).<  <  B  <  ' , (9SLOPESISR]) ,' 

a  IF  B  OUTSIDE  TBEC.I.  CALCULATE  NEW  EQUATION  BASED  ON  MEDIANS 
A+<  (BiSLOPESLRRlU  (BSSLOPEs\SP^p  ^ 

A  CHECK  TO  SEE  IF  TBE  SIZE  OF  SS  IS  EVEN  OR  ODD  FOR  FINDINC  MEDIANS 
*((2\NN)=U)/S1 

A  ,,,  COMPUTE  MEDIAN  FOR  ODD  CASE 

B*SLOPESl((NN+l)*2)l 
*S  2 

A  ,  COMPUTE  MEDIAN  FOR  EVEN  CASE 

SI :B*(SLOPESl(NN*2)l*SLOPESl((NN*2 )*2 )] )  +  2 
a  DO  TBE  SAME  FOR  TBE  X  AND  I  VECTORS 

S2  :♦(  (2lJV)  =  0  5/S3 

XBXB-*xt((iV+i)  +  2)] 

*OUT 

S3:IBAR*(ll(N*2)l*I(((N*2)*2)'i  )*2 
XBAR*(XtiN*2  JJ+ict  IN* 2  )  +  2  )]  )*2 
A  ,  COMPUTE  NEW  INTERCEPT  ' A • 

0UT:A*1BAR- (B*XBAR ) 

'TBE  LEAST  SQUARES  ESTIMATOR  OF  B  LIES  OUTSIDE  TBE  CONFIDENCE  INTERVAL.' 
J DISCARD  TBE  LEAST  SQUARES  EQUATION  AND  USE : 1 

;  t  i  -  ',(»*),'  +  ',< *s),'x  > 


<  '  ,(mSLOPES(RR1),'  <  B  <  '  .(vSLOPESlSR]),' 
.  CALCULATE  NEW  EQUATION  BASED  ON  MEDIANS 


1  =  ',(»*),'  +  ' , (*B  ) ,  '  X  ' 


'  TBIS  EQUATION  IS  BASED  ON  TBE  MEDIANS  OF  TBE  X  AND  I  DATA  AND  ' 

' THE  MEDIAN  OF  TBE  TWO-POINT  SLOPES  CALCULATED  FOR  THE  CONFIDENCE  INTERV 
AL  ON  3 . ’ 

f 

a  ALLOW  USER  TO  DO  SOME  ANAL2SIS  ON  SEW  EQUATION 
'DO  10U  WISH  TO  ENTER  SOME  X  VALUES  TO  GET  PREDICTED  1"S  FROM  TBE  NEW  E 
QUATION?  ( 1/N ).' 

*Cf}='1')/L2 
FF* ' 1 ' 

*L  1 


V  SIGN ;A:C:B:D xPVAL :X :M0 :N : CD? ; ALPHA: Cl : 7 : AA i SB ; CC ; DD ; PV ; PVI ; MW ; £POS ; ORD 
0  ■  ”  iviicTSIoh  bssk  Zu'C'o'ki^KR^il^TEST  TO  CALCULATE  TBE  K 

a  ORDERED  STATISTICS  WITH  CONFIDENCE  COEFFICIENTS . 

A  SUBPROGRAMS  CALLED  BY  TBIS  FUNCTION  INCLUDE :  BINOM,  NORMCDF ,  NORMPTB , 

A  INPUT, AND  QUANC. 

>  DID  700  ENTER  TBIS  PROGRAM  FOR  THE  SOLE  PURPOSE  OF  GENERATING  CON 

FIDENCE  INTERVALS  FOR  A  SPECIFIED  SAMPLE  SIZE  AND  QUANTILE?  ( Y/N).‘ 

WW*  E 

■*(UU~  »T*  )/B 4 

'  BC££  HYPOTHESIS  STATES  -  THE  POPULATION  MEDIAN  (M)  IS  EQUAL  T 

0  THE  HYPOTHESIZED  MEDIAN  (MO);  50:  M  =  MO.  « 

'  MATCH  ALTERNATIVE  DO  YOU  WISH  TO  TEST? 1 

S3:'  ENTER’.  1  F0H  31:  M  <  MO;  2  FOB  HI:  M  >  MO;  3  FOB  Hi:  M  *  MO.' 

c*g 

*(7c*l )a (£s2 )a (C=3 ) )/F3 

B  2  ENTER:  1  FOR  SINGLE-SAMPLE  PROBLEM’,  2  FOR  PAIRED-SAMPLE  PROBLEM.' 
AA+(3 

-C(AA=1;a(AA*2))/B2 
♦(AA=2 )/£9 

a  INPUT  DATA  FOR  SINGLE  SAMPLE  CASE 

X*INPUT  1 
NNN*oX 
MO*INPUT  3 
D*X-MO 


£9:B* INPUT  2 
3*1  aB 

2*1  +  (<3+1  )fR 

Y+lQ+D+R 

DD*X-7 

NNN*oDD 

MO*INPUT  u 

D*(X-Y)-MO 

£11:A*(£>=0 )/D 

A 

N*f>A 

*KPOS*+/(.A> 0  ) 
♦  (.N23  0  ) /NORM 
PVAL*BINQM  N 
*^KPOS>0)/P1 


PI :P7X*1 -PVALZKPOS] 
P2  ’.  PV*PVAL  t  (  EPOS + 1 )  ] 


PAIRED  SAMPLE  CASE 


COMPRESS  D  TO  REMOVE  ZEROS 
RECORD  LENGTH  OP  A  AND  ASSICN  TO  N 
KEEPINC  TRACE  OF  POSITIVE  SIGNS 


A  IF  N  IS  GREATER  THAN  30  USE  NORMAL  APPROX  V!  CONTINUITY  CORRECTION 
NORM’.  Z1*X  (.EPOS-0 .5)-(0.5xW))+(0.5*(M*0.5)) 
Z*((BPOS+0.5)-(0.5x/m+(0.5x<,V*0.5)) 


PV*NORMCDF  Z 
PVI*\- (NORMCDP  Zl) 


IF  PAIRED  SAMPLE  TEST  GO  TO  L 17  FOR  OUTPUT  STATEMENT 


L 6  :P73*2x  (L  /  (.PV , PVI )  ) 

*(Pr3Sl)/M5 

PV3*1 

NS : PVM*(3 ,l)p (PV.PVI .PV3 ) 

'COMPUTATIONS  ARE  BASED  ON  A  SAMPLE  SIZE  OF:  ' , («W) 

I  I 

1  THE  TOTAL  NUMBER  OF  POSITIVE  SIGNS  IS:  ',(.<•  EPOS) 

i  i 

*(AA=2)/£17 

'TEE  P-VALUE  FOR  HO  :  M  =  '  ,  ( *MO ) , '  VERSUS  HI:  M  '  ,  (.mLOGIClC:  1]  ) , 
f ' yIS: 1  ,(4*P7M[C;l] ) 

♦£18 

LIT. 'THE  P-VALUE  FOR  COMPARING  THE  MEDIAN  OF  TBE  POPULATION  OF  ' 

' DIFFERENCES  TO  THE  HYPOTHESIZED  MEDIAN , ' 


, ( *MO ) 


•30:  M( X-Y)  ■ 
:  .  .^tPVMCC: ; 


:,.*MO).'  VERSUS  31:  M(X-I)  '  ,  (9L0GICZC:  i: 


:  , >.  *MO  ' , 


L 18: 'WOULD  YOU  LIEE  A  CONFIDENCE  INTERVAL  FOR  TBE  MEDIAN?  (I/N). ' 
BB*g 

*(BB='7' )/L 16 
*QUANT 

a  INPUT  SIZE  OF  CONFIDENCE  INTERVAL 

£16  ’.CC+INPUT  5 
ALPHA* ( 100 -CC )♦  200 
*(NNN220)/NORM1 

a  COMPUTING  CONFIDENCE  INTERVALS  BY  EXACT  P-VALUE 

CDF*BINOM  NNN 

a  INDEX  POSITION  OF  CDF  FOR  ALPHA  ♦  2 

B*+/(CDFSALPBA) 

*(B>0)/SKIP 


B*  1 
+SKIP 

a  COMPUTING  CONFIDENCE  INTERVALS  BT  NORMAL  APPROX. 

NORM! :  RALPH  A*  ( ( 0 . 5  *  ( NNN*  0 . 5  ))  * ( NOP.MPTB  ALPHA )  )  +  ( 0 . 5  *  NNN  )  -0 . 5 
A  ROUND  RALPH A  DOWN  10  NEAREST  INTEGER  AND  INCREMENT  31  ONE 
B+lKALPBA+1 

A  IF  SINGLE  SAMPLE  CASE  GO  TO  LI 

SKIP:*(AA=l )/L 7 

A  CALCULATE  AND  PRINT  OUT  CONF .  INT .  FOR  PAIRED  SAMPLE  CASE 

LS-.OP.DD+UDl&DDh 
IT+oORDD 

CI*0RDDZ3\ ,0RDD;CIT-(3-1 ) )2  „„  „„„  , 

'4  ' .(*CC).'  PERCENT  CONFlbENCE  INTERVAL  FOR  TEE  MEDIAN  OF  THE  ' 

| POPULATION  OF  DIFFERENCES  IS 

'  (  ' . (9CIC13 ). '  £  MEDIAN(X-I )  £  ' , (9CIC21 ) , '  )' 

t  i 

*00 ANT 

a  CALCULATE  AND  PRINT  OUT  CONF.  INT.  FOR  ONE  SAMPLE  CASE 

L7 :0RDX*XCAX~ 

II*sORDX 

CI*0RDXZ5 }  ,0RDX((77-(3- 1 

M  '  ,  [aCC,  ,  1  PERCENT  CONFIDENCE  INTERVAL  FOR  THE  MEDIAN  OF  IRE  POPULATED 
N  IS:  • 

J  (  (  ' ,(*CICU ). '  £  MEDIAN  £  ',(nCIC23).'  )' 

QUANT: 'WOULD  70U  LIKE  CONFIDENCE  INTERVALS  FOR  A  SPECIFIED  QUANTILE?  (I/N 


Ql/ANT;' WOULD  TOO  LIKE  CONFIDENCE  INTERVALS  FOR  A  SPECIFIED  QUANTILE?  (7/N 
WV*  3 

*CWW='I'  )/Sl 

3U:< ENTER  DESIRED  SAMPLE  SIZE  (SINCLE  INTEGER  VALUE).' 

NNN*Q 

*i(pNNN)>l)/Eu 
*(  lljiV/ViV)*0  )/E^ 

31 : 'ENTER  DESIRED  QUANTILE;  FOR  EXAMPLE:  ENTER  20,  FOR  THE  20 13  QUANTILE . 


*°(\qUA£0  IvC  Off>I>100 ) )/El 
NNN  QUANC  QUA 


'  *****  this  TABLE  GIVES  CONFIDENCE  COEFFICIENTS  FOR  VARIOUS  INTERVALS  W 
ITH  ' 

J ORDER  STATISTICS  AS  THE  END  POINTS  FOR  THE  ' ,(vQUA) ,'TH  QUANTILE.' 

El? | ERROR :  THE  QUANTILE  VALUE  MUST  LIE  BETWEEN  0  AND  100;  TRI  AGAIN.' 

♦El 

E2:' ERROR:  IOU  HAVE  NOT  ENTERED  A  VALUE  OF  1  OR  2;  TRI  ACAIN . ' 

t  i 

*B2 

E3:’f  ERROR:  IOU  HAVE  NOT  ENTERED  A  VALUE  OF  1 ,  2 ,  OR  3 ;  TRI  AGAIN.  ' 

*B3 

EM :> ERROR:  IOU  HAVE  NOT  ENTERED  A  SINGLE,  INTECER  VALUE;  TRI  AGAIN.' 


7  SPEARMAN ;X:7;A:Q:R; CBA : PV 

A  THIS  Pl/N&Tloft  COMPOTES  THE  SPEARMAN  R  STATISTIC 

A  THE  DEGREE  OF  CORRESPONDENCE  BETWEEN  RANKINGS  OF  TVC 
A  VALUE  IS  GIVEN  FOR  TESTINC  ONE  AND  TWO-SIDED  B7PCTHE, 
A  SUBPROGRAMS  CALLED  37  THIS  FUNCTION  INCLUDE:  TIES.  T. 
A  INPUT,  SPAPROX,  INTERP,  AND  THE  VARIABLE  PMATSi  . 
R*INPUT  2 

I*(Q+1)*R 

a  CALL  SPMANP  TO  CALCULATE  THE  STATISTIC  AN V  ASS'~IATr 

-pmanp  : 


WHICH 
SAMFLl 
SIS  IF 
DESK.  S 


Ll:CHA*' DIRECT' 

12:PV*7*A17) 

*\PV£\)/ts 
PV*  1 

£3: 'SPEARMAN' 'S  R  EQUALS  ',(«•*: i: 3 

i  i 

'THE  P-VALUE  FOR  BO:  NO  ASSOCIATION  EXIST:  vrr_ 
;  El:  1 , ( *CBA ) , '  ASSCCIATI  a  F\ 

'THE  P-VALUE  FOR  THE  TWO-SIDED  TEST  CF  THF . 
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UNCLASSIFIED 


AN  APL  WORKSPACE  FOR  CONDUCTING  NONPARAAETRIC 
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F/G  12/3 


2/2 


LMiiUiiira 


f  fn  iHoH&KPBtfii  OWVT.'Br”"  mlu  0,1 "  T 

VtICt  ALTERNATIVE  DC  TOO  HZ  St  TO  TEST?  ’ 

t  1 

*|io  IWM:  1  ro*  B\i  R  *  MQi  2  rot  tu  H  *  N0|  9  /OX  Ml  s  t  •  mi.' 

n{  ,£r^ixtC*,iAISxsiij»Sx-siji/£X  problem i  2  tot  paired-sample  problem.' 


3iff5U2i""a),/w 


*rx/ar  : 


fO*r.ra 
3*a-mc 
•£  11 


2j:X»XX/0T  2 
0*1  ti 

/♦i*7q*i  )*X 
* 

NNN»oic 
*0*1*101  « 

3-  vi 

Ill:l*(0«O)/D 


:x/or  oxrx  /ox  sited  sample  case 


t AIMED  SAMPLE  CASt 


P0S*(A>01 


COMPRESS  D  TO  BEHOVE  ixtos 
tXCOMD  LMMGTi  Ot  A  AMD  ASSIGM  TO  X 
MEETING  TRACE  Ot  POSITIVE  SIGNS 
TAKE  TEE  ABSOLUTE  7 ALOE  OP  X:  ASSIGN  TO  B  AND  ORDER  3 


M  TIM 


(N»lt )/L3 


POSITIVE  SIGNS  TO  COINCIDE  NITt  PROPER  POSITIONS  IN  B 
CALL  PONCTION  TO  BREAM  TIES 
L^X  TPOS  ET  ADDING  ACROSS  ALL  POSITIVE  VALUES  OT  E 

mSIZE  OP  LETT  TAIL  OP  PRORARILITI  DISTRIBUTION 
'go' TQl STATEMENTS  EASED  ON  LMNGTE  OP  VECTOR  E 
GENERATE  NULL  DISTRIBUTION  POR  TPLUS 
1L£5  IN  LETT  SALT  OP  PROB  DIST  CALCULATE  PVALUS  AS  NORMAL 
nSMISE  USE  TEE  NEGATIVE  T  STATISTIC 

o/^J^XgXjJ£  USE  BOTE  TEE  INTEGER  ABOVE  AND  BELOV  AS  TPLUS 


*a£» *7rzos»o )/co 


ziic*i*(rposiscNM~i) : 


n*Py 

1*1*1 

rx-PVt 


}  ^illUOTT 'NORMAL  ^PPRQE.  if/CCNTINOITI  CORRECTION  FACTOR 

ENT  T  APPROXIMATION  VITB  CONTINOITI  CORRECTION  FACTOR 


CKr*(((X»(0p«2  )U(X-l))-( 

i  eWOttjfrfjlicf 0^  TO  AND 


- ( ( iNUM* 0 . 5 )*2  )♦ (X-l ) )  )*0 . 5 


gflj? Ki 

z:$l;!>^k/,h■Mi, 

’hBtiMM'lk'hgi  '‘oN  A  SAMPLE  SIZE  OP:  ’,(•») 


rCl))>*2 


M21 


in 

112 

113 

114 

.11! 

116 

Ilf 

[118 

S119, 
120, 
121 
122 
123 
124 
[125 

“■1  ''c 

■To'1 

[128] 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140, 

141 

142 

143, 

[144, 


;nrg  total  sum  op  positive  rapes  is:  ' . (•rposi ) 

■  IP  PAIRED  SAMPLE  TEST  GO  TO  £17  POR  OUTPUT  STATEMENT 

jrSS^/I^OS^^SO  j  f  a  <  ,(9MO),<  VERSUS  B Is  M  •  ,  (*£OC.rcCC;  1]  )  .  •  •  ,  ( *HO 

-cJa 

Ll?:'Z3E? -VALUE  POR  COMPARING  TBS  MEDIAN  OP  TBE  POPULATION  OP  ' 
'DIPPERENCES  TO  TBE  BIPOTBESIZED  MEDIAN ,< 

3C:  M(X-I)  s  \i»MO).’  7ERSUS  31:  M(X-I)  < . i •LOCIClC ; 1Z ) . ’  ',  (»*?).' .  ’ 
,(4*fK)|£C5l]  ) 

sou  like  a  conpidence  interval  por  tbe  median?  h/n).  < 

♦<SB*'T' )/£l6 
*0 

LlS-.CO+INPOT  * 

ALPBA- ( 1 0 C -CO ) ♦ 2 0 0 

«  ROUTE  TO  NORMAL  APPROX.  POR  CONP  INT  OP  LARGS  SAMPLE  SIZE 

~<MW>16)/£4 

a  COMPUTING  CONPIDENCE  INTERVALS  BS  EXACT  P- VALUE 

CVF’+*WILP  NNN 

«  .  ,  INDEX  POSITION  OP  CPP  POR  ALPBA  *  2 

TALPBA-(*/ ( CDPiALPBA ) ) 

-(TALPBA> 0 ) /JUMP 
TALPBA* 1 
-JUMP 

a  COMPUTING  CONPIDENCE  INTERVALS  BS  NORMAL  APPROX.  U/C.P. 

Lm:DBNOMZ-((2*NNNx<.NNN+1)x(<.2*NNN)+1  )  )  +  3  )«0.5 
^TALPBA*l ^NOHZxlNORMPTB  ALPBA) 1±InNNx{NNN*1  ) ) -2  )♦“ ) 

"  rAiBSl*  L  TALPBA  *  1 


NORMPTB  ALPBA) )*lNN 
\PBA  DOWN  TO  INTEGER 


)«-l))  +  3)*0.5  . 

UlNNNx (NNN+l ) )-2 )*4 ) 
TECER  VALUE  AND  INCREl 


■MENT  BS  ONE 


a _ _  JS  OW  SAMPLE  CASS  GO  TO  £7 

JUMP:*(AA*1 )/£ 7 

a  .  .  «£CO£S«  SSD  PATW  O0T  COW.  IW.  SOS  PWM5  SjU!P£S  OISS 

£  5 : Cl— T ALPS A  CQNFV  DD 

\xpM lmt*VXL  r0R  IBE  "EDIAN  0F  TBE ' 

;  (  (  ' . (»CJ[1] ) , '  S  MEDIAN (X-S)  S  ',(*C/[2]),»  )' 

♦0 

a  CALCULATE  AND  PRINT  OUT  CONP.  INT.  POR  ONE  SAMPLE  CASE 

£7 : CI+TALPBA  CONPW  X 

' A  ' , (»CC J , '  PERCENT  CONPIDENCE  INTERVAL  POR  TBE  MEDIAN  OP  TBE  POPULATIO 
If  IS: ' 

;  t  (  '.(•CJCl]).'  S  MEDIAN  £  ’.(wCIC 2]),'  )' 

♦0 

E2:> ERROR:  SOU  BAVE  NOT  ENTERED  A  VALUE  OP  1  OR  2;  TRS  AGAIN. ' 

♦52 

S3 : |  ERROR :  SOU  BAVE  NOT  ENTERED  A  VALUE  OP  1 ,  2 ,  OR  3 ;  TRS  AGAIN. ' 

♦S3 
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APPENDIX  F 


LISTINGS  OF  SUBPROGRAMS  BASIC  TO  BOTH  WORKSPACES 


•vdj 

a  TB 

’3' 

P+0 

X+0 

1 5  j 

CBN 

m 


il 

:8 

18 

lfl 

20 

21 

22 

23 

24, 

25 


CDF 


j 


A ' SUBPROGRAM  OF  TBS  SIGN  TEST  ( SIGN ).  IT  CALCULATES 
--'HAL  WHEN  PROBABILITY  =  .5.  N -SAMPLE - 


IE  SIZE. 


•  w 


E ! N ) x  (P* X) *  (  ( 1 -P ) *N- X ) 


7  U+L  CBIN  R 


A  SUBPROGRAM  OF  CONFIDENCE  INTERVAL  GENERATOR  FOR 
IT  RETURNS  THE  VALUE  OF  TEE  BINOMIAL  CDF  AT  R. 
SAMPLE  SIZE  AND  P-PR0BABILIT2 . 


j7-<*\<ri*ia«tfcij)i£U3  )*(ic2]*_ifti+£ci:  )x((i-lc2:  )*£Ci3-“i+u+£Ci:  ))w+ 

51 


7  CON*XX  CONFLR  Tl ;  BB  -,SS ;  AA  :A  ;XR  :2R;B:C:S 
A  _  ^BIS _ FUNCTION _ IS _A  SUBPROGRAM  OF _NONPARAMETRL 


a  (NPLR) 


LINEAR  REGRESSION 
IT ' CALCULATES' TBE~ TWO -POINT  SLOPE  FOR- EACH  PAIR  OF  POINTS 


a  RECORD  TBE  SIZE  OF  XX  AND  INITIALIZE  VARIABLES 
BB+pXX 

li:g° 

a  THIS  LOOP  COMPRESSES  XX  AMD  22  DOWN  TO  WHERE  TBE  XX  <  ALL  OTHER  XX <S 

TR+A/21 
B+pXR 
+(B=0)/L3 
C+  0 

a  TBIS  LOOP  CALCULATES  TBE  SLOPE  OF  EACH  PAIR  OF  PAIRED  DATA. 

LUQ+C+1 


•  C+C* 

S* (II CAA] -2RZC1 U (XXZAA1 -XRZC1 ) 


ss+ss.s 
+(C<B)/L 1, 
L3:+(AA<BB){L2 
l+(p SS). 


p SS),(SSZiSSl ) 


7  CONFM+AA  ^ONFHW  BB\A;B;C;D;E;FiC}B 


»  rj^AA  V  SStYt?  FAFT  W  •«  lO  •  W  •  V  lU  *  *  |  Vl  a 

TBIS  FUNCTION  IS  A  SUBPROGRAM  OF  TBE  MANN-WBITNEI  TEST 

OMPUTES  CONFIDENCE  IN. -  - -  - 

_ TION  BETWEEN  X  AND  2. 

A  BB-COMBINED  DATA  SAMPLES. 


_  _  _ _  (MANW). 

A  IT  COMPUTES  CONFIDENCE  INTERVAL  ENPOINTS  FOR  TBETA .  TBE  SHIFT  IN 
A  LOCATION  BETWEEN  X  AND  I.  AA-INDEX  POSITION  OF  C.I.  ENPOINT. 


A  ASSIGN  SIZE  OF  X  VECTOR  TO  il}  INDEX  POSITION  FOR  CONF  INT  TO  B. 
A+AAZt} 

3+AAZ21 

ASSIGN  X  VECTOR  TO  C\  2  VECTOR  TO  D 


C+A+3B 

D+A+BB 


KSEIS3 


REORDER  X  AND  2  VECTOR  VALUES  TO  ASCENDING  ORDER 


INITIALIZE  VECTOR  B  AND  VARIABLE  F 

B+  pO 
F*  1 

a INNER  AND  OUTER  LOOPS  CALCULATE  ALL  POSSIBLE  DIFFERENCES i  EVER!  2 
A  ELEMENT  MINUS  EVER 2  X  ELEMENT.  B  VECTOR  STORES  TBESE  DIFFERENCES. 
L2'.G+\ 

LUE+DZn-CZGl 
B+B  ,E 


98 


[23]  *7ci!pC))/£l 

[24,  1 

fll;  *"”(  *  1  VECTOR  VALORS  TO  ASCRRDZNC  ORDER 

Eh]  "r  m*'s  * 


COW  «?  71£crs  /AON 


»8f 


57AST  ff  ACCUMULATION  VECTOR  OFT  MITE  ORCIRAL  VECTOR  VALUES 


8*o  0 

£g.C 

a  OUTSIDE  LOOP  INCREMENTS  D  AND  RESETS  S  TO  D 
L2:D*D*l 

B*D 

A  INSIDE  LOOP  GENERATES  NEXT  SET  OF  A7 TRACES  AND  CONCATENATES  TO  ORIGINAL 
VECTOR 

L$*1cId]+CIE1  )*2 
3*8  ,F 

a  CONTINUE  INNER  LOOP  UNTIL  E  EQUALS  T8E  SIZE  OP  C 

*(E<(.oC))/LX 

a  CONTINUE  OUTER  LOOP  UNTIL  D  EQUALS  TEE  SIZE  OF  C  LESS  ONE 

*(D<UqC)-X))/L2 

a  ORDER  FINAL  ACCUMULATED  VECTOR  E 

8*8  LABI 

a  INDEX  CQNF  INT  VALUES  OUT  OP  8 


»  INDEX  CONI 

C0NP*8LA1  .aZ((o8)-U-l))l 


V  P*DF  FDISTN 


XAL-U ALLIS  TEST  {KRWL )  AND 
APPROX.  CUMULATIVE  PROBS. 
DEGREES  OF  FREEDOM. 


A  TAXES  VECTOR  ARC. 
A*Jt*X*N+XxH*2*DE  ,N*X+DP 

** TREAT* TBEN\ST.  FOUR  SPECIAL  CASES. 

Lli:P*l  loA»0 . 5  )*o*2 
♦0 

L12  :P*il*0 . 5 
♦0 

£2l:P*l-<W)*o.5 

*0 

L22'.P*A 

A  BEGIN , TEE , GENERAL , CASES .  INITIALIZ 
L:RH*Tl\H*2),0/RN*X\N*2 
SM*(RM=  0  )xRM£LM*\.M2*Q  .  S*M-2 
SN*lRN=0)»RNSLN*lN2*0 . 5*0-2 
A  TREAT  SN.CASE  N-EVEN  FIRST. 
*NODDx\l=i\N 
*MEVEN*\0=N2  „ 

SN*1 ++/x \ ( 1 -A )•.* 1 +<  1 +0 . 5 *M)* \N2 
*MEVEN 

NODD : *HEVEN* j  1  =J7 

a  S TREAT* IbE ' PORTION  OF  SN  TEAT  DOE SR 


INITIALIZE  TEE  QUANTITIES. 


*RMZLM*\.M2*0 .  S*M-2 
»RN£.LN*lN2*0.5*N-2 
.CASE  N-EVEN  FIRST. 


*S TREAT* IbE  PORTION  OF  SN  TEAT  DOESN'T  CBANCE  W/  M  ODD  OR  EVEN. 
*\2+X+QLC)*\0=LN 

SN*SNxl**'x\ (X-A)* .*\*2*M2*X+2*\LN 


,32J  a  TREAT  T3E  N-EVEN  SUBCASE 
*33*  *(  (M=2  ),2*2\M)/MEVEN.3*X*GLC 

k3u]  SN*SNxx/Z*+2*\ii2 

35  *MEVEN 

36  A  TREAT  TEE  M-ODD  SUBCASE 

37  SN*SNx 2*ol 

38  *END*\ X=M  . 

,39  SN*SNx*/l+*  l  +  2*»  0 .  S+JI2 
>0,  MEVEN:*M0DD*\\=2\M 
>1  *END*\2=M 

,>♦2,  SM*X++/x\A*.*l-*2*\H2 

43  *END 

>4  a  NEXT  TREAT  TBE  SPECIAL  CASES  FOR  ODD- M 
45  MODDiSM*(0*M-X)*{XA*0  ,5}+o*2l*3=M 
>6;  *END *i<l=#)v3=>l  , 

,47,  SM*CA*Q  .  5  )*  (* o*2  )*!++ /x\A»  .*1-*  1+2* \L1t 
[48]  END:P*LSN*A*M*2)*{'**RNxRN-RM)+URMxRN*  1 
7 


loA*0. 5  )*o*8  )-2*RN*SM*  (1-A  )«0  . 


XX*N XffOj 

S 

•  COMBINED . 


its or rai  samples  Btific  passed  in  b.  I*  all  data  samples 


SED  IS 


l l  data  Samples 


AA*oN 

*BB*C/N 


BBCOBD  SIZE  OP  N  TO  DETERMINE  NUMBER  OP  SAMPLES  IN  B 

LOCATE  LARGEST  SAMPLE  SIZE 
ORDER  B  SMALLEST  TO  LARGEST 


DD*B(kB 1 

%  SET  UP  MATRIX  OP  SIZE  REQUIRED  TO  STORE  SAMPLE  RANKING . 

xx*  (<u,  52)00 

•  CONCATENATE  0  AND  N 

S^’PIND  CUMULATIVE  SUMS  OP  SAMPLE  SIZES . 

NC**\NN 

SiBIS  LOOP  INDEXES  OUT  ORGINAL  SAMPLE  VALUES  FOR  TURTBER  CALCULATIONS. 
Ll:CC*CC*\ 

a  LOCATE  ;C'TB  ORIGINAL  SAMPLE  7 ALOES  IN  3 

»X1*aCi*C]iSSoNM$§IT}oNS  OP  ELEMENTS  OP  X  VECTOR  IN  B  AND  ASSICN  TO  C 
C*DD\Xl 


MrM 


»rblS  LOOP  DETECTS  TIED  INDEXED  POSITIONS  AND  INCREMENTS  TBE  INDEXINC  OP 
EACB  SOCCZSSI.lt  TIED  POSITION  37  ONE 
a  TSSi  DTB  ELEMENT  OP  C  AGAINST  REST  OP  C  POR  TIES 


L2:S~(c:D]*D*d 

a  SET  P  EQUAL  TO  TBE  APPROPRIATE  SI Z. 

P*DgO 

* pc$n%ATSNATS  P  AND  VECTOR  OP  0'S  AND  I'S  (1'S  APPEAR  WREN  TIE  OC CORED) 

n  ’  ADD  RESULTANT  P  VECTOR  TO  C 

C*C*F 
D*D  + 1 

«  CONTINUE  POR  ENTIRE  C  VECTOR 

*(Di  (of)  )/L2 

V 


TO  TBE  APPROPRIATE  SIZE  'DTB  SIZE)  ZERO  VECTOR 


V  I N+ 1  If  PUT  A  •  S  s  X  *  T 

J  PUNgribN  <fegTBE  INPUT  prompting  and  error  cbeckinc. 

s1gNB,DWISIG,  MANN,  KEN,  SPMAN,  AND  NPLR. 
L\\'\flTZkLThKbATA  MORE  TEAR  TWO  OBSERVATIONS  ARE  REQUIRED).  ' 


Ly/t^fiTzkLThK  b 

zCmm? 


LlilENTER  X  DATA  (.MORE  IRAN  TWO  OBSERVATIONS  ARE  REQUIRED ) .  ' 

♦TTpoX)30)/X1 
♦  (_(oX)*2)/ffl 

'ENTER  I  DATA  (NUMBER  OP  I  ENTRIES  MUST  EQUAL  NUMBER  OP  X  ENTRIES).' 

L3: ’ENTER  TBE  BIPOTBESIZED  MEDIAN.' 

IN*  0 

*((pIN)>l)/E3 

L*:' ENTER  TBE  BIPOTBESIZED  MEDIAN  POR  TBE  DIFFERENCES  OP  TBE  PAIRED  DATA.' 
IN*  U 

*((oIN)>D/E3 

*0 

L 5:'  ENTER  TBE  DESIRED  CONFIDENCE  COEFFICIENT :  ' 

i  f 

'  FOR  EXAMPLE:  ENTER  95,  POR  A  95  PERCENT  CONFIDENCE  INTERVAL . < 

El: 'ERROR:  TBE  SIZE  OP  JOUR  SAMPLE  IS  LESS  TBAN  TBREE ;  TRI  AGAIN.' 
*^A*l)/Ll 

E2:' ERROR:  SAMPLE  SIZES  ARB  NOT  EQUAL;  WANT  TO  TRI  ACAIN7  (I/N).' 


* ENTER' &ICBT  ARROW  *  TO  QUIT.' 


1  00 


[431  B*C 

,44  13:  TEE  BIPOTBESIZED  MIDI  AM  HOST  BB  A  SINGLE  VALUE  j  TXT  AGAIN .  ' 

5»:  :tf‘,)/£3 

,4f,  X4: 'ERROR:  TBIS  VALUE  HOST  LIE  BETVEEM  0  AMD  100 j  TXT  AGAIN. ' 

;jo;  *L S 

,|i,  Bi • \ ERROR:  TBIS  VALUE  HOST  BE  AM  INTEGEXi  TXT  AGAIN. < 


V  INTB 
*  f 
a  ||J»p 


JS  PASSED  AS 


RIGHT 


BE  IXUSEAL -WALLIS  ( KRWL ), 
PHAN  AND  SPMANXl .  WHEN 
ARGUMENT .  IT  CALCULATES  THE 
'THE  TABLE  VALUES  OP  STATS. 
T  ARGUMENT. 


a ^S^PAR+TE  TBS  CDP  TABLE  AND  STATS.  INTO  SINGLE  VARIABLES. 

nsc**)lu  uggftg  i  first  exceeds  ox  equals  one  or  TEE  table  values. 
^“Yr  rar  stat.  does  not  equal  ant  table  values  set  inter  *  ~i . 

*((*/Cl*0)/L2 

a  INDEX  LOCATION  OP  FIRST  OCCURENCE  OP  HATCH 

D*C  \  1 

a  IP  INDEXED  POSITION  EQUALS  ONE  INDEX  P-VALUE  OUT  OP  CG . 

*<D>1)/L1  , 

:ntsr*ggzdz 

a*°  OTHERWISE  CONDUCT  INTERPOLATION  TO  GET  PROPER  P-VALUE. 

G*GCtpZ-iGZD-\Z 
PL-lPxGi* I 
INTSR*GGCD-1]*PL  . 

ij: INTER*' 1 


■1 1*PL 


.  —  Wis  ruNctm •  n  X  kuBPm™»°r  UBbUP s  B  Ucen)  and  non- 

:immovM  mhsiMim-N.iT  calcMs  ibb  cuHOuh7E 

a1*  1 INITIALIZE  FREQUENCIES  POX  X  POX  SAMPLE  SIZE  NN. 

a  DETERMINE  SIZE  OP  RIGHT  PXOB .  TAIL  VECTOR 

C*(L  ((J»»(W-ll)*4l)+i 

a  Increments  tbroucb  the  n  sample  sizes  till  the  desired  one  is 

Ll:$*t>0 


0*0 

P*pX 

NN*NN+1 

B*((NN* 


ftiNN* (NN-1))*2)*1 

a  INNER  LOOP  GENERATES  NN* 1  FREQUENCIES  PROM  THE  VECTOR  OP  NN  FREQUENCIES 
L2:A*A-X 


L2:A*A-X 

aa*6 


OTHERWISE  CONTINU i 


NN  AND  INDEXES  OP  X  STILL  REMAIN  GO  TO  L 4 


_ ...  ISE  CONTINUE  to  increment  thru  old  preqs 

*((B-A)*F)/t2 

a  WHEN  LEFT  HALF  OP  NN*1  VECTOR  IS  COMPLETE  GO  TO  L 5 

*LS 

i  TRIG  LOOP  ALLOWS  1NLI  NN  TERMS  TO  3E  JSED  TO  GENERATE  N. 

/IA  *  <, 

£*£*xlIb-a)1-x:aa] 

Aui  (8-1  )*2  )♦!  )>  (B-A  ))/L<* 

a  Invert  vector  d  and  assign  to  e 

LUE*9D 

a  FOR  NN  OP  APPROPRIATE  SIZE  EITHER  DROP  FIRST  VALUE  OFF  E  OR 
*<M»*3,6.7,10,ll,14,lS,l8,19)/i3 

a3*j l:OHPjBTE  NEW  VECTOR  OP  PREQS  X  BI  CONCATENATING  D  WITH  E 

a  '  *  ’  CONTINUE  UNTIL  SIZE  OP  SAMPLE  N  IS  REACHED 

a*  GENERATE  VECTOR  OP  CORRESPONDING  P  STATS  OP  PROPER  SIZE  C 


3E  USED  TO  GENERATE  NEW 


TBER  DROP  FIRST  VALUE  OPT  E  OR  NOT 
3 


>4  »**  X kitCOLATJB  B  STATS  FROM  TBS  P 

45,  r*(((4*?2+7»>i(W- 1222*1' 

45  a _  TAKE  QNLI  6  ENTRIES  PROS  TBS  F 

47,  TPL*T*(Nx(N- 1))*2 

46  XX*GtX 

49  «„  CBARCS  FREQS  TO  CDF  VALUES  AND 

so  «ffp»(3,c)pr,((+\xx)*(ij»)5,rp£ 


ROM  TBE  P  VECTOR 

FROM  TBE  FRBQOSNCI  VECTOR 


OOTPOT  B  STATS  V! APPRO.  CDF  VALUES 


V  MAN+N  MANWP  M:P:QiPiS;TiO:ViB:NN:NiM:UO:MM 

a  TBIS  FUNCTION  IS  A  sbBpRdGRAM  Op  TpE  MANN-VBITNEI  TEST  (MAW), 

a  IT  GENERATES  TBE  CUMULATIVE  DIST.  FOR  TBE  U  STATISTIC.  N  =  SIZE  OF 
a  LARGER  SAMPLE-,  M  =  SIZE  OF  OTHER  SAMPLE. 

a COMPUTE  NUMBER  OF  TERMS  TO  BE  INCLUDED  IN  LEFT  TAIL  DISTRIBUTION  LESS  1 

MM*-(L(CWM)  +  2))*1 

a  SET  F  VECTOR  EQUAL  TO  1  CONCATENATED  WITH  MM  ZERO ' S 

P+l,MMp0 

a  SET  P  EQUAL  TO  TBE  MINIMUM  OF  N+M  OR  MM 

F*-l/((.N*M)  .MM) 

a  SET  Q  EQUAL  TO  TBE  MINIMUM  OF  H  OR  MM 

Q*l/M.MM 

a  66  TO  LINE  DENOM  IF  MM  IS  LESS  TEAM  N+USIZE  OP  X*l) 
*DENOM*\(.MM<N*-ll 

a  IF  MMiN+1  GENERATE  FIRST  BLOCK  OF  RECURSIVE  RESULTS  USING  NUM  LOOP 
NUM-.T+N+I 

1/2  *  #10*5 

a  PRIMARI  FORMULA  USED  IN  GENERATION  OF  FIRST  BLOCK  OF  RECURSIVE  RESULTS 
LI  :F7U1*-FrU •  -FrU-T"1 

a  assigns' nEv  Decremented  value  to  u  and  tests  if  t  <  iris  new  j 

~Ll*\(.T<U-U-l) 

+L2*\lPiT+T+t ) 

a  GENERATE  FINAL  RECURSIVE  RESULTS  USING  DENOM  LOOP 

DENOM: S*1 
L*:V~S+ 1 

a^  _  9rPRIJ£RJ  FORMULA  USED  IN  GENERATION  OF  FINAL  RECURSIVE  RESULTS 

■'♦£3*1 Um*i  '/zv+v+i ) 

»£4xi (fliS*S+12 

a  CONVERT  FREQUENCE  TABLE  TO  CDF  VALUES  FOR  FINAL  OUTPUT 

MAN*(.*\F)*(Ni\N*M)) 

V 


VZ*NORMCDF  X:A:BSC:D 

figures.  Ported  to  mainframe 

a  AVOID  mvEREZOV  PROBLEMS  WITH  *\ 
*({^A*yx<^)/Z*X*,X)sp,X)/3*QLC 


.  FOR  IX<4,  26.2.11 

i 


A*0 . 5+ ( *AX \ (*A*2 )*o2)*0.5 )x+/x\CA*2 )».*“l+2*\ (C+10rri0xr/|A)+l 

i[Yxnmn+A 

BlCsSi  16§89790  56295540  52050600  19934640  3680160  341952  15232  256 
D*  2027025  32432400  75675600  60540480  21621600  3843840  349440  15360  256 
B»l-(B*2x ( (02  )x*fl* 2 )*0 . 5  )*(  +  /( (0 . 5x ,B*2 )• .*0 , 1 7 )* ( (p ,B ) , 8 )pC)++/ ( (0 . 5* ,B* 


12  SW 


7  Z*NORMPTB  P :A  :B :C;D ; Q :T : S R :P 

»  iMLEMSJ^IS,AiG2Rll3MjSili6l  BEASLEJ  SPRINCER.  APPLIED  ST&T,  1977 

.JgyjJfggfe  --  32*  2UE  23  3233  -’fe Z22H- 


.  1977 
UAN- 


*<  v/ ( V ; Q+.P-l . 5 ,S3 . a 2 ) )/3 +CLC 

S+Z+.Q 


* EXT 

TWo.42ilC)/Z+,fl 

♦<M  ( p ,r]=p .p5 )/2+olc 

S*(0.42<Je)/Q 

A*  2-50662823884  18.6150006252  41.3 

fl*  8.4735109309  23.08336743743  21. 

r*rx( ( (r*2)».*o.i32+.*A)+i+( (r*2 )•. * 

zno.42iifl)/tp,fij*r 


18.6150006252  41.39119773534  25.44106049637 

3.08336743743  21.06224101826  3.13082909833 

)  +  . *A )  +  !  +  ( (T*2 )• . *  14  )♦. *B 


*(P=l)lQ 

EXT-.C*  2.78718931138  2.29796479134  4.8  501412713  5-2.321212768  58 

D*  3.54388924762  1.63706781897 

S4-(*S)*((B».*0.t3)'f.xC)*l  +  (((B*(l*0.5-|S)*0.5)».*  1  2  )  +  .*£) 


®3 

[22  IFF:' 


o.«»2<iG)/ipa3*s 


y-.'OSS  QS  ms  S  Z&ltBSS  ASS  qst  qs  SASQS • ' 


7  P*PERH  N-.X-.T-.Z 
a  TBIS  FONCtlbN  IS 
A  PERMUTATIONS  OF  N 
~Q*\ N*P*  1  1  pi 
Z+PSRH  N- 1 
J»*tX+0 

Ll:+0*\N<X+X+t 

I*(~(\N)tX)\Z 

p^\u*nin-i),n)p(.p 

*L1 

7 


a  TBIS  FOm 
a  PERMOTA  X. 


CALLS V  FT  SPBARP.  IT  GENERATES  ALL  POSSIBLE 
RANKS.  N*SAHPLE  SIZE. 


lxlN-l),N)e>(,P),.I 


7  M*N  QUANC  Q : I-.J:K:L iM; 0 

A  TBIS  FUNCTldN  GtvtS  A  CBOICE  OF  NONPARAHETRIC  CONFIDENCE  INTERVALS 
a  FOR  TBE  QTB  QUANTILE  OF  A  CONTINUOUS  POPULATION.  N -SAMPLE  SIZE, 
a  Q=QU ANVILS.  IT  CALLS  TBE  SUBPROGRAM  CBIN. 

a  CENTERS  CBOICES  AT  -TBE  ORDER  STATISTIC  NEAREST  ESTIMATE. 

I*lO.  5+N*Q*100 

a  I  IS  1.5  STANDERD  DEVIATIONS  (.APPROX.) 

J*'0. 5+1. 5* (1*1-4+100 )*0 . 5 
I~l+J*  1  1 

A  RQN  OUT  K  3  QTB  NATS. 

A  K CONFIDENCE  io\pFI Ci  E^S 


L*-/*(  2,  pF)o(JV  .0.01*0)  CSIAK  ( JCl]+F-l  )  .  (IC2]  -K)  ) 
M* ' ORDEk  STATISTICS  I  COEFFICIENTS ' 

A  PART  OF  FORMATTING  OUTPUT 
M*M,[0.l2  29o?(16o ),  » 7  *  1 

0*  J  0  S  0  ««(2,pJ:)p<lT<UC2]-2),-rCl>«-l)) 

0*0, ( (oF) ,4 )p 1  I1 
#*F.Cl!)  0 ,  12  6  «((p£),l)p£ 


)p  (1t(  (XC2 
i 

(p£),l)p£ 


7  SPN*N  SPAPROX  X-7 

a  TBIS  FONCTION  ts  A  SUBPROGRAM  OF  SPEARMAN ' S  R  ( SPMANP ) 

A  IT  APPROXIMATES  TBE  CUMMULATIVE  PROB  FOR  R  WREN  PASSED  TBE  SAMPLE 
A  SIZE  IN  TBE  LEFT  ARCUMENT  AND  TBE  ABSOLUTE  VALUE  OF  R  IN  TEE  RIGBT 
A  ARCUMENT.  SUBPROGRAMS  OP  TBIS  FUNCTION  INCLUDE'.  TDISTN 


A  CALCULATE  TBE  CONTINOITX  CORRECTION 
7*6 *Nx 1+0*2 

A  TRANSFORM  TBE  STATISTIC  R  INTO  ONE  TRAT  CAN  BE  USED  HITS  TBE  STUDENT 
A  T  DISTRIBUTION  %  % 

X*(X-3Ul<.N-2 )*1-Qr-T)*2)»0.5 
A  CALL  TBE  T  DIST  FUNCTION  TO  CALCULATE  TBE  P -VALUE 
^PN*\-(N-2)  TDISTN  X 


7  SPEAR*SPEARP  NsCl:A:B:C:D:E:M:NiDliD2:LIM:R:CDFt 
A  TBIS  FUNCTION  IS  A  sbBpRbckAH  bP  SPEARMAN'S  R  (.SPMANP).  IT 
a  CALCULATES  TBE  EXACT  CUMULATIVE  DIST.  FOR  R  FOR  TBE  SAMPLE  SIZE 
A  PASSED  AS  TBE  RIGBT  ARGUMENT.  BECAUSE  OF  TBE  LARGE  COMPUTER  MEMOR1 
A  REQUIREMENTS ,  N  IS  LIMITED  TO  SIX  ON  TBE  PC  AND  7  ON  TBE  MAINFRAME. 
a  SUBPROGRAMS  tALLED  B1  TBIS  FUNCTION  INCLUDE-.  PERM 
A*  B*0 

a  INITIALIZE  VARIABLES ,  VECTORS,  AND  MATRICES . 

Cl*01*oO 

M*i.N,N)oO 

c* o.i (n-i ) 

a  TBIS  LOOP  GENERATES  AN  :V*V  A.RRAT  OF  TBE  POSSIBLE  VALUES  OF  DIFFERENCES 

A  BETWEEN  ANT  TWO  PAIRED  RANKS  BETWEEN  SAMPLES. 

Ll:F*F+l 

M(B-,)*C-A 

4*4+1 

*(B<N)/L1 

a  NOW  CALCULATE  TBE  SQUARES  OF  ALL  POSSIBLE  DIFFERENCES . 

M*M*  2 

A  CALL  PERM  TO  LIST  ALL  POSSIBLE  PERMUTATIONS  OF  N  NUMBERS 
D*PERM  N 

a  CALCULATE  SIZE  LIMIT  OF  FINAL  VECTOR  OF  R  STATS., 

LlM*l  +  U\N*3l-N)*l2l 

A  INITIALIZE  VALUES  BEFORE  INDEXING  OUT  COMBINATIONS  OF  ALL  POSSIBLE 


PASSED  AS  TBE  RIGBT  ARGUMENT.  BECAUSE  OF  TBE  LARGE  COMPUTER  MEMORT 
REQUIREMENTS .  N  IS  LIMITED  TO  SIX  ON  TBE  PC  AND  7  ON  TBE  MAINFRAME. 
SUBPROGRAMS  tALLED  BT  TBIS  FUNCTION  INCLUDE:  PERM 


A  SQUARED  VALUES. 

B*(N,'.N)oO 

a  TBIS  LOOP  CALCULATES  ALL  POSSIBLB  COMBINATIONS  OF  TBB  SQOAP.ED  VALUES . 
L2:A*A*1 

n*ADD  DOWN  ALL  ROWS  FOR  BACB  COLUMN  TO  SUM  UP  SQUARES  COMBINATIONS. 
D2*WE 

a  ADD  UP  NUMBER  OF  DUPLICATED  SUMS  OF  D-SQUARED  VALUES  AND  COMPRESS 
a  VECTOR  DOWN  TO  UNIQUE  VALUES. 

L3:C1*C1.(L/D2 ) 

Di*Dl,U/l(L/D2)=D2)) 


°im 


A  TRANSFORM  SUM  OF  SQUARES  VALUES  TO  SPEARMAN'S  R  STATISTIC. 
R*1  -  ( 6 *C1 ) ♦  ( JVx  (  < N* 2  J  - 1 ) ) 

a  CALCULATE  CDF  VALUES  ASSOCIATED  WITB  TBE  R  STATISTIC. 
CDF*(+\D1)*  '.N 

a  FORM  TWO  ROW  MATRIX  FOR  OUTPUT  OF  R  STATS  AND  CDF  VALUES. 
SPEAR*  ( 2  ,  ( p  Cl ) )  p  <J? ,  CDF ) 

7 


7  SPM*X  SPMANP  T;C;D;DD;Dl:D2;N;DEN0MR;XXi71;NS;NUMR;PiPV;PVAL:SU;SV;RB0 

n .  y .  n  duq  •  .  y^  .  x .  y .  jyy 

A  ’  ' THIS ’ FUNCllof/  Is  A  SUBPROGRAM  OF  NONPAR  LINEAR  REGRESSION  ( NPLR ) 
a  AND  SPEARMAN'S  R  (SPMAN).  IT  COMPUTES  TBE  SPEARMAN  3  STATISTIC 
a  AND  ASSOCIATED  ? -VALUES.  TBE  LEFT  ARGUMENT  TBAT  IS  PASSED  IS  TBE  X 
a  SAMPLE i  TBE  RIGHT  ARGUMENT  IS  TBE  T  SAMPLE. 

a  SUBPROGRAMS  OF  TBIS  FUNCTION  INCLUDE:  TIES,  TTESK,  SPEARP ,  SPAPROX . 
a  INTERP,  AND  TBE  VARIABLE  PMATSP . 

A  ORDER  I  IN  INCREASING  ORDER  OF  X 


V*IUX1 

>o*xzbxi 

i 

(j 


ORDER  X  IN  INCREASING  ORDER 
COMPOTE  CURRENT  RANKING  OF  X 
NOW  ORDER  I  RANKS  IN  INCREASING  ORDER 


D1*VZ6V] 

A  IF  TIES  EXIST  IN  EITBER  X  OR  7  RANKED  VECTOR  USE  MID -RANK  METHOD 
DD*  1  TIES  D 1 
XX+1  TIES  U 

A  FIND  ORIGINAL  RANKING  OF  7  WITB  TIES  RESOLVED 

ll*DD£Cl 

A  RECORD  SIZE  OF  INPUT  VECTOR 

a  11  CALCULATE  DIFFERENCES  BETWEEN  RANKS  OF  X  AND  7  VECTORS 
D+XX-X 1 

A  DETERMINE  THE  SUN  OF  SQUARES  OF  TBE  DIFFERENCES 
D2*+/(D*2) 

a  OBTAIN  TBE  NUMBER  OF  TIES  IN  EACH  VECTOR  USING  TEE  TIESK  FUNCTION 
U\*TIESK  U 
V\*TIESK  DA 

SU*U+/lUl*3})-t+/UAU*12 
SV*U*/lVA*3))-l*/Vl))*12 
NS*N* ( (N* 2 )-l5 

A  CALCULATE  TBE  R  STATISTIC  INCLUDING  TBE  CORRECTION  FOR  TIES 
NUMR*(NS)  +  a  S)*D2)  +  a  6>(SC/+SV)) 

DENOMR*UNS- (12*SU)  )«0 . 5  )*  (  (.NS- (12*SV)  )*0 . 5  ) 

RBO*NUMR* DEN OMR 
ARBO*\RHO 

a*<J726  /l1CALIi  SpEARp  F0  CALCULATE  TBE  RIGHT  TAIL  OF  TBE  CDF  OF  R 
P*SPEARP  N 
*L2 

LA:*(N>A0}/L3  „ 

P*PMATSP( (N-5): :] 

A  CHANGE  SIZE  6P  P  TO  AN  M*N  MATRIX 
P*PZ 1: s* 

a  CALL  INTERP  TO  CALCULATE  P-VALUE  B 7  INTERPOLATION 

L2:PVAL*&RH0  TNTERP  ? 

-<.PVAL-  ii/L7 
PVAL* 0 . 5 
*C  7 

A  CALCULATE  P  VALUE  USING  STUDENT  T  APPROX . 

L3:PVAL*N  SPAPROX  ARHO 

LI :SPM*(RBO) , PVAL 
7 


7  P*K  TDISTN  X: 7 
a  IBIS  FUNCTION  IS  A  SUBPROC. 
B  FOR  SPEARMAN ' S  R  JSPEARP) . 
b  THE  STUDENT'S  T  VIST  HITS  . 
a  THIS  FUNCTION  CALLS  ON  TBE 
V*(X*,X)iO 

P*-O.S*(lJt)  FDISTN  X*2 


SUBPROCRAM  OF  TBE  CUMULATIVE  PROBABILITY  CSNERATOR 
PEARP).  IT  CALCULATES  TBE  CDF  AT  X  USING 
T  WITS  K  DECREES  OF  FREEDOM. 

ON  TBE  1  ?'  DISTRIBUTION  FUNCTION  (.FDISTN). 


P* O.S*(l  JO 
♦8*io=v/r 
plv/  tpjn+o.i 


PC7/tpX3+0.5+7/P 
»0*  l  1  =  a/7 

PC(~7)/ip*>0.S-(~!O/P 

7 


7  TI*BB  TIES  B:C:DiI:NiTiIiZiK:M:L : PP : NR :MM 
B  THIS  FUNCTiof/  is  A  sbBt>RbckAk  bF  kENb ALL'S  B  (KEN).  SPEARMAN'S 
a  R  ( SPHANP).  XRUSXAL-WALLIS  ( KRWL ).  MANN -WHITNEY  ( MANW ) 
a  AND  VILCOXbN  (WISIG) .  IT  CHECKS  THE  RIGHT  ARC.  VECTOR  FOR  TIES  AND 
a  CHANCES  TBE  TIED  POSITIONS  OP  THE  LEFT  ARC.  BY  TBS  MIDRANK  METHOD. 

N+cB 

A  IF  NO  VECTOR  OF  RANKS  IS  PASSED ;  CENERATE  ONE 

*UppBB)*0)/LB 
BB*\N 
L6 :Z*pO 
L+NpO 
T*  1 

a  CHECKING  FOR  TIES  BY  INCREMENTING  THRU  TBS  VECTOR 

LB  *  **  ('w+3)2Sr'T,*' 

a  '  COUNT  NUMBER  OF  TIES ;  IF  NO  TIES  GO  TO  L 2 
D**/C 
*(D=0)/L2 

a  RECORD  WHERE  TIBS  STARTED  AND  BOW  MANY  RANKS  INVOLVED 
I*I,I,(D*\ ) 

A  INCREMENT  NEXT  T  BY  TBE  NUMBER  OF  TIES  ENCOUNTERED  PLUS  1 

L2:T+T+(D±1) 

A IF  T  LESS  THAN  SIZE  OF  ORIGINAL  VECTOR  GO  TO  L 3  AND  START  ACAIN  AT  NEW 

-( T<N)/L2 
Y*pi 

a  ASSICN  TBE  RANKS  OF  TBE  LEFT  ARC.  TO  TI 

TI*BB 

A  IF  NO  TIES  FOUND  QUIT 

*(Y=0)/0 

z*o 

a  LOCATE  TBE  INDEXED  POSITIONS  OF  TIED  RANKS 
LS:PP*(  ( jTl+ZJ j-l ( (iTl+Zj +J[2+Zj')-1  ) 
a  FIND  THE  MIDRANK  VALUE  OF  THESE  RANKS 
NR+1+/TIIPP1 )+I[2+Z] 

A  SET  UP  VECTOR  WITH  ZEROS  AND  ONES i  ONES  WBERE  TIE  RANKS  INVOLVED 
K*0 
MM-t-iVpO 
L*:K*K+1 

M*(r/=ncppcjm ) 


A  LOCATE  TBE  INDEXED  POSITIONS  OF  TIED , 
£5:pp-»aiti+zl)-i)+i(aj!i+z]+j[2+ziVi: 
a  FIND  THE  MIDRANK  VALUE  OF  THESE  RANK i 


HM+HM+M 
*{k*IL2+Z1)/LU 
b  SET  UP  A  VECTC 


A  SET  UP  A  VECTOR  WITH  ZEROS  WBERE  TIED  RANKS  OCCUR 
L*~MM 

A  TRANSFORM  ONES  OF  MM  VECTOR  TO  MIDRANK  VALUE 
MM*NR*MM 

a  TRANSFORM  ONES  OF  L  VECTOR  TO  REMAINING  UNCHANGED  RANKING  VECTOR 
L*TI*L 

A  FILL  IN  MIDRANK  VALUES 
TI*L*MM 
Z+Z+2 

A  DO  TBE  SAME  FOR  ANY  OTHER  TIES  INVOLVED  BUT  WITH  NEWLY  COMPUTED  TI 
*(Y>Z)/L 5 
7 


7  TIE-TIE SK  AA:AA:B:CiD;I:N:T 

a  T HI 5  FUNCTION  IS  A  SUBPROGRAM  OF  ".END ALL'S  3  (REN'.  SPEARMAN'S 

a  R  (SPMAN)  AND  .SPMAN 1).  AND  KRUSKAL -WALLIS  .KRWL).  FT  ORE  IKS  IRE 
a  .'.TOST  ARGUMENT  FOR  TIES  AND  RECORDS  FEE  CUMBER  OF  OCCURENCES  OF  IACB 
fl  TIE  AND  THE  TOTAL  NUMBER  OF  TIES  IN  IRE  VECTOR. 


ASSICN  ORDERED  VECTOR  TO  B  AND  INITIALIZE  VALUES 


10]  I* 1 


B+AAtbAAl 
TIE* pO 


L3:C+(T*B)*BlT 


D**/C 


'BACKING  FOR  TIES  BZ  INCREMENTING  TBRO  TRE  VECTOR 
•otlNT  NUMBER  Of  TIES’.  IF  NO  TIES  GO  TO  11 


IS  *<5sO)/L2 

it; 

IS  L2’.r+T+(D+il 
19  *IF  T  LESS, TB AN  , 


•NT  NEXT  T  BZ  TBE  NUMBER  OF  TIES  ENCOUNTERED  PLUS  1 


tIF  T  LESS  TBAN  SIZE  OF  ORIGINAL  VECTOR  CO  TO  L 3  AND  START  ACAIN  AT  NEW  T 
*(.T<(.N-l))/L3 
1 


V  VAR*VARMV  B:C;D:E’.Dl;Sl 

a  TBIS  FUhdtlbN  IS  A  SUBPROGRAM  OF  TBE  MANN-WBITNEZ  TEST  W). 

•  IT  GENERATES  TBE  RANKING  SCBEME  USED  IN  CALCULATING  TBE 

a  DIFFERENCES  IN  SCALE  (1  ASSICNED  SMALLEST.  1  ASSICNED  LARCEST .  3 
a  ASSIGNED  NEXT  LARGEST,  u  SECOND  SMALLEST , Pi C .  31  TWOS  TILL  PP.tPSR 
a  SAMPLE  SIZE  IS  REALS E 5).  SAMPLE  SIZE  IS  PASSED  IN  TBE  RIL3T  ARC. 

B+E l*p 0 
5*0 

a  FIND  FLOOR  OF  MIDPOINT  OF  VECTOR  AND  ASSIGN  TO  C 

C+LCS+2) 

•  LOOPS  GENERATE  RANKING  VALUES  LEFT  HALF  FIRST 
L2 :5*5+ 1 

B+E.D 

*UoE)=C)/L3 
5*5  +  3 
B+-S  ,D 

+  U3S)<C)/L  2 

a  NOW  GENERATE  RICBT  BALE 

-3:51*1 
51*51+1 
E\*E\ .51 
*t(oEi  )=C)/55 
16:51*51+1 
£l*ffl,5l 
*((oBi)-C)/L3 

51*51^3 

21+51,51  » 

♦<(pri )<C)/L 6 

a  IF  SIZE  OF  VECTOR  IS ’ODD  VALUE  CONCATENATE  MIDDLE  RANK  IN  BETWEEN  BALES 
£3:+((2m*0)/57 

a  IF  SIZE  IS  EVEN  CONCATENATE  LEFT  BALE  WITS  TBE  REVERSE  OF  TBE  RICBT 
VAR*E.l»B 1) 

♦0 

L7-.VAR+E.B.  CE1  ) 

V 


V  VIL-WILP  NN \N :A:P :T :NN \W :PP :NM 

a  TBIS  FUNCTION  Ip  A  SUBPROGRAM  OF  TBE  WILCOXON  SIGNED  RANK  TEST 
a  IWISIC).  IT  GENERATES  TBE  CUMULATIVE  DIST.  FOR  TBE  TEST  STATISTIC 
a  TBE  GENERATOR  USES  A  RECURSIVE  FORMULA.  NN-SAMPLE  SIZE. 
J»«*(LU+/iAW)t2))+l 


J»«+(l((+/iAW)  +  2))+l 
N-2 

a  SET  P  EQUAL  TO  PROS.  DIST.  WBEN  N  EQUALS  2. 

P*  upl 
L3:N*N*1 
A*1 

a  SET  T  VECTOR  TO  PROPER  SIZE  OF  ZEROS. 

T+(+/(\N))pO 

a  IF  ASN  USE  TRUNCATED  RELATION  TO  COMPUTE  OCCURRENCES. 

L*i  +  UiN)/Ll 

a  IF  A>N  USE  FULL  FORMULA  TO  COMPUTE  OCCURRENCES. 

*(A>N)/L 2 

a  WBILE  Ta-N )  IS  NECATIVE  TRUNCATE  FORMULA  TO  AVOID  A  NEGATIVE  INDEX. 
LX’.TlA-\*PlA1 

A-A*  1 
*5“ 

a  IF  A  IS  LARCER  TBAN  TBE  LENCTB  OF  P  GO  TO  56 . 

L2:*(A><s?) :/L 6 

a  ONCE  ’A-N)  BECOMES  POSITIVE:  TBE  RECURSIVE  FORMULA  CAN  3E  USED. 

ClAl  *~P  [  \Amii  I  j  *P  CA  l 

a  ONCE  A  IS  LARCER  TBAN  TBE  LENCTB  OF  P  TRUNCATE  FUUCTION  AGAIN. 
56:rU>PCU-JO] 

L7:A+A+1 

a  IF  A  AS  AN  INDEX  BAS  NOT  EXCEEDED  N(N+ l)/2  GO  ACAIN. 

+  (Ai(.+n\N)))/L'i 

a  CONVERT  T  INTO  P  AND  CONCATENATE  1  FOR  USE  IN  NEXT  ITERATION  OR  OUTPUT. 


»  IF  A  AS  AN  INDEX  BAS  NOT  EXCEEDED  N(.N+ 1)/2  GO  ACAIN. 
*(Ail*/(\N)l)/L  i* 

»  CONVERT  T  INTO  P  AND  CONCATENATE  1  FOR  USE  IN  NEXT  ITERATION  OR  0U 
lS:P+((.  +  /(\N))pT)  ,1 
PP+NM+P 

»VILCBECK  IF  LENCTB  OF  INPUT  VECTOR  EXCEEDS  NUMBER  OF  N<S  GENERATED. 
*(,NN>N)/L3 


I,1  i, 


APPENDIX  G 

LISTINGS  OF  PROGRAMS  USED  TO  GENERATE  C.D.F.  COMPARISON 
TA3LES 


EENTEST :  N :  ALP  BA :  S  :  A  :  C ;  TAD ;  P ;  NUM ;  DEN ; 
■,rSiS;J;EkRZZCitc-.Kk 


Z  :  22  \ERRZZ  ;M;ZZC;D;AA  iPP;NUMC;I  ;B;F 


aTSIS  PROCRAM  ~ZNE?.ZTS:  TABLES  OF  C.D.F.  COMPARISONS  FOR  RENTALS '  S  3 
a  SET  SAMPLE  SIZE  AND  ALPBA  VALUES. 


: :  . .  - :  :  .  . :  : . :  ?  : . :  a  ; . : :  : . :  3 

raps  pocp  increments  sample  size. 


£1:PP-oC 

#*o0 

XJt-oO 
PS-o  C 


a  COMPUTE  CUMULATTVE  LIST  AND  ASSOCIATED  STATS. 

P-KENCALP  V 

a  TS.'S  LOOP  CALCULATES  ALPBA  7ALUES  AND  APPORXIMATIONS . 

£2 : 5* 2*1 
C-ALFSA'S: 

A*'*'  Plz  '.sc : : 

TAVzV-ii- 

XX*  XF.,  I  AC 

AA*6*N «  1*W»2 
PP-PP.P:2;X: 

a  COMPUTE  NORMAL  APPROXIMATION . 

NCJM*3*  (TAU )*  (  (JV*  (tf-1  )  )«0 . 5  ) 

DBN*(.2*  (2*/»)-5  )  )«0.5 
Z*NUM*DEN 

ZZ*ZZ  Al-NORMCDF  Z ) 

a  COMPUTE  NORMAL  APPROXIMATION  WITS  CONTINOIT7  CORRECTION. 

NUMC*3*(.TAU-AA)*{  (.N*  («- 1  )  >0 . 5  ) 

ZC*NUMC*DEN 

ZZC*ZZC .  ( 1  -NORHCDP  ZC) 

~(B<7)/l2 

a  COMPOTE  ERROR  DIFFERENCES . 

?RRZZ*PP-ZZ 
RRZZC*PP-ZZC 

a  pfljwr  oar  rxpir  or  radars . 

1*0.1 

.7-195 

5* 1  TEST  ST AT  VALUE  1 

M*H,'  '  PROBIB  2  fi];  FOR  SAMPLE  SIZE  EQUAL  TO  ',(2  0  *N),' 

*L3 

£5:P-19u 

L3:M*M.LI1  CU7CC 18o 197) ,.7.6l0<80l97),.7] 

*(£=>*  )/H* 

B*  1  18  oS 
-(P*0)/£6 

F*>K^G<\  9.9999  >'  GPMT(1  7  oKK) 

F*  1  62  pF 
FS*B,F 
*11 

£6 : PS-oO 
ff-17oS 
:'S*'S  .  3 


F*'  I  '  .  (7  5  aPPCC] ) 
PS-PS.P 
*lC<7 J/£9 
L 7 : M*M ,  [  1 J  FS 
D*D*  1 
P-1 

*(tfl.r2.P3,£5)C0] 
rirs-'mcr  c.d.f. 

J*  198 
*L3 

N2:S*' ERROR;  NORMAL 
PP*ERRZZ 
*13 


1  07 


t.  i.  -u  -ta  .k'-U  err? 


S3: S*< ERROR;  NORM.  N/CC  ' 
PP*ERRZZC 
*LZ 
£4:J* 

*(N<l<t)/Ll 

7 


7  KWTSST  -.A-.B-.C-.  PP  -.PC :  AW ;  t J» :  K R :  RVALUE ;  PV»£  ;  ».»:  ;  PF ;  =7F  ;pVFl  ;C  ;  r.W 

fl  S  ERR? ; ERRFl ; CDp ;CC:klC;I  iJ;FF;S:M 

a  THIS  PROGRAM  GENERATES  TABLES  OP  C.D.P.  COMPARISONS  TOR  TBE  ZRUSKAL- 
A  WALLIS  TEST. 

A*  4  4  4 
CC*  2 

ALPRA *  0.01  0.02  0.03  0.05  0.08  0.13  0.18 
*L  2 

Ll:A*NN 

a  CALL  KROWALP  TO  GENERATE  EXACT  DISTRIBUTION  FOR  SELECTED  SAMPLE  SIZE. 
L2:P*KRUWALP  A 
8*0 
N*+/A 
CC*CC+ 1 
K*pA 
M*p0 
PP*  oO 
XX*  0  0 
PC*  oO 
PF*  o  C 
FBI-00 
CDF*? 12; 2 
£4:B-»-B+l 

b  LOCATE  POSITION  OF  EXACT  CDF  7 ALOE  &  ALPRA . 

D*  (.+/(. CDFZALPBA LB2 ) ) 

*<0*0)/£3 

£>♦1 

a  DETERMINE  CORRESPONDING  TEST  STATISTIC  7 ALOE . 

L2:B*?:i;D2 
KK+KK.  3 

a  RECOkD  EXACT  VALUE  OF  CDF. 

PVALOE*PX 2  :D2 

pp*pp  .pvAlIje 

A  COMPOTE  CORRESPONDING  P-VALOE  USING  CBISQ  APPROX. 

PVAL*{K-ll  CBISQ  B 
PC*PC‘(1-PVAL) 

nCOMPutE  CORRESPONDING  P-VALUE  USING  F  APPROXIMATION 
F*CiN-K)*B)*  <  C/C-l  )*<  {N-l  )-B) ) 

PVP*UK-llAN-K2)  FDISTN  F 


PVP*UK-1),(.N-K2)  FDlSTN  F 
PF*PF. (l-PpF) 

A  COMPUTE  CORRESPONDING  P-7ALUE  USING 
Pm+UM).((W)-D)  FDISTN  F 
PF1*PP1. (1-P7F1) 

»(B< 7)/i4 


F  APPROX  V/  1  LESS  D.F.  IN  DENOM 


D*  0 

ERRB*EP-PC 

ERRF*PP-PF 

ERRP1*PP-PF1 

a  PRINT  OUT  TABLE  OF  VALUES. 

1*0.1 

J*'-< 

S* ' TEST  STAT.  VALUE  ' 

M*M,  1  PROBIH  i  <  .CU7C46]  .  ']  :  FOR  A  CROUP  OF  3  SAMPLES  CONSISTING  OF  ** 
,  4,  AND  ' , («CC). '  OBS.  ' 

*L  8 

L2:J*'-' 

L8 :M*M .112 (17p ) ,J,61p (8p ' - ' ) ,J 
*(.D=5)/L8 
F*  pO 
fl-16pS 
P*F.B 
C*  0 

*(D*0 )/£l0 
£11: 

r?t'  -t,'  ■  '7  u  vZL'.C'.  ) 

i(C<7)/£U 

+L12 

£10:C+C+1 

FF-'  I' , (7  5  9KKLC2 ) 

F*F,FP 
♦(C<7 )/L10 
£12 :H*M, Cl3  F 
D*D*1 
1*1 

*(.N1,N2.N2,NU.LS)ID2 
N1 lS*' ALPHA  VALUE  ' 

KK*PP 
J*'  +' 


£gSUMIM>  CtISQ0JLU' 
#1:5*  ERROR-.  r  r  1ST 


pj  "  1  "■ 

ii]  *Ii 
i;  ;j:# 


ao!  -^.rco o: 

v 


■.U-.a-.ir-.  PVAL  ;PPl.BB;Mt;*-.X  I EPRC  :  ERRFl  ;  f 


m  *f  *  *  r/?  ■*  J  ~JIM  JrVr^  -  ~  r  '  ~  f4  ~  r  r?"  ^  '  *  "r*- 

•  iisfr  jit  'i::oc"*l)cd#~ vtci'oks  or  •z"luii'iKcs"uh::a  "s'imvlate  roitfa.*:*: 
«  «  SAMPLES  TP  «  trsertatioms  EASE 

4*1 

#*ioocc 

5i*o° 

#**»0 

»#*' 

.-■i 

<W*- 

a  TEUERATE  «  lANDCM  'SOTOS  Of  SAMEL SOS. 

4**o,>*: 

77*0  3 

•  SOTf  rtr  .340*5  OP  fid  54JIPd. 

:::»**  ■*#. .«♦*;: 

*T^  *»  * 

i*c  * a ’ 

■^rI2caiif?2fsl  s  srarosTrc  pop  0400  yrrro*  op  ramkimos 

■  *^a5c)  °&UrDO/i ' S£^£  * 1°30:  ITERATIONS. 

OK*  C  *0  *C  1CCC  a  *<♦175 

L»:-(NN*V)/Ll 

■  jRD 4gf  ff'S  700#  LARCEST  TO  SMALLEST. 


•S  PROM  LARCEST  TO  SMALLEST. 


AL^BA*  C.C1  C  .  02  0.03  C.CS  0.01  0.13  C.H 


«•»<  *<  . 

Z|*oC 
?Pl*oC 
>C*eO 
13 : 1*1*1 

■  DETERMINE  crsdfr 
B^ai\(.RHM  ALPS  Alt.  j; 


0  74~0r  CORRESPOMCIHC  TC  ALPBA  VALVE. 


EX~Kk,B 


ess  z.r.  :u  venom 


PVPl-((t-i).(CR- 
pfi-ppi^.i-prp: ) 

ri0'*ic?04 -pc 

ERRPl-ALPBA-Pr: 

SB  *  »* 


P0.*0f  oot  rasor  op  valors. 


^P  •  OSS.  t-Oa  . 

Li:J»- -■ 

£•:#*#, t  ne 1  -  •  ,  .o.eic 
*(5**. '  72S 


c*c 

• (5«C  )  Lie 
L 11 : C^v ♦ 1 

.  V 7  3 


w  W  _  ) 


*”  <* 


?! 

8 

?s 

is 

:«i 

•  5 

IS 

:«o 

.31 

•flu1 


*(C<7)/Lll 

•£l2 

U0:C*C*1 


r*> 


r*P.rp  '(7  5  *“cc3) 

6+P+i 


&'! hSHUH1 


H+alpSa 

j+'  + 


'iihw0*'  '«««"*• 


-;i 

23:il 


O-  20000  30000  1 
V*?*2 ; 


V 


E  3 


tl 

■! 

i6] 

11 

12 
13 
1«* 

15 

16 
12 

11 

12 
20 
21 
22 
23 
2“ 

25 

26 
22 
21 

29 

30 

31 

32 

33 
3“ 

35 

36 

37 

38 

39 

40 
“1J 

u?3 


7  MANTES!  :0:NN  :MM :  ALPHA  :3:A:C:  !:?:NUM :  VEN :  :  :  ZZ  \ERRZZ  -M-ZZC :  DiPPlNUMCi 
?  •  *  S  :  J  : :  TT :  * ;  PTC :  TC :  NUitC ;  NUVli ;  t/ENOM ;  OSNOMC ;  SRnA  Vi  :  iRP.TT ;  iP.RTC ;  i.-.R. 
C-.ERRTOZC-.IZ, 

a  THIS  PROGRAM  GENERATES  TABLES  OP  C.D.P.  COMPARISONS  FOR  TRZ  MANN- 

«  nbItnet  test. 

SET  SAMPLE  SIZE  AND  ALPHA  VALVES. 


NN+7 

MM*7 

ALPHA*  0, 


01  0.02  0.C3  7.05  0.08  0.13  0.18 
THIS  LOOP  INCREMENTS  SAMPLE  SIZE. 


COMPOTE  CUMULATIVE  VIST  AND  ASSOCIATED  STATS. 

P*NN  MANUP  MM 

•  THIS  LOOP  CALCULATES  ALPHA  VALUES  AND  APPORXIMATIONS . 
L2:B*B+1 

C+ALPBAIB1 
A*(+/(.PiC n 
0-i-l 
KX+KK.U 
PP*PP.PIA1 

•  COMPUTE  NORMAL  APPROXIMATION . 

NUM*U-((NN*MM)+2) 

DEN*UNN*MM*  ( NN+MM+1 )  )*12  )*0 . 5 
2 *NUM+DEN 
22*22  NORMCDP  Z 

•  ’  COMPUTE  NORMAL  APPROXIMATION  WITH  CONIINUITI  CORRECTION. 
NUMC*(U+0.i)-<.(.NN*MM)  +  2) 

ZC*NUMC+DEN 
ZZC+ZZC .NORMCDP  ZC 

•  COMPOTE  STUDENT  T  APPROXIMATION 
DENOM*  ( ( INN+MM- 1 )  -  ( Z«2 ) )  ♦  (.NN+MM-  2  )  )»0 . 5 
l*Z+ VENOM 

.r-rr.  nn+mm-z  '  totstn  * 

«  COMPOTE  CTO LENT  0  APPROXIMATION  WITH  CONTINUITT  CORRECTION. 

VUMT-  COM  *•.  .  . 

LENOMC+.  ... NN+MM-  1  ;«  t.LEN+2  i  )♦  i.NN+MM-2)  >-  v  (  '.NUH+C  .  5  )  +  2  )  +  •.NN+MM- 2  )  )  ;*C  . 

tc+numt+vEnomc 

T^+T^C^<.NN+MM-2)  TDISTN  TC) 

•*  COMPUTE  ERROR  DIFFERENCES. 

ERRZZ*PP-Zl 
■RRZlc+PP-ZZC 
RRTT+PP-TT 
'RRTC+PP-TTC 
•RRAVE+PP-X\ZZ+TT)+2) 

ERRTCZC+PP-UZZC+TTCI+  2  ) 

PRINT  OUT  TABLE 


I*  3.1 
•7*195 


OP  VALUES. 


1  1  0 


S* ' TEST  ST AT.  7 ALUS  ' 

M*M'  7 ROB  [0  i  0] i  FOR  SAMPLE 

J^^QUAL  TO  «.(2  0  •«*),*.  ' 

05:5*194 

Di7C(1Spl97)*«7*61p(8pl97)«*r3 

B+  1  18  pS 
♦(i?*0  ;/I6 

F^*5C<1  ZZ9  >•  QPMT(1  7  o KK) 

P+-  1  5  4  pP 

Iff-* 

I6s?S*oO 


j  FOR  SAMPLE  SIZES  N  EQUAL  TO  ',(2  0  *NN ) . '  AND 


FS*FS,B 
C*  0 

09:C*C+1 


9:C*0+1 , 

rF*s*2y7  5  > 

*(C<7)/09 
7 :M*H, tlJ  FS 


*(C<  7 

L7;M*H 

0*0+1 

X+l 


*<«1  .£2  ,£3  ,W4  *5  ,J96  ,N7  ,0 5  )  CO] 
ffl:S+'0Xilcr  C.D.F.  ' 


,7*198 

*03 

N7:S*' ERROR;  NORMAL  ' 

PP-ERRZZ 
♦03 

NS :S~< ERROR;  NORM.  W/CC ' 
PP-ERRZZC 
*03 

NU:S+< ERROR;  T  DIST  ' 

PP+ERRTT 
*03 

NS :S*' ERROR;  T  W/CC  ' 

PP*ERRTC 
*03 

ff6:S*’£PP0P;  A7P  T/Z  ' 

PP-ERRAVE 
♦03 

N7:S*< ERROR;  AVB  TC/ZC  ' 
PP+ERRTCZC 
♦03 
04  :N 

*(iW<9  )/01 
V 


NORM.  W/CC ' 


T  DIST 


AVB  T/Z 


v  sxOTrssr;jy;AOPMi3M;<:ji:;P;Z8ZZsttH?zZj#;zzc;0;PP5X;P;P;PSsS;.r;£ja?zzc;z<; 

«  TPXS  PROGRAM  GENERATES  TABLES  OF  C.D.F.  COMPARISONS  FOR  TEE  SIGN  TEST. 

«  SPT  SAMPLE  SIZE  AND  ALPEA  VALUES. 

N*  23 

AOPffA*  0.01  0.03  0.06  0.12  0.22  0.35  0.5 
«  „„  X8XS  O0OP  INCREMENTS  SAMPLE  SIZE. 

01:PP*p0 

H*p0 

AX*p0 

ZZ*p0 

ZZC*p0 

W*W+1 

0*0 

fl*0 

«  COMPUTE  CUMULATIVE  DIST  AND  ASSOCIATED  STATS. 

P+BINOM  N 

*  XffXS  0O0P  OAOCOOArffS  ALPHA  VALUES  AND  APPORXIMATIONS . 

02:8*5+1 

0+AOPSACS] 

-i*\+/  (PSC:  ) 

X**  " 

KK*KK,K 
PP+PP  PCA] 

«  ,  COMPUTE  NORMAL  APPROXIMATION. 

Z*(.K-(q.S*N))*(0.5x(N*Q.S)) 

22+22  iNORMCDF  Z) 

«  .  ’  „  COMPUTE  NORMAL  APPROXIMATION  WITS  CONTINUITI  CORRECTION. 

« _  COMPUTE  ERROR  DIFFERENCES. 

ERRZZ+PP-ZZ 

ERRZZC*PP-ZZC 

•  .  PPXOT  OOT  TABLE  OF  VALUES. 

X*  0.1 


PC/XO  NORMAL  APPROXIMATION  WITB  CONTINUITI  CORRECTION. 
(0 . 5x0 ;  )♦ (0 . 5* (0*0 . 5 ) ; 


s+*risz  ST  AT. 
M*M,  ' 


VALUE 


0A7C (18pl97),J,61p (8pl97 },J3 


•PPCC3  ) 


♦03 

LS:J*l<an 

L3}M+M.Zn  QAF[ (18pl97  }  , 
♦(£=4)/£4_ 

♦a>*o^/£§s 

F+'SC<1  Z29  >'  OPMZCl 

F+  1  82  pP 
FS+S.F 
♦£7 

£6 :PS+p0 
S+17pS 
FS+FS.B 
C+O 

£9 :OC+l 

P-fJ.  (7  5  -PPW3) 

♦<C<75/£9 

£7:M+M,[1]  PS 
£♦0+1 

♦Tifl  ,i72  ,JV3  ,£S  )  C03 
ffl:S+*SXicZ  C.D.F. 

J*  198 
♦£3 

N2:S*< ERROR:  NORMAL 
PP+-ERRZZ 
♦£  3 

N3:S*< ERROR:  NORM.  W/CC 
PR+ERRZZC 
♦£3 
£4  :M 

♦UV<26)/£1 
7 


^POSCF  S  JP:  FOR  SAMPLE  SIZE  EQUAL  TO  ',(2  0  •*), 


7  SPMZPSZ :  J7 ;  4£?.M  ;  S :  >1  s  C ;  AA  :  ? :  Z ;  Z  2 :  PPF  Z  Z : « :  Z  ZC ;  D ; 
ZCj SP/M  PS  s  ERRTT \ ERRiC ; ER&ZiC \ZC\PC\ RBO : KK 


}PP;I;B;P;FSiS;J;TT;T;TTC; 


n  TBIS  PROGRAM  GENERATES  TABLES  OF  C.D.F.  COMPARISONS  FOR  SPEARMAN'S  R. 

«  SET  SAMPLE  SIZE  AND  ALPBA  VALUES. 

ALPBA*  0.01  0.02  0.03  0.05  0.08  0.13  0.18 
5  TBIS  LOOP  INCREMENTS  SAMPLE  SIZE. 

£l:PP-+pO 

M+pO 

FJ&pO 

FS+pO 

ZZ+pO 

ZZC+pO 

zz*pO 

ZZC+pO 
N*N+ 1 
£♦0 
B+0 

PP+6*Wx-i+ff*2 

a  COMPUTE  CUMULATIVE  VIST  AND  ASSOCIATED  STATS. 

Pi:P*PHATSPtN-S ::1 


£2 : B*B* 1 

lit®’ 

KK+KK.R&O 
PP+PP.PC2 : 
a  .  C 


CALCULATES  ALPBA  VALUES  AND  APPORXIMATIONS . 


PP+PP , PC  2 : X ] 

,,  COMPUTE  NORMAL  APPROXIMATION. 
Z*RBO*  (  (J7-1  )*0 . 5  ) 

ZZ+ZZ , ll-NORMCDF  Z ) 


a  COMPUTE  NORMAL  APPROXIMATION  WITB  CONTINUITI  CORRECTION. 

ZC*(RBO-AA  )x( (fl-i  )»o  .  3  ) 

ZZC+ZZC , l 1 -NCRMCDF  ZC ) 
a  CONFUTE  STUDENT  T  APPROXIMATION 

TT*TT.N  SPAPROX1  RBO 

a  COMPUTE  STUDENT  T  APPROXIMATION  WITB  CONTINUITI  CORRECTION. 


T^C+T^C^N  SPAPROX  RBO 


a  COMPUTE  ERROR  DIFFERENCES. 

ERRZZ*PP~ZZ 

ERRZZC+PP-ZZC 

fRRTT*PP-TT 
RRTC*PP-TTC 

a  PRINT  OUT  TABLE  OF  VALUES. 

I*  0.1 

,7+195 

S* ' TEST  ST AT.  VALUE  ' 

M*M,'  PROBIR  >  FOR  SAMPLE  SIZE  EQUAL  TO  ',(2  0  9N),< 


52  £5:.T-*-194 

53  Z3:H*H.lll  QdVrC(18pl97),^,61p(8pl97),^] 

.54,  ~(D=6 )/£4 

55  S*  1  M  55 

56  *(D* 0)/L6 

;S7i  {♦'X4C<|  9.9999  >'  DW(l  7  pU) 

,58,  ?♦  1  62  pP 

59  PS*B,P 

60  *1 7 

61  £6 : FS+o 0 

62  S*17pS 

63  PS*PS,B 

64  O  0 

65  £9:OC+l  , 

66  £♦'  T  1  .  (7  5  *PPCC]) 

67  ps*ps.p 

68  ♦CC< 7j/£9 

69  £7:#*#,tl]  PS 

70  ><-0+1 

71 3  !♦! 

|72;  --<jVl.JV2.fl3,JV4JV5,£S)[in 
73  N\iS*' EXACT  C.b.P.  1 

.74  .7*198 

[75  *£3 

[76;  N2:S*'ERR0R}  NORMAL  ' 

[77  PP*ERRZZ 
[78  *£3 

[79  N3-.S*’ ERROR}  NORM .  V/CC' 

[80  PP-ERJRZZC 
[81  *£ 3 

[ 8 2 2  J74:5*'5H50fi;  »  OPST  ' 

[83  PP+ERRTT 

84  -13 

[ 8 5 [  NilS*' ERROR}  T  W/CC  ' 

86  PP-SRRTC 
[87  *L3 

88  £4 :M 

[89 j  +(N<10 )/£l 
7 


TBIS^PROCRAM^ENERATES  TABLES  Of  C.B.P.  COMPARISONS  FOR  TBE  WILCOXON 


AZPBA*  0.0 

£1 :PP*pO 
M*p0 
KK*  pO 
PS*pO 
Z2*p0 

2zc*po 

MVpO 

rrc-po 

N+N+l 
D*  0 
B*  0 

B  C 

P+WILP  N 


SET  SAMPLE  SIZE  AND  ALPS A  VALDES. 

01  0.02  0.03  0.05  0.08  0.13  0.18_ 
TBIS  LOOP  INCREMENTS  SAMPLE  SIZE. 


COMPOTE  COMOLATIVE  BIST  AND  ASSOCIATED  STATS. 

TBIS  LOOP  CALCULATES  ALPBA  VALVES  AND  APPORXIMATIONS . 


L2:B*B*1 

r-i-i 

KK-KK.T 

PP*PP.PZA1 

a  COMPUTE  NORMAL  APPROXIMATION . 

WOT**r-((N*(«*l)  )  +  4) 

DEN-*-  (  (W*  {N*2  jx  (  (  2*/V)+l )  )*2U  )*0 . 5 

Z~NUM*DEN 

ZZ+ZZ .SORMCDP  Z 

a  COMPUTE  NORMAL  APPROXIMATION  WITB  CONTINUITI  CORRECTION. 

«a»e*(r+o. 5)-((w*<«+i5 >4) 

ZC+NVMC*DEN 
ZZC*ZZC .NORMCDF  ZC 

a  ,  COMPUTE  STUDENT  T  APPROXIMATION  , 

DENOM*((lN>‘  (DEN+2j)*Xn-1)  )-{(NUH*  2  )♦  CiV-1 ) )  )*0. 5 

T+NUM*DENOM 

TT*TT.(N-lj  TDISTN  T 

a  COMPUTE  STUDENT  T  APPROXIMATION  WITB  CONTINUITI  CORRECTION. 
NOHT*XnUM)+0 .5 

( (W*  (CffW*2  )  )♦  (P-1 )  )- ( ( ( (WWn  +  0 . 5  )*2  )♦  (#-l  )  )  )*0 . 5 

Tc+numt*denomc 

TTC+TTC. (Xn-1)  TDISTN  TC) 


;(AW#*2)*(S-1)))*0.5 


♦(P<7)/£2 

ft  COMPUTE  ERROR  DIFFERENCES • 

MPZZ*PP-Z2 
\RRZtc+PP -ZZC 

_  . 


|Sfi^S*PP-ffzZ+IT)*2 ) 

PJMTCZOPP-  ( <ZZC+ZTC)*2  ) 


55]  •  PRINT  00T  TABLE  OF  VALDES. 

56]  JV0.1 
57 j  «/»19  5 

58]  S*'TEST  ST AT.  VALUE  ' 

59]  M*M, 1  PROBIN  S  JT]j  FOR  SAMPLE  SIZE  EQUAL  TO  ',(2  0  *N) 

60]  +£3 

61]  £5:.7*194 

62]  Z3:H*#,CT]  QA7C(18pl97),«7t61p(8pl97)t,7] 

63]  ♦(P=8)7iu 
6«]  B *  1  18  oS 

65]  *(D*0)/Z6 

66]  P*'BC< I  ZZ9  >'  QPMTCl  7  p/t/t) 

67]  F *  1  62  pP 

68]  ?S*B.P 

69]  *L 7 

70 3  L6:PS*pO 

T 1 ]  B*l7oS 

72]  FS-FS.B 

73]  .00 
7u]  £9:001 

75]  ?*'  I' 1.(7  5  9PPZC2) 

-6]  PS*PSJ 

77]  *(C<75/£9 

78]  L7 :M*H, ZlZ  PS 
79.  OOl 

80  Ol 

81,  *IN  1 . P2  ,P3  ,P4  , /VS  ,P6  ,P7  ,15  )  CP] 

82  Pl:5*rPXicr  C.b.F.  f 

83,  0198 

84,  *13 

85]  N2:S*' ERROR;  NORMAL  ' 

86]  PP*ERRZZ 

|  '"K' 

|  *  VZSI  ' 

|  iUosn0”  r  */ec  • 

Hi  iikm“  m  r/z  ■ 

>9!  ,  *L 3 

100]  N7 :S*' ERROR;  AVE  TC/ZC  ' 

101]  PP*ERRTCZC 
L02]  *L 3 


103]  £4:# 


»(P<10)/£1 


niumnwMWJumwwwwpwwwMM 
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